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PHYTOMICS: A GENOMIC-BASED APPROACH TO HERBAL COMPOSITIONS 



This invention relates to herbal compositions. More specifically, this invention 
provides tools and methodologies for improving the selection, testing, quality control and 
manufacture of herbal compositions, and to help guide the development of new herbal 
compositions and identify novel uses of existing herbal compositions. 

BACKGROUND OF THE INVENTION 

All publications and patent applications herein are incorporated by reference to the 
same extent as if each individual publication or patent application was specifically and 
individually indicated to be incorporated by reference. 

Herbal medicine has been in use for centuries by people of Asia and Europe. In the 
United States (US), herbs have become commercially valuable in the dietary supplement 
industry as well as in hohstic medicine. Approximately one third of the US population has 
tried some form of alternative i medicine at least once (Eisenberg et al., 1993, N. Engl. J. Med. 
328:246-252). 

Botanicals, including herbs,- have also become a focal point for the identification of 
new active agents to treat diseases. Active compounds, derived from plant extracts, are of 
continuing interest to me 'pharmaceutical industry. For example, taxol is an antineoplastic 
drug obtained from the bark of the western yew tree. It is estimated that approximately 50 
percent of the thousands of drugs commonly used and prescribed today are either derived 
from a plant source of Contain chemical imitations of a plant compound (Mindell, E.R., 1992, 
Earl Mindell's Herb Bible, A Fireside Book). 

Chirrently, a number of medicinal formulations, food supplements, dietary 
supplements and the like contain herbal components or extracts from herbs. Herbal 
medicines have been used for treating various diseases of humans and animals in many 
different countries for a very long period of time (see, e.g., LA. Ross, 1999, Medicinal Plants 
of the World, Chemical Constituents, traditional and Modern Medicinal Uses, 
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Humana Press; D. Mqlqny,,1998, The American Association of Oriental Medicine's Complete 
Guide to Chinese Herbal Medicine, Berkley Books; Kessler et al., 1996, The Doctor's 
Complete Guide to Healing Medicines, Berkley Health/Reference Books); Mindell, supra). 

Herbal Medicines. There are many branches of herbal medicine around the world, 
such as Ayurveda, Unani, Sida and Traditional Chinese Medicine (TCM). While modern 
Western medicine typically consists of adniinistering a single chemical entity capable of 
intervening specific biochemical pathways, each formula of TCM typically contains hundreds 
of chemical entities from several herbs which are designed to interact with multiple targets in 
the body in a coordinated manner. . Although empirical practice contributed in a significant 
way to the herbal composition and prescription of these ancient herbal medicines, they are also 
supported, to a varying degree, by a set of theories which all are distinct from that of modern 
Western medicine in terms of | anatomy, pharmacology, pathology, diagnosis treatment, etc. 
Among the different herbal medicine fields, TCM has developed a more complete set of 
theories over several cenmries which have been well documented and practiced by local 
physicians caring for a huge population (>1.3 billion people) in greater China and in East Asia 
including Korea and Japan; 

Western medicine generally uses purified compounds, either natural or synthetic, 
mostly directed towards a single physiological target. However, the compositions used in 
TCM are usually composed of multiple herbs and compounds which are aimed at multiple 
targets in the body based on' unique and holistic concepts. TCM mainly used processed crude 
natural products, with various combinations and formulations, to treat different conformations 
resulting in fewer side eff ects. The great potential of TCM has yet to be realized for the 
majority of the world's people,;: ' : ' ' 

The herbs in a typical TCM prescription are assigned roles as the principal herb and the 
secondary herbs, mcluding assistant, adjuvant and guiding herbs. The principal herb produces 
the leading effects in freatihg me cause or the main symptom of a disease. An assistant herb 
helps to strengthen the effect of the principal herb and produces leading effects in the treatment 
of the accompanying symptoms. There are three types of adjuvant herbs: 1) those that enhance 
the therapeutic effects of the principal and assistant herbs or treat tertiary symptoms, 2) those 
that reduce or ehminate uie toxicity and other side effects of the principal and the assistant 
herbs and 3) those which act bh complementary target tissues not specifically affected by the 
principal herb. A guiding herb directs the effect of other herbs to the affected site and/or 
coordinates and mediates the effects of the other herbs in the prescription or formulation. In 
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contrast to most of the herbal medicines or supplements that consist of one or more parts of a 
single plant, the intended effects of TCM. are directed at multiple tissues. 

For example, a well-known TCM recipe, "Ephedra Decoction" used for treating asthma 
is composed of ephedra, cinnamon twig, bitter apricot kernel and licorice. Ephedra, is the 
principal herb, which expels cold, induces diaphoresis and facilitates the flow of the Lung Qi 
to relieve asthma, the main symptom. Cinnamon twig, as the assistant herb, enhances 
ephedra's induction of diaphoresis and warms the Channels to ensure the flow of Yang Qi for 
redueing headache and pantalgia. Bitter apricot kernel, as the adjuvant herb, facilitates the 
adverse flow of the Lung Qi and strengthens the asthma relief by ephedra. Licorice as the 
guiding herb moderates the effects of both ephedra and cinnamon to ensure a homeostasis of 
the vital Qi. While each of the four herbs clearly exhibits its respective activity, they 
complement as well as supplement each other when they are combined. In practice, the 
principal herb can be prescribed with one or more secondary herbs, depending on the 
symptoms at a patient's presentation (Prescriptions of Traditional Chinese Medicine, Chapter 
One, pplO-16, E. Zhang, editpr in Chief Publishing House, Shanghai University of Traditional 
Chinese Medicine, 1998). "... : / .. - : v. 

The main theories, of TCM that guide the treatment of sickness with herbal medicine 
and other means, such as acupuncture, are 1) the theory of Yin and Yang, 2) the theory of Five 
Elements, 3) the theory of Viscera and Bowels, 4) the theory of Qi, Blood and Body Fluid, and 
5) the theory of Channels and Collaterals. 

In TCM, the first important aspect of making the proper diagnosis is to ascertain 
whether the disease is Yin of Yang. For example, those patients who have a fever, are thirsty, 
constipated or have a rapid pulse condition are of Yang character. Those individuals who have 
an aversion to cold, are not thirsty, and diarrhea and a slow pulse condition are of Yin 
character. The property, flavor and function of herbs can also be classified according to Ying 
and Yang theory. For example, herbs of cold and cool nature belong to Ying, while herbs 
which are warm and hot in nature belong to Yang. Herbs with sour, bitter and salty flavor 
belong to Ying, while herbs with pungent, sweet and bland flavor belong to Yang. Herbs with 
astringent and subsiding function belong to Yin, while herbs with dispersing, ascending and 
floating function belong to Yang. In TCM, the principles of treatment are based on the 
predominance or weakness of Yin and Yang. Herbs are prescribed according to their property 
of Ying and Yang and their function for restoring the imbalance of the Ying and Yang. In so 
doing, the benefit of treatment is achieved. 
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According to the theory of Five Elements there are five basic substances that constitute 
the material world (i.e., wood^ fire, earth, metal and water). In TCM, this theory has been used 
to explain the physiology and pathology of the human body and to guide clinical diagnosis and 
treatment. Herbal physicians, have applied the laws of generation, restriction, subjugation and 
reverse restriction of the five elements to work; put many effective and specific treatment 
regimens, such as reinforcing earth to generate metal (strengthening the function of the spleen 
to benefit the lung), replenishing water to nourish wood (nourishing the essence of the kidney 
to benefit the liver), supporting, earth to restrict the wood (supplementing the function of the 
spleen to treat the hyperactivity of the liver), and strengthening water to control fire 
(replenishing the essence of the kidney to treat hyperactivity of the heart). Specifically, the 
property of some herbs is assigned to each of the five Elements for the purposes of guiding the 
prescription of a TCM recipe. 

In TCM, the internal organs of the human body are divided into three groups: five 
Viscera (the Heart, the Liver;: the Spleen, the Lung and the Kidney), Six Bowels (the Gall 
Bladder, the Stomach, the.'large Intestine, the Small Intestine, the Urinary Bladder, and the 
Triple Warmer), the Extraordinary Organs (the Brain, the Medulla, the Bone, the Blood 
Vessel, the Gall Bladder, arid the Uterus). In TCM, the Viscera or the Bowel are not only 
anatomic units, but are also concepts of physiology and pathology about interactions between 
different organs. For example, the heart also refers to some of the mental functions and 
influence functions of blood, hair, tongue and skin. Ying-Yang and the Five Elements 
influence the interactions among' these Viscera, Bowels and Organs. The complexity of 
interplay of the theories is used Uo explain the pathology of diseases to which herbs are 
prescribed, as discussed below.: 

The prescription of herbal medicine in TCM starts with the diagnosis, which consists of 
four main items: mterrogatioh/iiispection, auscultation and olfaction, pulse taking and 
palpation. During the interrogation phase, much information is gathered, including the 
characteristics of the main symptoms. For instance, if the main symptom is characterized by 
dull pain of epigastric region,' which may be relieved by wanning and pressing, this suggests 
the insufficiency of the Spleen- Yang. Soreness and weakness of the loins and knees, 
intolerance of coldness with cold extremities manifests a weakness of the Kidney-Yang. 
During inspection, observations are made for vitality, skin color and the general appearance 
and the condition of the tongue. For example, a pale complexion corresponds internally to the 
Lung of autumn, whose Qi is i dry. This may occur when Yang Qi is lacking and the circulation 
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of Qi and blood is impeded, or when the coldness in the channels and collaterals causes them 
to contract 

in TCM, it is from Qi, blood and body fluid that come energy needed by the Viscera 
and Bowels, Channels and Collaterals, tissues and other organs for carrying-out their 
5 physiological functions; and on which the formation and metabolism of Qi, blood and body 
fluid depend. Prescriptions of TCM consider the herbal effects on Qi and blood for treatments. 

TCM holds that Channels, Collaterals and their subsidiary parts are distributed over the 
entire body. It is through them that herbs exert influence on pathological targets and achieve 
the improvement of sickness. For example, ephedra acts on the Channels of the Lung and 

10 Urinary Bladder so as to induce sweat for reheving asthma and promoting diuresis. As noted 
above, clinical applications of acupuncture are also guided by the theory of Channels and 
Collaterals. '/ .'-. 

In summary, while the nature or property of each herb in TCM may be assigned as Yin 
or Yang, and to one of the Five Elements, they act through Channels and Collaterals and are 

15 mediated via Qi, Blood and Fluid to yield therapeutic effects on targets, such as Viscera and 
Bowels. Pathogenic factors may bedisguised as decoy through the very same systems of 
Channels and Collaterals to adversely affect the functions of Viscera and Bowels and thus 
cause sickness. . • 

From the foregoing discussion, it is clear that the TCM terminology is as much of a 

20 philosophical concept as an anatomical one. For example, the Heart represents a host of 

tissues, organs or systems in the body that contribute to a function described in TCM. Thus, 
the concept of the Heart requires a multiple dimension data set to describe each concept of 
TCM. Once this is accomplished, a molecular holistic medicine can be developed. 

U.S. Regulatory J?r6cess. in the US, dietary supplements (such as botanical products, 

25 vitamins and minerals, amino acids and tissue extracts) are regulated under the Dietary 
Supplement Health and Education Act of 1994 (the DSHE Act). This Act removed the 
ingredients of dietary supplements from regulation as food additives under the Federal Food, 
Drug, and Cosmetic Act: In addition, the DSHE Act requires that The Food and Drug 
Admimstration (FDA) bear the burden of proof that a marketed dietary supplement presents a 

30 serious or unreasonable risk under the conditions of use on the label or as commonly 

consumed. Thus, there are currently no federal regulations that establish specific criteria for 
purity, identification and manufacturing procedures for dietary supplements. In addition, few 
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published papers on herbal quality have resulted from the establishment of the Office of 
Alternative Medicine by Congress in 1992 (Angell et al, 1998, N. Engl. T Med . 339:839-841). 

At the present time, the FDA must approve each one of the chemical entities in a drug 
composition or cocktail, and then clinical trials must be undertaken so as to obtain separate 
5 FDA approval for marketing the drug. This process is extremely tedious and costly. A 

molecular holistic mecucine may require a less arduous evaluation since the previous use of a 
particular herbal composition as a.botanical drug permits clinical trials with multiple chemicals 
at the outset {i.e., clinical trials using the herbal composition or specific components of the 
herbal composition). Recently, the FDA has approved the testing of some herbal medicines in 

10 clinical trials as botanical drugs (FDA Guidance on Botanical Drugs, April, 1997). While 

these events represent a positive development for health care in general, it also raises important 
issues regarding the formulation, manufacturing and quality control of herbal medicines and 
dietary supplements, mcluding the. traditional Chinese medicines. 

The multitude of relevant biological responses induced by the multiple chemicals in 

1 5 herbs are not currently available and will be increasingly important to support marketing 
approval by the FDA Y~:. 

Herbal-based industries are coming under increasing pressure to upgrade their current 
practices (see, e.g., Angell et al., supra). The need to apply scientific testing to the preparation 
and administration of herbal medicines , and food supplements has been highlighted by several 

20 recent reports of toxicity resulting from ingesting herb-based formulations. For example, one 
patient who took an herbal 7 based dietary supplement experienced digitalis toxicity (Slifinan et 
al., 1998, N. Engl. J. Med. 339:806-81 1). It was subsequently determined that the herb 
ingredient labeled as plantam in the supplement was actually contaminated with Digitalis 
lanata, an herb known to contain at! least 60 cardiac glycosides. In another instance, an herbal 

25 preparation was found to be the cause of chronic lead intoxication in a patient (Beigel et al., 
1998, N. Engl. J. Med. 339:827-830}. This is not a completely unexpected occurrence since 
contamination of traditional Asian herbal remedies by lead and other heavy metals is well 
documented (Woolf et at, 1994, Ann. Intern. Med. 121:729-735). 

Characterization of Botanicals. It is well known that the genetic identity (e.g., 

30 genera, species, cultivar, variety, clone), age of herbal growth, harvest time, the specific plant 
part utilized, processing method,'ge6graphical origin, soil type, weather patterns, type and rate 
of fertilizer, and other growth factors have a great impact on the particular chemical 
composition of any particular herb "harvested" from any particular area. 



WO 01/66803 PCT/US01/07608 

-7- 

Increasing numbers of various types of tests have been instituted to assure the 
consistent quality of herbs used in'medicine and as dietary supplements; including inspections 
at the macro- and microscopic levels as well as a variety of chemical analyses. Recently, high 
performance liquid chromatography (HPLC) profile of marker molecules in an herbal extract 
has become one reference standard. However, there are problems with this approach, 
including that some of the bioactive molecules may not adsorb UV or the visible lights for 
HPLC detection, and the amount of a chemical is not necessarily proportional to its biological 
potency. For these reasons, herbal manufacturers resort to a practice of mixing raw herbs from 
different sources to mini mize chemical variations. 

Mass spectrometry (MS) is an analytical method for determining the relative masses 
and relative abundances of components of a beam of ionized molecules or molecular fragments 
produced from a sample in a high vacuum. MS, unlike HPLC, is not optical density- 
dependent. In practice it is i^ed in conjunction with HPLC or capillary electrophoresis (CE): 
the HPLC separates the chemicals and the MS then can be used to identify what they are. 
Commercial systems are available which integrate MS and HPLC for biological uses. Mass 
spectrometry is limited to samples that are gaseous or volatile at low pressure, or that can be so 
rendered by derivatization. . 

These steps are no longer adequate. Recent publications report a greater variation in 
the quality of herbs by specific suppliers, and the difficulty of providing biological equivalence 
of herbal extracts. Furthermore, the correlation between safety and efficacy and chemicals in 
an herb is not well defined in most cases. Recently, in response to complaints from consumer 
groups and regulatory agencies (Federal Register, February 6, 1997, Volume 62, No. 25, 
Docket No. 96M-0417, cGIMDP'in Manufacturing, Packing or Holding Dietary Supplements, 
Proposed Rules), some herbal manufacturers have begun to implement Good Manufacturing 
Practice (GMP) which requires stringent controls at all levels. 

Chemical and spectroscopic methods have been used to characterize the components of 
herbal medicines and food supplements. For example, three new hederagenin-based acetylated 
saponins were isolated from the fruits of Gliricidia sepium using these two methods (Kojima et 

1998, Phytochemistry 48(5):885-888). The botanical sources of Chinese herbal drugs in a 
number of commercial samples were inferred by comparing the contents of some characteristic 
constituents which were analyzed With high-performance chromatography (HPLC) or capillary 
electrophoresis (CE) (Shuenn-Jyi Sheu, 1997, Journal of Food and Dru g Analysis 5(4):285- 
294). For example, the ratio of ephedrine/pseudoephedrine was used as a marker to 
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differentiate Ephedra intermedia from other species; total alkaloid contents were used to 
distinguish between species of Phellodendron; and the contents of ginsenosides were used to 
differentiate between species of Panax. However, these methods do not provide a direct 
measurement of the effect of the various herbs on the molecular, physiological or 
morphological responses following human treatment with the herbs. 

Using gas cbxomatography-mass spectrometry and atomic-absorption methods, the 
California Department of Health Sciences, Food and Drug Branch, recently tested Asian 
medicines obtained from . herbal stores for contaminants (R. J. Ko, 1998, N. Eng l. J. MeH 
339:847). Of the 260 products they tested, at least 83 (32 percent) contained undeclared 
pharmaceuticals or heavy metals, and 23 had more than one adulterant. Using high- 
performance liquid chromatography, gas chromatography, and mass spectrometry, a 
commercially available combination of eight herbs (PC-SPES), was found to contain 
estrogenic organic compounds (DiPaola et al, 1998, N. Engl. J. Med. 339:785-791). The 
researchers concluded that PC-SPES has potent estrogenic activity and that prostate cancer 
patients that took PC-SPES could confound the results of standard therapies and may 
experience clinically significant adverse effects. Gas chromatography data was also collected 
for different samples of the traditional Chinese medicine 'wei ling xian' and correlated to the 
antiinflammatory activity of the samples (Wei et al., Study of chemical pattern recognition as 
applied to quality assessment of the traditional Chinese medicine "wei ling xian," Yao Hsueh 
Pao 26(10): 772-772 (1991)): This study did not provide relevant HBR Array data, such as 
time course, dose dependent response, control samples to substantiate the differential power of 
the biomarkers, nor it utilize a reiterative type of data construction process to establish a 
comprehensive database for characterizing effects of the herbal composition. 

Changes in protein levels have also been used to characterize the effects of herbal 
compositions or specific components of herbs. , For example, the production of granulocyte 
colony-stimulating factor (G-CSF) from peripheral blood mononuclear cells was found to vary 
depending on which specific Chinese herb was added to the culture (Yamashiki et al., 1992, L 
Clin. Lab. Immunol. 37(2):83-90). Expression of interleukin-1 alpha receptors was markedly 
up regulated in cultured human epidermal keratinocytes treated with Sho-saiko-to, the most 
commonly used herbal medicine in Japan (Matsumoto et al, 1997, Jpn. J. Pharmar.nl 
73(4):333-336). The expression of Fc gamma 11/111 receptors and complement receptor 3 of 
macrophages were increased by treatment with Toki-shakuyakusan (TSS) (J. C. Cyong, 1997, 
Nippon Yakurjgaku Zasshi 1 10(Suppl.. l):87-92). Tetrandrine, an alkaloid isolated from a 
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natural Chinese herbal medicine, inhibited signal-induced NF-kappa B activation in rat 
alveolar macrophages (Chen et dl., ,'1997, Biochem. Biophvs. Res. Commim . 231(1):99-102). 
The herbs Sairei-to, alismatis rhizoma (Japanese name "Takusha") and hoelen (Japanese name 
"Bukuryou") inhibited the synthesis and expression of endothelin-1 in rats with anti-glomerular 
5 basement membrane nephritis (Hattori et al., 1997, Nippon Jinzo Gakkai Shi 39(2):121-128). 

The increase or decrease m mRNA levels has also been used as an indicator of the 
effect of various herbs and herbal components. Intraperitoneal injection of Qingyangshen 
(QYS), a traditional Chinese medicine with antiepileptic properties, and diphenylhydantoin 
sodium reduced alpha- and beta-tublin mRNAs and bippocampal c-fos mRNA induction 

10 during kainic acid-induced chronic seizures in rats (Guo et al., 1993, J. Tradit. Chin. Merl 
13(4):281-286; Guo et al. T 1995. J. Tradit. Chin. Med . 15(4):292-296; Guo etal., 1996, L 
Tradit. Chin. Med . 16(l):48r51): Treatment of cultured human umbilical vein endothelial cells 
(HUVECs) with the saponin. astragaloside IV, a component purified from Astragalus 
membranaceus, decreased plasminogen activator inhibitor type I (PAI-1) specific mRNA 

15 expression and increased tissue-type plasminogen activator (t-PA) specific mRNA (Zhang et 
al., 1997, J. Vase. Res . 34(4):273-280). One component isolated from the root of Panax 
ginseng was found to be a potent inducer of interleukin-8 (IL-8) production by human 
monocytes and by the human monocytic cell line THP-1, with this induction being 
accompanied by increased EL-8 mRNA expression (Sonoda et al., 1998, Immunopharmacology 

20 38:287-294). ■ . 

Recent advances in nucleic acid microarray technology enable massive parallel mining 
of information on gene expression; This process has been used to study cell cycles, 
biochemical pathways, genome-wide expression in yeast, cell growth, cellular differentiation, 
cellular responses to a single chemical. compound, and genetic diseases, including the onset 

25 and progression of the diseases (M. Schena et al., 1998. TIBTECH 16:301). Because cells 

respond to the nticro^envifohment changes by changing the expression level of specific genes, 
the identities of genes expressed in a cell determine what the cell is derived of and what 
biochemical and regulatory systems are involved, among other things (Brown et al. 1999, 
Nature genet , 21 (1) supplemeht:33). Thus, cellular gene expression profiles portray the 

30 origin, the present differentiation of the cell, and the cellular responses to external stimulants. 
No researchers to date, if any, have attempted to apply these new technologies to study the 
molecular effects of whole herbal treatments and supplements. 
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Some researchers have, attempted to characterize the effects of the major active 
constituents isolated from selected herbs. For example, treatment of HUVECs with 
notoginsenoside Rl (NR1), purified from Panax notoginseng, resulted in a dose- and time- 
dependent increase in TPA synthesis (Zhang et al., 1994, Arteriosclerosis and Thromobosis 
5 14(7):1040-1046). Treatment with NR1 did not change urokinase-type plasminogen activator 
and PAI-1 antigen synthesis, nor did it effect the deposition of PAI-1 in the extracellular 
matrix. TP A mRNA increased as much as twofold when HUVECs were treated with NR1, 
whereas expression of PAI-1 -specific mRNA was not significantly affected by NR1. Since 
most studies oh P. notoginseng have involved its mixture with other herbs, the researchers 

10 noted that it was difficult to assess how their results relate to the situation in vivo when is used 
therapeutically in humans (M., at 1045, second column, first paragraph). In addition, since the 
researchers only studied one major componenfof the herb, it is not possible to ascertain the 
molecular effect of the whole herb or the interactions among components of the herb from this 
study. 1 

15 Dobashi et al. (1 995: Nteuroscience Letters 197:235-238) studied the effect of two of 

the main components of saiko agents, a Chinese herbal drug used to treat nephrotic syndrome, 
bronchial asthma and chronic rheumatoid arthritis. Adrninistration of SS-d increased plasma 
adrenocorticotropin (ACTH) levels; proopiomelanocortin mRNA levels in the anterior 
pituitary and the CRT mRNA level in the rat hypothalamus in a dose dependent manner. In 

20 contrast, treatment with SS-a failed to affect the levels of these molecular markers. While this 
study indicates that admimstratipn of SS-d may have an important role in saiko agents-induced 
CRF release and CRF gene expression in rat hypothalamus, it fails to address the molecular 
effect of the herbal medication as a whole. 

Kojima et al (1998! Biol. Pharm. Bull. 4:426-428) describe the utilization of 

25 differential display of mRNA to isolate and identify genes transcriptionally regulated in mouse 
liver by sho-saiko-to, an herbai medicine used for treating various inflammatory diseases in 
Japan. These researchers limited their study to the use of mRNA differential display 
techniques in investigating the molecular mechanisms of herbal medicine. It also failed to 
address effects in multiple organs 6f treated animals and did not provide any guidance for 

30 quality control, new use, and standardization of effects. In addition, the study failed to analyze 
the individual components of the herb and compare the individual results with the results 
obtained using the whole herbal mixture. 
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Ma Ji et al. (1 998; Chinese Medical Journal 1 1 1 (1): 1 7-23) investigated the therapeutic 
effect of the herb Astragali membranacetis on sodium and water retention in rats experiencing 
aortocaval fistula-caused experimental congestive heart failure. Chronic heart failure rats with 
and without Astroglia treatment were compared for changes in various morphological 
5 characteristics (e.g., body weight, serum sodium concentration); physiological characteristics 
(e.g., mean arterial pressure, heart rate, hematocrit and plasma osmolality); mRNA expression 
levels (e.g., hypothalamic, arginine vasopressin (AVP), AVP Via receptor, renal AVP V 2 
receptor, aquaporin-2 (AXP2)) and protein excretion (e.g., plasma atrial monophosphate 
peptide (ANP) and urinary cyclic guanidino monophosphate (cGMP)). The researchers found 

1 0 that treatment with Astroglia improved cardiac and renal functions, partially corrected 

abnormal mRNA expressions of the AVP system and AQP2, and improved the renal reaction 
to ANP. This study did hot address using the collected data to guide the development of new 
formulations or for elucidating the synergistic or other interactions among various herbs in a 
formula, or validate the differential power of the effects for quality control purposes. 

15 As shown by the above ; review of relevant scientific articles, molecular-based 

technology has not been used to explore and validate cellular and molecular responses in 
biological systems that are treated - or challenged with, multiple chemicals at the same time, 
such as herbal medicines aridTCM. Furthermore, these recent advances have not been 
integrated with other technologies and methods to produce a process for the systematic 

20 exploration of biological effects of herbal medicines and TCM. 

SUMMARY OF THE INVENTION 
This invention provides' the tools and methodologies for creating, mamtaining, 
improving and utilizing Herbal BioResponse Arrays (HBR Arrays), wherein the HBR Arrays 
constitute data sets associated; with particular herbal compositions. The HBR Arrays of the 

25 present invention may include' information on'the plant-related parameters of the herbal 

constituents, marker information collected following the exposure of a biosystem to the herbal 
composition, and biological response information collected following the exposure of a 
biosystem to the herbal composition. 

The present invention' provides the tools and methodologies necessary for estabhshing 

30 standardized HBR Arrays for particular herbal compositions, wherein the standardized HBR 
Arrays are used as benchmarks by which to evaluate batches of similar or different herbal 
compositions. The present invention further provides the tools and methodologies necessary 
to update and maintain the standardized HBR Arrays. Particular embodiments of the present 
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invention involve iterative processes whereby data for additional batches of the herbal 
composition, additional plant-related data, additional marker information, and/or additional 
BioResponse information is periodically added to the standardized HBR Arrays. Thus, the 
present invention provides the tools and methodologies for creating, mainta inin g, updating and 
5 using HBR Arrays on an ongoing basis. 

The present invention provides the tools and methodologies necessary to guide the 
standardization of herbal .compositions; to deteraiine which specific components of herbal 
compositions are responsible for particular biological activities; to predict the biological 
activities of herbal compositions; for the development of improved herbal therapeutics; for 

10 adjusting or modifying an herbal composition; for measuring the relatedness of different herbal 
compositions; for identifying specific molecules in the batch herbal composition which retain 
the desired biological activity; for detennining which herbal components of a known herbal 
composition can be eliminated from the known herbal composition while mamtaining or 
improving the desired biological activity of the known herbal composition; for identifying new 

15 uses and previously unlmbwh biological activities for the batch herbal composition; and for 
using the predicted biological activity of the batch herbal composition to aid in the design of 
therapeutics which include herbal components and synthetic chemical drugs, including the 
design of therapeutics using ithe methods of combinatorial chemistry. 

More specifically, the present invention provides methods of establishing standardized 

20 Herbal BioResponse Arrays (I3BR Arrays) for herbal compositions, wherein the methods 
comprise: • 

a) selecting a characterized herbal composition; 

b) exposing a biosystem td a batch of the characterized herbal composition and collecting data 
on two or more markers, wherein one of the markers is a change in gene expression 

25 determined through the ,use ; df a nucleic acid microarray, produced by the steps comprising: 

i) producing a cell banking system; 

ii) profiling trie- gehe' expression pattern of cells from the cell banking system 
before and after exposure to the herbal composition; 

iii) selecting as markers those genes whose expression levels are changed by 
30 exposure to the' herbal composition; 

c) storing the marker data 6f step b) as a standardized HBR array. 

The present invention further provides such methods which further comprise exposing a 
biosystem to on£ or moreTjatcties of the herbal composition, collecting the data on one or more 
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BioResponses, and adding the collected BioResponse data to the standardized HBR Array for 
that herbal composition. . 

The present invention provides methods of evaluating herbal compositions, wherein the 
methods comprise exposing a biosystem to a batch of the herbal composition and collecting 
5 data on two or more markers; and comparing the collected marker data with a standardized 
( HBR Array for the same or a substantially same herbal composition as that of the batch herbal 
compositions. 

- The present invention provides a system for predicting the biological activity of an 
herbal composition comprising: 
10 1). a biosystem comprising one or more different types of cells, tissues, organs or 

in vitro assays; ... r\.- 1 

2) . a batch herbal composition; 

3) . two or more molecular markers; 

4) . a means for exposing the biosystem to the batch herbal composition and 
1 5 measuring the differential responses of the molecular markers; 

5) . a computer processor, including memory, for analyzing and storing the 
differential response measurements of the molecular markers so as to create an Herbal 
BioResponse Array (HBR Array) data set for the batch herbal composition; 

6) . a computed processor, including memory, for comparing the HBR Array of the 
20 batch herbal composition to one or more previously-stored HBR Arrays so as to predict the 

biological activity Of the batch, herbal composition, wherein the biological activities of the 
herbal compositions used to generate the one or more previously-stored HBR Arrays are 
known. . :'\ ' 

BRIEF DESCRIPTION GF THE DRAWINGS 

25 Figure 1 . Figure 1 provides a schematic of the basic method steps for constructing a 

Standardized Herbal BioResponse Array (HBR Array) for any selected herbal composition. 
The figure is shown in its rhost basic form for ease of understanding. As discussed herein, 
each of the pathways of the schematic can be done iteratively. Furthermore, any information 
contained in one box can be used to guide decisions regarding gathering information for any 

30 other box. In this way, numerous feedback loops are also possible throughout the scheme. 

Figure 2. Figure 2*prpviides a- schematic of the basic method steps for constructing a an 
Herbal BioResponse Array (HBR Array) for any batch herbal composition and for comparing 
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this batch HRB Array to a selected subset of information from the Standardized HBR Array. 
The figure is shown in its most basic form for ease of understanding. As discussed herein, 
each of the pathways of the schematic can be done iteratively. Furthermore, any information 
contained in one box can be used to guide decisions regarding gathering information for any 
5 other box. In this way, numerous feedback loops are possible throughout the scheme. 

Figure 3. Figure 3 provides a schematic of the basic method steps for estabushing and 
using a major data set. The figure is shown in its most basic form for ease of understanding. 
As discussed herein, each of the pathways of the schematic can be done iteratively. 
10 Furthermore, any information contained in one box can be used to guide decisions regarding 
gathering information for any other box. In this way, numerous feedback loops are possible 
throughout the scheme. 

Figure 4. Western blot for various herbal compositions. 
15 A. No herbal composition. .. 

B. Huang Qihg Tang A (HQT A) (0.2 mg/ml). 

C. HQT A (4 mg/ml). 

D. HQT B (0.2 mg/ml). 

E. HQT B (4 mg/ml). 
20 F. Scute (0.2 mg/ml). 

G. Scute (4 mg/mlj. . . 

Figure 5. HPLC for Paeonie lactiflora pallus. 

25 Figure 6. HPLC toxZiziphifructus. 

Figure 7. Figure 7 provides a schematic for estobhshing a bio-response data set for an 
herbal composition. The data set is based on differentially expressed gene induced by the 
herbal medicine for more . than three different concentrations in a mammalian cell culture. 



30 



Figure 8. Figure 8 provides a schematic for establishing a characteristic expression 
profile database or HBR Array for an herbal medicine or a complex herbal preparation. 



;il I'M! iiji »•■»••":':: : , if .i'i ii" i; im) «":!' iiji « ;J; 



WO 01/66803 PCT/US01/07608 

-15- 

Figure 9. Figure .9 provides a schematic for identifying an unknown herbal 
composition. The expression profiles induced by the unknown herbal medicine are aligned 
with the expression profile database and statistical method is employed to score the possible 
identities of herbal medicines archived in the database. 

5 '.' /'•'; [ '. ; 

Figure 10. Figure 10 provides a schematic for extracting signature genes for an herbal 
composition or a complex herbal preparation. 

Figure 11. Figure 1 1 provides a schematic for extracting signature genes for individual 
chemical constituents in ah herbal medicine or a complex herbal preparation. 

1 0 Figure 12. Clustered display of gene expression data from cells treated with three types 

of single-element herbal extracts (Cordyceps Sinensis Mycelium(CSM), ST024, ST117) with 
high and low concentrations ^(indicated with H and L, respectively). 

(A) Cluster analysis was performed by the program "Cluster" (Eisen et al., 1999) with 
492 selected genes (see text). 

1 5 (B) Enlarged image of genes up-regulated by ST 1 1 7 treatment but down-regulated by 

other herbal extract treatments. The clone ID and putative gene name are indicated. 

(C) The clustering algorithm separated CSM, ST024 and ST1 17 into 3 distinct clusters. 
The distance between each cluster as displayed by the hierarchical dendrogram can be viewed 
as the difference between'the expression profiles of the three herbal extracts treated cells. 

20 • . -' ... '.: 

Figure 13. .? - '■ 

(A) Pseudo-color encoded display of clustering results as calculated based on the 
selected 492 genes. The boxes in (A) indicate the positions of the three clusters of genes 
described above. 

25 (B) Enlarged image Of genes down-regulated by the CSM but up-regulated by the 

others. 

(C) Genes up-regulated by all kinds of herbal treatments. 

(D) Genes down-regulated by CSM and up-regulated by the others. The IMAGE clone 
ID and putative gene name are mdicated. I 

30 '• \ .:- 7 -M-. ■ 

Figure 14. Clustered display of expression data from 2 batches of multi-element herbal 
preparations of the Huang' Chin Tang (PHY906-303503 (#11) and PHY906-284003 (#12)) 
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treated cells with High and low concentrations (indicated with H and L, respectively). The data 
were averaged based on three' repeated experiments on three different dates. Cluster analysis 
was performed based on the selected 500 genes (see text). (B) The clustering algorithm 
separated #11-L, #11H and (#12-H and #12-L) into 3 distinct clusters. Distance between 
5 clusters or resemblance coefficient is indicated by the hierarchical clustering dendrogram. 

Figure 15. Enlarged image of (A) averaged and (B) individual gene expression levels 
measured by three independent experiments. Boxl encloses genes that were down regulated in 
#1 1-L treated cells but up regulated in others, Box2 encloses the genes that were up regulated 
10 by all the herbal treatments. Box3 enclosed the genes that showed no response by #1 1-L 
treatment but were down regulated by the others. Box4 encloses the genes highly down 
regulated by low concentration herbal treatments but show mild response at high concentration 
herbal treatments. The clone ID arid putative gene name are indicated beside each gene. 

15 Figure 16. Classification of gene expression profiles in the cells treated by herbal 

medicines. Hierarchical clustering of (A) the data sets normalized with the expression data of 
the untreated control cells and (B) data sets standardized to have zero-mean and unit-variance. 
(C) The result of a non-hierarchical clustering by the self-organizing maps algorithm. 

20 Figure 17. Candidate class predictors for the classification of herbal medicines based 

on the gene expression profiles induced by the medicines. 50 class predictors with their 
expression profiles for discrirrunating #1 1 and #12 herbal preparations are shown in this figure. 
The IMAGE clone TD and putative gene name are indicated beside each gene. 

25 Figure 18. The gene expression profiles induced by a batch of a complex herbal 

preparation of five different concentrations. A 6x4 clustering of expression profiles is shown 
in (A), and the details of the gene expression profiles for the selected clusters are shown in (B). 

Figure 19. Figure 19 illustrates how the expression profiles in Figure 18 are 
categorized into three different ^groups for subsequent hamming distance calculation. 

30 Fi gure 2Q . Figure 20 shows the analysis results of gene expression profiles induced by 

five batches of a complex herbal preparation. The numbers in the table are hamming distance. 
The smaller the distance, the more similar are the expression profiles. 
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Figure 21. Shown in . (A) is a table of integrated peak intensities of 4 chemical 
constituents in HPLC analyses , of five batches of a complex herbal preparation. Two additional 
parameters, BG+B and BG/B are introduced to the table and a 6 parameter radial plot is shown 
in (B) to illustrate that one batch is more similar to a second batch #18 than to the other batches 
by the HPLC analysis. . 

Figure 22. A display of the signature genes induced by a complex herbal preparation, 
the Huang Chin Tang, in Jurkat T cells. 

Figure 23. Figure 23 illustrates the principle of identifying signature genes induced by 
individual chemical constituents in a mix of herbal medicines. The signature genes are those 
whose expression levels correlate with the amount of chemical constituents in the herbal 
medicine and that the correlation.coefficient is larger than 0.99 or smaller than --0.99. (A) 
shows that the R value between the, gene and Glycyrrhizin was 0.998, and (B) shows that the 
gene whose expression levels. :increase with the decrease of Wogonin has an R value of -0.997. 

Figure 24. The signature genes induced by the chemical constituent Albiflorin in a 
complex herbal preparation, Huang Chin Tang, in Jurkat T cell. (A) show the genes that were 
positively correlated with Albiflorin, and (B) shows the genes that were negatively correlated 
with Albiflorin. ./ 

Figure 25. Correlation of gene expression profiles to a control group. (A) is the gene 
expression profile of a control group, and (B) is the gene expression profile of a sample group. 
(C) shows the number of genes with a differential expression ratio having greater than 2-fold 
increase with concentration of herbal treatment. 

Figure 26. Clusters of expression profiles clustered by a non-hierarchical analysis 
program, wherein the program is based on a self-organizing map (SOM) principle. The X-axis 
represents the herbal concentration from low to high and the Y-axis is the gene-expression 

'VI <-i . ' • 

ratio. 

Figure 27. Figure 27 shows the induced and repressed genes commonly found in two 
batches of Huang Chin Tang. 
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Figure 28. SOM clustering results for two batches of Huang Chin Tang. (A) shows the 
SOM clustering results for the expression profiles of two batches of Huang Chin Tang. (B) 
shows that ten genes have similarly responded to the two batches, and (C) shows how the 
weighing factor decreases as cluster I and cluster j become more different. 

5 Figure 29. Calculation of S score between pairs of herbal preparations in cluster 

analysis. (A) is a tabulation of the scores, and (B) is demonstrates how 5 batches of similar 
herbal preparations are related. 

DETAILED DESCRIPTION OF THE INVENTION 

Unless defined otherwise, all technical and scientific terms used herein have the same 

10 meaning as commonly understood by one of ordinary skill in the art to which this invention 
belongs. Although any methods and materials similar or equivalent to those described herein 
can be used in the practice or testing of the present invention, the preferred methods and 
materials are described. . 
Overview of the Invention 

15 As set forth above, the present invention is directed to tools and methods useful for 

predicting the biological response Of art herbal composition. More particularly, this invention 
provides methods of creating Herbal BioResponse Array (HBR Array) databases as well as 
methods for using such databases to improve the design of effective herbal-based therapeutics. 
The goal of the present invention is the overall design, creation, improvement and use of HBR 

20 Arrays for the preparation,' testing and adrhinistration of herbal compositions, and guide 
development of new herbal compositions and novel uses of existing herbal compositions. 

Phytomics. As used herein, depending on the context in which it is used, 
"phytomics" refers to using bioinfofmatics and statistical approaches to address the qualitative 
and quantitative aspects of the components of herbal compositions or to the actual data bases 

25 wMch are developed for addressing such aspects. 

Herbal BioResponse Array. As used herein, an HBR Array constitutes a data set of 
two or more observations' or measurements associated with an herbal composition. The HBR 
Array may include qualitative arid quantitative data on the plants in the composition (plant- 
related data), marker information obtained after exposure of a biosystem to the herbal 

30 composition including a dose dependent study, and a database of BioResponse data obtained 

after exposure of a biosystem to the herbal composition. The data in any particular HBR Array 
can be statistically analyzed ■ineither .2- or 3 L dimensional space. 
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HBR Arrays may be designated as batch HBR Arrays and standardized HBR Arrays. 
Batch HBR Arrays are arrays of data associated with specific batches of an herbal 
composition. Standardized HBR Arrays are arrays of data associated with a standardized 
herbal composition. 

5 Major Data Set. As used herein, the term "major data set" refers to the data set which 

acts as the baseline set of data by which various other sets of data are compared or otherwise 
analyzed for the same or different herbal compositions. Generally, the major data set is created 
using biotechnological techniques to ascertain some genetic or protein aspect of the herbal 
compositions. Thus, the major data set will usually, but not always, be based on a genomic or 

10 proteomic set of data. For example, nucleic acid microarray results could be the major data set 
which is used to compare to other, dependent or minor data sets. 

Minor or Dependent Data Set. As used herein, the "minor data set" or "dependent 
data set" refers to one or. more data sets which are used for comparing to the major data set. 
Generally, but riot always, the minor data set will consist of information on an herbal 

1 5 composition which are collected by more traditional methods. For example, the minor, or 
dependent, data set may consist of a collection of plant-related data obtained by more 
conventional means. Examples of plant-related data include, but are not limited to, the 
genus/species of the herb(s) in the herbal composition, the particular plant parts of the herb(s) 
in the composition and the gebgraphic location where the herb(s) were located. Another 

20 example of a minor data set might consist of a set of biological responses of a cell, tissue, 
organ or organism after treatment with, one or more different amounts of the herbal 
composition. Examples of such biological data or a whole organism may include, but are not 
limited to, cell toxicity studies^ /enzyme treatment studies, growth rates, weight gain or loss, 
changes in motor skills arid changes in mental abilities. 

25 Herb. Technically speaking an herb is a small, non-woody (i.e., fleshy stemmed), 

annual or perennial seed-bearing plant in which all the aerial parts die back at the end of each 
growing season. Herbs are valued for their medicinal, savory or aromatic qualities. As the 
word is more generally used and as the word is used herein, an "herb" refers to any plant or 
plant part which has a food supplement, medicinal, drug, therapeutic or life-enhancing use. 

30 Thus, as used herein, an herb- is not limited to the botanical definition of an herb but rather to 
any botanical, plant or plant part used for such purposes, including any plant or plant part of 
any plant species or subspecies of the Metaphyta kingdom, including herbs, shrubs, subshrubs, 
and trees. Plant parts used in herbai compositions include, but are not limited to, seeds, leaves, 
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stems, twigs, branches, buds, flowers, bulbs, conns, tubers, rhizomes, runners, roots, fruits, 
cones, berries, cambium and bark. 

Herbal Composition. As used herein, an "herbal composition" refers to any 
composition which includes herbs, herbal plants or herbal plant parts. Thus, as used herein, an 
5 herbal composition is any herbal preparation, including herbal food supplements, herbal 

medicines, herbal drugs and medical foods. Examples of herbal compositions include, but are 
not limited to, the following components: a whole plant or a plant part of a single plant species; 
whole plants or plant parts of multiple plant species; multiple components derived from a 
single plant species; multiple components derived from multiple plant species; or any 

10 combination of these various, components. For a thorough review of various herbal 

compositions, see, for example, Kee. Chang Huang, The Pharmacology of Chinese Herbs . CRC 
Press (1993), herein mcoiporatefrin its. entirety. Representative examples of various herbal 
compositions are provided in the following paragraphs. 

Herbal compositions wMch include the bark of the willow tree have been used to treat 

1 5 fever since the mid-eighteenth century in England. The active ingredient in willow bark is a 
bitter glycoside called sahcin, which on hydrolysis yields glucose and salicylic alcohol. 
Aspirin (acetylsalicylic acid) and aspirin-like drugs (e.g., ibuprofen), all of which are often 
called nonsteroidal antiinflammatory drugs (NSATDs), are frequently used to treat pain, fever, 
and inflammation. Meadowsweet is another herb that contains salicylates. Treatment of 

20 arthritic and arthritic-like symptoms with willow bark or meadowsweet requires the 

consumption of large quantities of herbal teas made from these plants. The entire Populus 
species (i.e., poplar trees and shrubs) also contains salicylate precursors and poplar-buds have 
been used in antiinflammatory, antipyretic and analgesic medications. 

U.S. Patents have been issued for herbal compositions used for the treatment of various 

25 diseases and other health-rreiated problems afflicting humans and animals. For example, U.S. 

Patent No. 5,417,979 discloses a composition comprising a mixture of herbs, mcluding species 
of Stephania and Glycyrrhiza, as well as their extracts, which is used as an appetite stimulant 
and for the treatment of pain. Herbal compositions which include Glycyrrhiza uralensis have 
been found useful for treating eczema; psoriasis, pruritis and inflammatory reactions of the 

30 skin (U.S. Patent No. 5,466,452). U.S. Patent No. 5,595,743 discloses various herbal 

compositions which include licorice extract (Glycyrrhiza) and siegesbeckia, sophora, stemona 
and tetrandra herbs used for the treatment of various mammalian diseases, including 
inflammation and rheumatoid. arthritis. Ocular inflammation can be treated with a 
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pharmaceutical composition containing the plant alkaloid tetrandrine (U.S. Patent No. 
5,627,195). 

U.S. Patent No. 5,683,697 discloses a pharmaceutical composition having anti- 
inflammatory, anti-fever, expectorant or anti-tussive action, wherein the composition includes 
5 plant parts from the. species Meliq, Angepica, Dendrobium, Impatiens, Citrus, Loranthus, 

Celosia, Cynanchum and Glehnia. An herbal composition which includes extracts of the roots, 
rhizomes, and/or vegetation ofAlphinia, Smilax, Tinospora, Tribulus, Withania and Zingiber 
has been found to reduce or alleviate the symptoms associated with rheumatoid arthritis, 
osteoarthritis, reactive arthritis and for reducing the production of proinflammatory cytokines 

10 (U.S. Patent No. 5,683,698). 

Herbal compositions are available in many forms, including capsules, tablets, or coated 
tablets; pellets; extracts or tinctures; powders; fresh or dried plants or plant parts; prepared 
teas; juices; creams and ointments; essential oils; or, as combinations of any of these forms. 
Herbal medicines are adrmnistered by any one of various methods, including orally, rectally, 

15 parenterally, enterally, transdermally, intravenously, via feeding tubes, and topically. 

Herbal compositions encompassed by the present invention include herbal 
compositions which also contain non-herbal components. Examples of such non-herbal 
components include, but are-not limited to, whole insects and insect parts, worms, animal or 
insect feces, natural or petroleum oils, carbonate of ammonia, salt of tartar, liquor, water, 

20 glycerin, steroids, pharmaceuticals, vitamins, nutrient extracts, whey, salts, and gelatin. 

For oral administration, the herbal compositions disclosed may take the form of, for 
example, tablets or capsules prepared by conventional means in admixture with generally 
acceptable excipients sudhas binding agents (e.g., pregelatinised maize starch, 
polyvinylpyrroUdone or hydrbxypropyl methylcellulose); fillers {e.g., lactose, microcrystalline 

25 cellulose or calcium phosphate); lubricants {e.g., magnesium stearate, talc or silica); 

disintegrants {e.g., potato starch or sodium starch glycolate); or wetting agents {e.g., sodium 
lauryl sulphate); glidants; artificial' and natural flavors and sweeteners; artificial or natural 
colors and dyes; and sblubilizers'. The herbal compositions may be additionally formulated to 
release the active agents in a time-release manner as is known in the art and as discussed in 

30 U.S. Patent Nos. 4,690,825 and'5,055,300. The tablets may be coated by methods well known 
in the art. 

Liquid preparations for oral administration may take the form of, for example, 
solutions, syrups, suspensions', or slurries (such as the liquid nutritional supplements described 
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in Mulchandani et al, 1992 U.S. Patent No. 5,108,767), or they may be presented as a dry 
product for reconstitution with water or other suitable vehicles before use. Liquid preparations 
of folic acid, and other vitamins and minerals may come in the form of a liquid nutritional 
supplement specifically designed for ESRD patients. Such liquid preparations may be 
5 prepared by conventional means with pharmaceuticaUy acceptable additives such as 
suspending agents {e.g., sorbitol syrup, methyl cellulose or hydrogenated edible fats); 
emulsifying agents (e.g., .lecithin or acacia); non-aqueous vehicles (e.g., almond oil, oily esters 
or ethyl alcohol); preservatives (e.g., methyl or propyl p-hydroxybenzoates or sorbic acid); and 
artificial or natural colors and/or sweeteners. 

1 0 For topical adrninistration, herbal components may be combined in admixture with at 

least one other ingredient constituting an acceptable carrier, diluent or excipient in order to 
provide a composition, such as a. cream, gel, solid, paste, salve, powder, lotion, liquid, aerosol 
treatment, or the like, which is most suitable for topical application. Sterile distilled water 
alone and simple cream, omtriient and gel bases may be employed as carriers of the herbal 

15 components. Preservatives and buffers may also be added. The formulation may be applied to 
a sterile dressing, biodegradable, absorbable patches or dressings for topical application, or to 
slow release implant systems With a high initial release decaying to slow release. 



biosystem. Standardized herbal compositions are usually standardized by chemical tests well 
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Standardized Herbal Composition. As used herein, a "standardized herbal 
composition" or a "characterized herbal composition" refers to a particular herbal composition 
which is chosen as the staridard^herbal composition for evaluating batch herbal compositions 
which have the same, similar or different components as the components of the standardized 
herbal composition. Sometimes herein also, referred to as the "master herbal composition." 
Standardized herbal compositions, aire generally herbal compositions which have been well 
characterized and which demonstrate the desired biological responses in a particular 
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known to one skilled in the art and are properly stored for long term usage and reference. The 
standardized herbal composition is used to establish a standardized HBR Array based on 
observations and measurements for the plants (i:e., plant-related data), markers and 
BioResponses so as to characterize the herbal composition. 
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Batch Herbal Composition. As used herein, a "batch herbal composition" refers to 
any test herbal composition which is used to establish a HBR Array based on observations and 
measurements for the plants and markers so as to characterize the herbal composition. 
Sometimes herein also referred to as a "test" or "batch" herbal composition. Observations and 
5 measurements of BioResponses may or may not be included. The herbal compositions used to 
establish the standardized herbal composition may also be referred to as "batch herbal 
compositions" until designated as "standardized herbal compositions." 

" - Batch. As used herein, a "batch" refers to a particular quantity of an herbal 
composition which can be identified as to some particular attribute so as to distinguish it from 
10 any other particular quantity of that same herbal composition. For example, one batch of an 
herbal composition may differ from another batch of that same herbal composition in that one 
of the batches was harvested . at a different time, or in a different geographical location than the 
other batch. Other differences that distinguish particular batches may include, but are not 
limited to, the following:' 1) the particular plant part used (e.g., the root of an herb was used in 
15 one batch while the leaves of that same herb were used in a different batch); 2) the post-harvest 
treatment of the mdividual herbs or herbal composition (e.g., one batch may be processed with 
distilled water while a different batch may be processed with Hydrogen Chloride to simulate 
the acidity of the human stomach); and, 3) the relative proportions of the individual herbs in an 
herbal composition (e.g., one batch may have equal parts by weight or volume of three 
20 different herbs while another batch has proportionally more of one herb than the other two). 

Biosystem. As used herein, a "biosystem" refers to any biological entity for which 
biological responses may be observed or measured. Thus, a biosystem includes, but is not 
limited to, any cell, tissue; organ, whole organism or in vitro assay. 

Biological Activity; As used herein, the "biological activity" of an herb refers to the 
25 specific biological effect peculiar to an herbal composition on a given biosystem. 

Plant-Related Data; As used herein, "Plant-related data" refers to the data collected 
on the herbal composition including, but not limited to, data about the plants, their growing 
conditions and the handling of me plants during and after harvesting. The plant-related data 
also includes the relative proportions of the components in an herbal compositions, wherein the 
30 components may be different plant parts, different plant species, other non-plant ingredients 
(e.g., insect parts, chemical drugs) or any combinations of these variables. 

Plant-related data which may be gathered for an herbal composition includes, but is not 
limited to, the following: !) the plant species (and, if available, the specific plant variety, 
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cultivar, clone, line, etc.) and specific plant parts used in the composition; 2) the geographic 
origin of the herbs, including the longitude/latitude and elevation; 3) the growth conditions of 
the herbs, including fertilizer types and amounts, amounts and times of rainfall and irrigation, 
average microEinsteins received per day, pesticide usage, including herbicides, insecticides, 
5 miticides and fungicides, and tillage methods; 4) methods and conditions used for processing 
the herbs, including age/maturity of the herbs, soaking times, drying times, extraction methods 
and grinding methods; and 5) storing methods and conditions for the herbal components and 
the final herbal composition. 

Additionally, the standardized herbal composition may be analyzed chemically. 
10 Chemical characterization may be accomplished by any chemical analysis method generally 

known by one skilled in the art. Examples of applicable chemical analyses include, but are not 
limited to, HPLC, TLC, chemical fmgerprmting, mass spectrophotometer analyses and gas 
chromatography. 

Cell Banking System. As used herein, a "cell banking system" includes a Master Cell 
15 Bank (MCB) and a Working Cell Bank (WCB) of cells. The use of a cell banking system 
minimiz es cell variability for herbal medicine testing, and is used for all types of cells in 
nucleic acid microarray Studies. . 

Bioinformatics. As used herein, "bioinformatics" refers to the use and organization of 
information of biological interest. Bioinformatics covers, among other things, the following: 
20 (1) data acquisition and analysis; (2) database development; (3) integration and links; and (4) 
further analysis of the.resulrihg database. Nearly all bioinformatics resources were developed 
as public domain freeware until the early 1990s, and much is still available free over the 
Internet. Some companies have developed proprietary databases or analytical software. 

Genomic or Genomics. As used herein, the term "genomics" refers to the study of 
25 genes and their function. Genomics emphasizes the integration of basic and applied research 
in comparative gene mapping,;molecular cloning, large-scale restriction mapping, and DNA 
sequencing and computational analysis. Genetic information is extracted using fundamental 
techniques, such as DNA sequencing, protein sequencing and PCR. 

Gene function is determined (1) by analyzing the effects of DNA mutations in genes on 
30 normal development and health of the celL tissue, organ or organism; (2) by analyzing a 

variety of signals encoded m the DNA sequence; and (3) by studying the proteins produced by 
a gene or system of related genes. 
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Proteomic or Proteomics. As used herein, the term "proteomics", also called 
"proteome research" or "phehome", refers to the quantitative protein expression pattern of a 
genome under defined conditions. As used generally, proteomics refers to methods of high 
throughput, automated analysis using protein biochemistry. 
5 Conducting proteome research in addition to genome research is necessary for a 

number of reasons. First, the level of gene expression does not necessarily represent the 
amount of active protein in a cell. Also, the gene sequence does not describe post-tranlsational 
modifications which are essential for the function and activity of a protein. In addition, the 
genome itself does not describe the dynamic cell processes which alter the protein level either 
10 up or down. 

Proteome programs seek to characterize all the proteins in a cell, identifying at least 
part of their amino acid sequence of an isolated protein. In general, the proteins are first 
separated using 2D gels or HPLC and then the peptides or proteins are sequenced using high 
throughput mass spectrometry. Using a computer, the output of the mass spectrometry can be 

15 analyzed so as to link a geriejand the particular protein for which it codes. This overall process 
is sometimes referred to as "junctional genomics". A number of commercial ventures now 
offer proteomic services (e.^i, Pharmaceutical Proteomics™, The ProteinChip™ System from 
Ciphergen Biosystem; PerSeptive Biosystems). 

For general irdbrmation on proteome research, see, for example, J.S. Fruton, 1999, 

20 Proteins. Enzvmes. Genes: The Interplay of Chemistry and Biology . Yale Univ. Pr.; Wilkins et 
al., 1997, Proteome Research: New Frontiers in F unctional Gen omics (Principles and Practice) , 
Springer Verlag; A J. Link, 1999, 2-D Proteome Analysis P rotocals (Methods in Molecular 
Biology. 112. Humaiia'Pr-ffeimp et al:, 1999, Proteome an d Protein Analysis. Springer 
Verlag. 

25 Signal TransductionlyAs used herein, "signal transduction", also known as cellular 

signal transduction, refers tb the pathways through which cells receive external signals and 
transmit, amplify and direct mein internally. Signaling pathways require mtercommunicating 
chains of proteins that transniit the signal in a stepwise fashion. Protein kinases often 
participate in this cascade of -reactions, since many signal transductions involve receiving an 

30 extracellular chemical signal, which triggers the phosphorylation of cytoplasmic proteins to 
amplify the signal. . . 

Post-translational Modification. As used herein, "post-translational modification" is 
a blanket term used to coyer the alterations that-happen to a protein after it has been 
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synthesized as a primary polypeptide. Such post-translational modifications include, but are 
not limited to, glycosylation,. removal of the N-terminal methionine (or N-formyl methionine), 
signal peptide removal, acetylation, formylation, amino acid modifications, internal cleavage 
of peptide chains to release smaller proteins or peptides, phosphorylation, and modification of 
5 methionine. 

Array or Microarray.. As used herein, ah "array" or ,, microarray ,, refers to a grid 
system which has each position or probe cell occupied by a defined nucleic acid fragment. The 
arrays themselves are sometimes referred to as "chips", "biochips", "DNA chips" or "gene 
chips". High-density nucleic acid microarrays often have thousands of probe cells in a variety 

10 of grid styles. - 

Once the array is fabricated, DNA or protein molecules derived from a biosystem are 
added and some form . of chemistry occurs between the DNA or protein molecules and the array 
to give some recognition pattern that is particular to that array and biosystem. 
Autoradiography of radiolabeled batches is a traditional detection strategy, but other options 

15 are available, including fluorescence, colorimetry, and electronic signal transduction. 

Markers. As used herein, the term "markers" refers to any biological-based 
measurement or observation for a particular herbal composition that is characteristic of a 
particular biosystem which is being exposed to a particular batch of an herbal composition. 
The term "marker" encompasses both qualitative and qualitative measurements and 

20 observations of a biosystein. : The marker database constitutes a data set which characterizes 
gene expression patterns in response to herbal therapies, wherein the patterns show which 
genes are turned on, off, up or down in response to specific herbal compositions. Thus, 
"markers" refers to any biologically-based measurement or observation whose up- and down- 
or temporal regulations, or quahtative or quantitative changes of expression levels in a 

25 biosystem are used to characterize differential biological responses of a biosystem to an herbal 
composition. 

The particular batch of an herbal composition to which the biosystem is exposed may 
be an unknown herbal composition, a known herbal composition, or a standardized herbal 
composition. Examples of markers useful in accomplishing the present invention include, but 
30 are not limited to, molecular markers, cytogenetic markers, biochemical markers or 

macromolecular markers. Macromblecular markers include, but are not limited to, enzymes, 
polypeptides, peptides, sugars', antibodies, DNA, RNA, proteins (both translational proteins 
and post-translational proteins),:niicleic acids, polysaccharides. 
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Any marker that satisfies the definition of "marker" herein is appropriate for conducting 
the present invention. The term "markers" includes related, alternative terms, such as 
"biomarker" or "genetic marker" or "gene marker." There may be one or more primary 
markers along with secondary markers, or a hierarchy of markers for achieving the purposes of 
5 increasing the mscrimmatihg power of a HBR array. Thus, selected molecular markers may be 
combined with various other molecular, cytogenetic, biochemical or macromolecular markers 
to enable an even more accurate, extended HBR Array. 

A molecular marker comprises one or more microscopic molecules from one or more 
classes of molecular compounds, such as DNA, RNA, cDNA, nucleic acid fragments, proteins, 

10 protein fragments, lipids, fatty acids, carbohydrates, and glycoproteins. 

The establishmenti. generation and use of applicable molecular markers are well known 
to one skilled in the art. Examples of particularly useful technologies for the characterization 
of molecular markers include differential display, reverse transcriptase polymerase chain 
reactions (RT-PCR), large-scale sequencing of expressed sequence tags (ESTs), serial analysis 

15 of gene expression (SAGE), Western immunoblot or 2D, 3D study of proteins, and microarray 
technology. One skilled in the art. of molecular marker technology is familiar with the methods 
and uses of such technology: (see, e.g:,' Bernard R. Glick and Jack J. Pasternak, Molecular 
Biotechnology. Principles and Applications of Recombinant DNA, Second Edition , ASM Press 
(1998); Mathew R. Walker and Ralph Rapley, Route Maps in Gene T echnology, Blackwell 

20 Science (1997); Roe et al:, DNA Isolation and Sequencing . John Wiley & Sons (1996)James 
D. Watson et al.; Recombinant DNA. Second Edition. Scientific American Books (1992)). 

DNA, RNA and protein isolation and sequencing methods are well known to those 
skilled in the aft. Examples' of such well known techniques can be found i n Molecular 
Cloning: A Laboratory Manual 2nd Edition . Sambrook et al., Cold Spring Harbor, N.Y. 

25 (1989); Hanspeter Saluz and J. P, Jost, A Laboratory Guide to Genomic Sequencing: The 

Direct Sequencing of Native Uncloned DNA (Biomethods Vol 1) . Birkhauser (1988); and B. 
Roe et al., DNA Isolation and Sequencing . Wiley (1996). Examples of conventional molecular 
biology techniques include, but are hot limited to, in vitro ligation, restriction endonuclease 
digestion, PGR, cellular triuisformation, hybridization, electrophoresis, DNA sequencing, cell 

30 culture, and the like. Specific kits and tools available commercially for use in the present 

invention include, but are not limited to, those useful for RNA isolation, PCR cDNA library 
construction, retroviral expression libraries, vectors, gene expression analyses, protein 
antibody purification, cytotoxicity assays, protein expression and purification, and high- 
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throughput plasmid purification (see, e.g., CLONTECHniques product catalog, XIII(3), 1-32 
(1998) or www.clontech.com; Atlas™ cDNA Expression Assays product catalog (1998); 
SIGMA® product catalog ( 1997)). 

For discussions, methodologies and applications of oligonucleotide arrays, microarrays, 
5 DNA chips or biochips, see, for example, U.S. Patent Numbers 5,445,934, 5,605,662, 

5,631,134, 5,736,257, 5,74/1,644, 5,744,305, 5,795,714; Schenaet al., Parallel human genome 
analysis: Microarray-based expression monitoring of 1000 genes, Proc. Natl. Acad. Sci. USA 
93, 10614-10619 (1996); DeRisi et aL, Exploring the Metabolic and Genetic Control of Gene 
Expression on a Genomic Scale. Science 278, 680-686 (1997); Wodicka, et al., Genome-wide 
10 Expression Monitoring in Saccharomyces cerevisiae, Nature Biotechnology 15, 1359-1367 

(1997) ; Pardee, Complete Genome Expression Monitoring: The Human Race, Nature 
Biotechnology 15, 1343-1344 (1997); Schafer et al., DNA Variation and the Future of Human 
Genetics, Nature Biotechnology 16, 3.3-39 (1998); DeRisi et al., Use of a cDNA Microarray to 
Analyze Gene Expression Patterns in Human Cancer, Nature Genetics 14, 457-460 (1996); 

15 Heller et al., Discovery and Analysis of Inflammatory Disease-Related Genes Using cDNA 
Microarrays, Proc. Natl. Acad., Sci. USA 94, 2150-2155 (1997); Marshall et al., DNA Chips: 
An Array of Possibilities. Nature Biotechnology 16, 27-31 (1998); Schena et al., Microarrays: 
Biotechnology's Discovery Platform for Functional Genomics. Tibtech 16, 301-306 (1998); 
Ramsay, DNA Chips: State-bf-the-aVt, Nature Biotechnology 16, 40-44 (1998); Chee et al., 

20 Accessing Genetic Information with High-Density DNA Arrays, Science 274, 610-614 (1996); 
and Chen et al., Profiling Expression Patterns and Isolating Differentially Expressed Genes by 
cDNA Microarray System with Colorimetry Detection, Genomics 50, 1-12 (1998); P. Andrew 
Outinen et al., Characterization of the stress-inducing effects of homocysteine, Biochem. J. 
332, 213-221 (1998); and Gelbert et al., Will genetics really revolutionize the drug discovery 

25 process. Curr Opin Biotechnol 8C6Y 669-674 (1997). 

Other, more specific, references applicable to the instant invention include, but are not 
limited to, those addressing the expression technologies, such as ESTs (see, e.g., Michael R. 
Fannon, Gene expression in normal and disease states - identification of therapeutic targets, 
TIBTECH 14, 294-298 (1996)); the generation of protein profiles (see, e.g., Robinson et al., A 

30 Tyrosine Kinase Profile of Prostate Carcinoma, Proc. Natl. Acad. Sci. USA 93, 5958-5962 
(1996)); chemical and spectroscopic methods for identifying components of herbal 
compositions (Kojima et al, Saponins from Gliricidia sepium, Phvtochemistry 48(5), 885-888 

(1998) ); the determination of Afunctional antigens (see, e.g., Aris Persidis, Functional 
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antigenics. Nature Biotechnology 16, 305-307 (1998)); HPLCs (see, e.g., Milton T. W. Hearn 
(Editor), HPLC of Proteins. Peptics, and Polynucleotides: Contemporary Topics and 
A pplications f Analytical Techniques in Clinical Chemistry and Laboratory Manual! VCH 
Pub. (1991); electrophoresis (see, e.g., Westenneier et al., Electrophoresis in Practice: A Guide 
5 to Methods and App lications of DNA and Protein Separations. John Wiley & Sons (1997)); 
and cross-reactivity marker assays (see, e.g., Irving Millman et al., Woodchuck Hepatitis 
Virus: Experimental Infection and Natural Occurrence, Hepatology 4(5):8 17-823 (1984)). The 
use of structural genomics for solving the structures of all the proteins encoded for in 
completed genomes, wherein the. methodology includes high-throughput direct structure 

10 determinations and computational methods, is discussed by Terry Gaasterland, Structural 
genomics: Bioinformatics in the driver's seat, Nature Biotechnology 16, 625-627. For 
bioinformatics methodologies, see, for example, Andreas Baxevanis (Editor), Bioinformatics: 
A Practical Guide to the Analysis of Genes and Proteins . John Wiley & Sons (1998) and Luke 
Alphey, DNA Sequencing: From Experimental Metho ds to Bioinformatics (Introduction to 

15 Biotechniques Series) . Springer Verlag (1997). 

Cytogenetic parameters include, but are not limited to, karyotype analyses {e.g., 
relative chromosome lengths* centromere positions, presence or absence of secondary 
constrictions), ideograms; (i. e.; a magraminatic representation of the karyotype of an organism), 
the behavior of chromosomes during mitosis and meiosis, chromosome staining and banding 

20 patterns, DNA-protein interactions (also known as nuclease protection assays), neutron 
scattering studies, rolling circles (A.M. Diegelman and E.T. Kool, Nucleic Acids Res 
26(13):3235-324l (1998); Backert et aL, Mol. Cell. Biol. 16(ll):6285-6294 (1996); Skaliter et 
al., J. Viol. 70(2):1 132-1 136 (1996); A. Fire and S.Q. Xu, Proc. Natl. Acad. Sci. USA 
92(10):4641-4645 (1995)), and autoradiography of whole nuclei following incubation with 

25 radiolabelled ribonucleotides; 

Biochemical parameters include, but are not limited to, specific pathway analyses, such 
as signal transduction, protein- synthesis and transport, RNA transcription, cholesterol synthesis 
, and degradation, glucogenesis and glycolysis. 

Fingerprinting. As used herein, the term "fingerprinting" as used herein refers to the 

30 means of making a characteristic, profile of a substance, particularly an herb, in order to 
identify it. The term " fingeipririt" as used herein refers to the display of the result of the 
particular means employed for the fingerprinting. 
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Examples of the various types of fingerprinting means include, but are not limited to, 
DNA fingerprinting, protein fingerprinting, chemical fingerprinting and footprinting. 

DNA fingerprinting, or profiling, refers to a way of making a unique pattern from the 
DNA of particular biological source (e.g., a particular plant, plant species, genus of plant, plant 
5 part or plant tissue). The DNA fingerprint, or profile, can be used to distinguish that particular 
biological source from a different biological source. The patterns obtained by analyzing a 
batch using microarrays, ohgohucletide arrays, DNA chips or biochips are also referred to as 
"fingerprmts". 

Protein fingerprinting refers to generating a pattern of proteins in a cell, tissue, organ or 
10 organism, such as a plant, which provides a completely characteristic "fingerprint" of that cell, 
tissue, organ or organism at that time. 

Chemical fmgerprinting refers to the analysis of the low molecular weight chemicals in 
a cell and the resulting pattern used to identify a cell, tissue, organ or organism, such as a plant. 
The analysis is usually done using Gas Chromatography (GC), HPLC or mass spectrometry. 
1 5 Footprinting refers to a method of fmding how two molecules stick together. In the 

case of DNA, a protein is bound to a labeled piece of DNA, and then the DNA is broken down, 
by enzymes or by chemical attack. This process produces a "ladder" of fragments of all sizes. 
Where the DNA is protected by the bound protein it is degraded less, and so the "ladder" 
appears fainter. Footprmtingis a common technique for homing in on where the proteins that 
20 regulate gene activity actually bind to the DNA. 

The means, or methods, used to accomplish each type of fingerprinting are described in 
detail elsewhere herein. ; : 

BioResponses. As-iised '-herein, a "BioResponse" refers to any observation or 
measurement of a biolo gical response of a biosystem following exposure to an herbal 
25 composition. Sometimes herein also referred to as a "biological effect." A BioResponse is a 
qualitative or quantitative data point for the biological activity of a particular herbal 
composition! BioResponse data includes both dosage and temporal information, wherein such 
information is well known' to one skilled in the art of measuring responses of biosystems to 
various treatments. Thus, BioResponse data includes information on the specific biological 
30 response of a specific biosysfein to a specific dosage of herbal composition achriinistered in a 
particular manner for a specific period of time. 

BioResponses include, but are not limited to, physiological responses, morphological 
responses, cognitive responses, motivational responses, autonomic responses and post- 
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translational modifications, such as signal transduction measurements. Many herbal 
compositions demonstrate more than one BioResponse (see, e.g., Kee Chang Huang, The 
Pharmacology of Chinese Herbs . GRC Press (1993)). Some particular BioResponses may be 
included in more than one of the delineated groups or have aspects or components of the 
5 response that encompass more than one group. BioResponses applicable to the instant 

invention are well known to one skilled in the art. The following references are representative 
of the state of art in the field: Kee Chang Huang, The Pharmacology of Chinese Herbs, CRC 
Press (1993); Earl Mindell . Karl Mindell's Herb Bible . Simon & Schuster (1992); Goodman & 
Gilman's The Pharmacological Basis of Therapeutics- Ninth Edition. Joel G. Hardman, et. al. 

10 (eds.), McGraw Hill, Health Professions Division (1996); P. J. Bentley, Elements of 

pharmacology. A primer on drug action. Cambridge University Press (1981); P. T. Marshall 
and G. M. Hughes, Physiology of mammals and other vertebrates. Second Edition , Cambridge 
University Press (1980); Report of the Committee on Infectious Diseases , American Academy 
of Pediatrics (1991); Knut Schmidt-Nielsen, Animal Phv siologv: Adaptation and Environment. 

15 5th Edition . Cambridge University' Press (1997); Elain N. Marieb, Human Anatomy & 
Phvsiologv . Addison-Wessley Pub. : Co . (1997); William F. Ganong, Review of Medical 
Physiology (1 8th Ed\ Appletbn & Lange (1997); Arthur C. Guyton and John E. Hall, 
Textbook of Medical Physiolog y. W.B. Saunders Co. (1995). 

A "physiological response" refers to any characteristic related to the physiology, or 

20 functioning, of a biosystemu Physiological responses on a cellular, tissue or organ level 

include, but are not limited to, temperature, blood flow rate, pulse rate, oxygen concentration, 
bioelectric potential pH. value, cholesterol levels, infection state {e.g., viral, bacterial) and ion 
flux. Physiological responses ori a whole organism basis include gastrointestinal flmctioning 
{e.g., ulcers, upset stomach, mdigestion, heartburn), reproductive tract functioning {e.g., 

25 physiologically-based impotence; uterine cramping, menstrual cramps), excretory functions 
{e.g., urinary tract problems, kidney ailments, diarrhea, constipation), blood circulation {e.g., 
hypertension, heart disorders'), oxygen consumption, skeletal health {e.g., osteoporosis), 
condition of the cartilage and connective' tissues (e.g., joint pain and inflammation), 
locomotion, eyesight (e.g.* myopia, blindness), muscle tone {e.g. wasting syndrome, muscle 

30 strains), presence or absence of pain, epidermal and dermal health {e.g., skin irritation, itching, 
skin wounds), functioning of the endocrine system, cardiac functioning, nervous coordination, 
head-related health (e.g., headaches, dizziness), age (e.g., life span, longevity) and respiration 
{e.g., congestion, respiratory ailments). 
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A "morphological response" refers to any characteristic related to the morphology, or 
the form and structure, of a biosystem following exposure to an herbal composition. 
Morphological responses, regardless of the type of biosystem, include, but are not limited to, 
size, weight, height, width, color, degree of inflammation, general appearance (e.g., 
5 opaqueness, transparency, paleness), degree of wetness or dryness, presence or absence of 
cancerous growths, and the presence or lack of parasites or pests (e.g., mice, lice, fleas). 
Morphological responses on a- whole organism basis include, but are not limited to, the amount 
and" location of hair growth (e g!; hirsutism, baldness), presence or absence of wrinkles, type 
and degree of nail and skin growth, degree of blot clotting, presence or absence of sores or 

10 wounds, and presence or absence of hemorrhoids. 

A "cognitive response" refers to any characteristic related to the cognitive, or mental 
state, of a biosystem foUowing exposure to an herbal composition. Cognitive responses 
include, but are not limited fo, ! perceiving, recognizing, conceiving, judging, memory, 
reasoning and imagining: 

15 A "motivational response" refers to any' characteristic related to the motivation, or 

induces action, of a biosystem following exposure to an herbal composition. Motivational 
responses include, but are not limited to, emotion (e.g., cheerfulness), desire, learned drive, 
particular physiological needs, (<?.£., appetite, sexual drive) or similar impulses that act as 
incitements to action (e.g., stamina, sex drive)! 

20 An "autonomic response" refers to any characteristic related to autonomic responses of 

a biosystem following exposure to an herbal composition. Autonomic responses are related to 
the autonomic nervous system of the biosystem. Examples of autonomic responses include, 
but art not limited to, involuntary functioning (e.g., nervousness, panic attacks), or 
physiological needs (e.g., respiration, cardiac rhythm, hormone release, immune responses, 

25 insomnia, narcolepsy). ; 

BioResponses of cells, tissues, organs and whole organisms treated with various herbal 
compositions or herbal, components are well known in the herbal arts. For example, the herbal 
compositions Sairei-tb (TJ-1 14), alismatis rhizoma (Japanese name 'Takusha') and hoelen 
(Japanese name 'Bukuryou') were each found to inhibit the synthesis and expression of 

30 endothelin-1 in rats ; (Hattori et a'L, Sairei-to may inhibit the synthesis of endothelin-1 in 

nephritic glomeruli. Ni ppon Jinzo Gakkai Shi 39(2), 121-128 (1997)). Interleukin (IL)-1 alpha 
production was significantly promoted by treatment of cultured human epidermal keratinocytes 
with the herbal medicine Shb-saiko-to (Matsumoto et al., Enhancement of interleukin- 1 alpha 
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mediated autocrine growth of cultured human keratinocytes by sho-saiko-to, Jpn J. Pharmacol 
73(4), 333-336 (1997). Adding ShO-saiko-to to a culture of peripheral blood mononuclear 
cells obtained from healthy volunteers resulted in a dose-dependent increase in the production 
of granulocyte colony-stimulating factor (G-CSF) (Y amashiki et aL, Herbal medicine "sho- 
5 saiko-to" induces in vitro granulocyte colony-stimulating factor production on peripheral blood 
mononuclear cells, J Clin Lab Immunol 37(2), 83-90 (1992)). These researchers concluded 
that the admimstration of Shorsaiko-to may be useful for the treatment of chronic liver disease, 
malignant diseases and acute infectious diseases where G-CSF is efficacious. Plasminogen 
activator inhibitor type 1 (PAI-l)-specific mRNA expression decreased and tissue-type 

1 0 plasminogen activator (t-PA)-specific mRNA increased after treatment of human umbilical 

vein endothelial cells (HUVECs) with the saponin astragaloside IV (AS-IV) purified from the 
Chinese herb Astragalus membranaceus (Zhang et al., Regulation of the fibrinolytic potential 
of cultured human umbihcal -veih endotheUal cells: astragaloside TV down regulates 
plasminogen activator inhibitor-1 and up regulates tissue-type plasminogen activator 

15 expression, J Vase Res 34(4>, 273-280 (1997)). One component out of four components 

isolated from the roots of Panax ginseng was found to be a potent inducer of IL-8 production 
by human monocytes and THP-1 cells, and this induction was accompanied by increased IL-8 
mRNA expression (Sonoda et al.; Stimulation of interleukin-8 production by acidic 
polysaccharides from the root of panax ginseng, Immunopharmacolo trv 38(3), 287-294 

20 (1998)). By flow cytometric analysis, the expression of Fc gamma 11/111 receptors and 

complement receptor 3 (CR3) oh macrophages were found to be increased by treatment with 
the Kampo4ierbal medicine : Toki-shakuyakusan (TSS) (Cyong, New BRM from kampo-herbal 
medicine, Nippon Yakurigaku ; Zasshi 110 Suppl 1, 87P-92P (1997)). Using computer image 
analysis, Chen etal. (Image analysis for intercellular adhesion molecule- 1 expression in 

25 MRI/lpr mice: effects of Chinese herb medicine, Chung Hua I Hsueh Tsa Chih 75(4), 204-206 
(1995)) found that the distribution intensity of intercellular adhesion molecule-1 (ICAM-l), 
immunoglobulins and C3 were significantly decreased in MRL/lpr mice after treatment with 
the Chinese herb stragalih. Western blot analysis showed that tetradrine, isolated from a 
natural Chinese herbal medicine, inhibited signal-induced NF-kappa B activation in rat 

30 alveolar macrophages (Chen et al., Tetrandrine inhibits signal-induced NF-kappa B activation 
in rat alveolar macrophages, Biochem Biophvs'Res Commun 231(1), 99-102 (1997)). 

Algorithm. As used herein, an "algorithm" refers to a step-by-step problem-solving 
procedure, especially an established, recursive computational procedure with a finite number 
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of steps. Appropriate algorithms for two- and three-dimensional analyses of the plant-related, 
marker and BioResponse data sets are well known to one skilled in the computational arts. 
Such algorithms are useful in constructing the Herbal BioResponse Arrays of the present 
invention. For general information on algorithms, see, for example, Jerrod H. Zar, 
5 Biostatistical Analysis, second edition, Prentice Hall (1 984); Robert A. Schowengerdt, 

Techniques for image processing and classification in remote sensing, Academic Press (1983); 
Steven Gold et al., New Algorithms for 2D and 3D Point Matching: Pose Estimation and 
Correspondence. Pattern Recognition . 31(8):1019-1031 (1998); Berc Rustem. Algorithms for 
Nonlinear Programming and Multiple-Objective Decisions , Wiley-Interscience Series in 

10 Systems and Optimization, John Wiley & Sons (1998); Jeffrey H. Kingston, Algorithms and 
Data Structures: Design. Correctness. Analysis. Internati onal Comp uter Science Series . 
Addison-Wesley Pub. Co. (1 997); Steven S. Skiena, The Algorithm Design Manual . Springer 
Verlag (1997); and Marcel F.Neuts, Algorithm Probability: A Collection of Problems 
(Stochastic Modeling) . Chapman & Hall (1995). For information more specific to the 

15 application of algorithms to genetic-based data, see, for example, Dan Gusfield, Algorithms on 
Strings. Trees, and Sequences: Computer Science and Computational Biology . Cambridge 
University Press (1997); Melanie Mitchell, An Introduction to Genetic Algorithms (Complex 
Adaptive Systems) , MIT Press (1996); David E. Goldberg, Genetic Algorithms in Search. 
Optimization and Machine Learning , Addison- Wessley Pub. Co. (1989); Zbigniew 

20 Michalewicz, Genetic Algorithms + Data Structures = Evolution Programs . Springer Verlag 
(1996); Andre G. Uitterhnderi : arid Jan Vijg . Two-Dimensional DNA Typing: A Parallel 
Approach to Genome Analysis , Ellis Horwood Series in Molecular Biology, Ellis Horwood 
Ltd. (1994); and Pierre Baldi and Soiren Brunak, Bioinformatics: The Machine Learning 
Approach (Adaptive Computation and Mach ine Learning ), MIT Press (1998). 

25 Combinatorial Chemistry. As used herein, "combinatorial chemistry" refers to the 

numerous technologies used to create hundreds or thousands of chemical compounds, wherein 
each of the chemical compounds differ for one or more features, such as their shape, charge, 
and/or hydrophobic characteristics; Combinatorial chemistry can be utilized to generate 
compounds which are chemical Variations of herbs or herbal components. Such compounds 

30 can be evaluated using the methods <)f the present invention. 

Basic combinatorial chemistry concepts are well known to one of ordinary skill in the 
chemical arts and can also be found in Nicholas K. Terrett, Combinatorial Chemistry (Oxford 
Chemistry. Masters^ . Oxford Univ. Press (1998); Anthony W. Czarnik and Sheila Hobbs 
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Dewitt (Editors), A Practical. Guide to Combinatorial Chemistry . Amer. Chemical Society 
(1997); Stephen R Wilson (Editor) and Anthony W. Czamik (Contributor), Combinatorial 
Chemistry: Synthesis and Application . John Wiley & Sons (1997); Eric M. Gordon and James 
F. Kerwin (Editors), Combinatorial Chemistry and Molecular Diversity in Drag Discovery . 
5 Wiley-Liss (1 998); Shmuel Cabilly (Editor), Combinatorial Peptide Library Protocols 

(Methods in Molecular Biology) . Human Press (1997); John P. Devlin, High Throughput 
Screening . Marcel Dekker (1998); Larry Gold and Joseph Alper, Keeping pace with genomics 
through combinatorial chemistry; Nature Biotechnology 15, 297 (1997); Aris Persidis, 
Combinatorial chemistry. Nature Biotechnolog y 16, 691-693 (1998). 

10 EXAMPLES 

Example 1. Establishing a Standardized HBR Array for Selected Herbal 
Compositions. : - ;; t ;V.' 

The basic scheme for estabhshing a Standardized HBR Array is provided in Figure 1 . 
Definitions of each component of the schematic are provided above. 

15 Following selection of an herbal composition of interest, data is collected for various 

traits associated with the herbal composition, including, but not limited to plant-related 
characteristics and marker and BioResponse information. 

Plant-related data includes, but is not limited to, the plant species, specific plant parts, 
geographic origin of the plants in the herbal composition, the growth conditions of the plants, 

20 the processing methods used to 'prepare the herbal components, storage methods and 

conditions, and various chemical analyses of the herbal composition. Marker information 
includes qualitative and quantitative data for markers collected after exposure of a biosystem to 
the herbal compost. Apphcable markers include, but are not limited to, molecular markers, 
cytogenetic markers, biochemical'markers and macromolecular markers. BioResponse 

25 information includes quahtativle' and quantitative data for biological responses collected after 
exposure of a biosystem .to : the herbal composition. 

Each type of data (e.gt, chemical, marker, BioResponse) can be obtained using one or 
more assays "on the same, similar, substantially similar, or different batches of the herbal 
composition of interest. Such different assays can be conducted at the same or different times. 

30 In addition, data can be collected for the same or different markers at the same or different 
times. Similarly, BioResponse* data can be collected for the same or different biological 
responses at the same or different times. Thus, collection of the data for the HBR Array is 
either collected at one time' or collected on an on-going basis. Where a biosystem is exposed 
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to an herbal composition so as to collect data, information is recorded on the aolministered 
dosages of the herbal composition as well as treatment times. BioResponse data may also 
consist of post-translational modifications, such as measurements of signal transduction. 

After collection of two or more types of data (e.g., data for two or more markers and a 
5 BioResponse; data for plant-related traits and data for a BioResponse), the data is analyzed 
using algorithms so as to create 2- and/or 3-dimensional Herbal BioResponse Arrays. 

Various statistical parameters may be calculated for the HBR Array and may become 
part pf the HBR Array data set. These statistical parameters may include, but are not limited 
to, means, standard deviations, correlation or match (or mismatch) matrices, ratios, regression 
10 coefficients, and transformed values (e.g., arcsin percentage transformations of the raw data). 
Thus, the HBR Array may consist of the raw data as well as certain calculations, distributions, 
graphical presentations and other data manipulations associated with the raw data. Particular 
examples of such information include, but are not limited to, digital images, scatter graphs, 
cluster analyses and large scale gene expression profiles for marker data. 
15 The total accumulated data and resultant analyses constitute a standardized HBR Array 

for the particular herbal composition used to establish the HBR Array data set. Due to the 
iterative nature of the process 'used to establish and maintain an HBR Array for an herbal 
composition, such arrays can be viewed as either static at any one point in time or dynamic 
over time. : 
20 The resulting analyses can identify subsets of the standardized HBR Arrays which are 

correlated (positively or negatively) or associated (i.e., showing a general trend) with one or 
more specific biological activities of any particular herbal composition. 

Example 2. Establishing a Batch HBR Array for Batch Herbal Compositions. 

The basic scheme for establishing a HBR Array for a batch of an herbal composition is 
25 provided in Figure 2. Definitions of each component of the schematic are provided above. 
The procedure for estabhshing such an array is the same as that set forth immediately above 
for the standardized HBR Artay; 

Generally,- the amount, of data collected for a batch HBR Array will be less than that 
collected to establish a standardized HBR Array. However, data collected for a batch herbal 
30 composition may be added to an .established HBR Array or used to establish a new 
standardized HBR Array. . : - 

Generally, the only data collected for a batch herbal composition is that data which has 
been found to be highly -correlated or associated with the desired biological activities of the 
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herbal composition beingitested. For example, if it has been determined that a particular subset 
of plant-related and marker data is highly correlated to a desired biological activity of a 
particular herbal composition (based on the standardized HBR Array data and analyses 
discussed above), it is only necessary to test the batch herbal composition for that subset of 
5 traits in order to determine v/hether or not the batch has the desired biological activity. By 
comparing the data obtainedfor that subset of traits obtained from the batch (i.e., the batch 
HBR Array) with me standardized HBR Array . for that particular herbal composition, one 
skilled in the art can determine whether or not that particular batch has the desired biological 
activity. ■ " 

10 Example 3. Establishing and Using a Major Data Set. 

The basic scheme for.estabhshing and using a major data set for an herbal composition 
is provided in Figure 3 . DefMtiohs of each component of the schematic are provided above. 

The first step is the estabhshment of a major data set for a selected herbal composition 
or batch herbal compositions This is accomplished by exposing a biosystem to the herbal 

1 5 composition and collecting the resultant marker information which will constitute the major 
data set. In most, but not instances, the major data set will consist of genomics and/or 
proteomics data in the form of an array, such as an array obtained with a DNA biochip. 

Next, the major data set is analyzed to see if differential expression/results have been 
obtained for the tested herbal composition. Differential expression/results are necessary in 

20 order to generate meariingful algorithms in the next step. Examples of such differential 
expression/results include, but are hot limited to, indications that certain genes are up- or 
down-regulated in response to exposure to the herbal composition or that the levels of certain 
proteins have been increased or decreased in response to the exposure. 

If no meahingmi or useful differential expression/results are obtained, then it is 

25 necessary to repeat the exposure and marker collection step. If it is believed that experimental 
error lead to the lack of a adequate result the first time then the exposure/data collection step 
can be repeated with all of the variables the same as the first time (e.g., same biosystem, same 
marker set, same experimental protocol, etc.). However, it may be necessary to vary the 
biosystem sampling (e.g:, type of cell utilized, stage of cell growth), use a different marker set 

30 and/or change the experimental protocol in order to get differential expression/result.. 
Example 4. Using 1 HBR Array Information. 

The HBR Array information discussed herein can be used for many different purposes 
including, but not limited to, 'the following: 1) evaluating the components of an herbal 
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composition; 2) predicting the BioResponse of an herbal composition; 3) determining which 
marker information is most highly correlated with a particular BioResponse of an herbal 
composition; 3) determining what data set of information (i.e., plant-related data, marker data, 
and BioResponse data) is most correlated with a particular BioResponse of an herbal compost; 
5 4) determining which type of biosystem is best for evaluating the biological activity of an 

herbal composition; 5) adjusting or changing the components of a herbal composition so that 
the HBR Array of that herbal composition corresponds to a standardized HBR Array for the 
same or substantially the same herbal composition; 6) adjusting or changing the components of 
an herbal composition so that the herbal composition will have the desired biological activity; 

10 7) measuring the relatedness of different herbal compositions; 8) creating and updating 
standardized HBR Arrays; 9) identifying specific components (e.g., plant parts, proteins, 
molecules) which retain the desired biological activity of an herbal composition; 10) 
deterrnining which components of an herbal composition can be eliminated while mamtaining 
or improving the desired biological activity Of the herbal composition; 1 1) identifying one or 

15 more previously unknown biological activities for an herbal composition; 12) aiding in the 
design of therapeutics whichinclude herbal and non-herbal components, such as chemically- 
synthesized drugs or pharmaceuticals and 13) utilizing the HBR Array information to 
complement combinatorial chemistry methods of designing therapeutics. Each of these 
embodiments of the present invention can be accomplished by one skilled in the applicable art 

20 using the methods and tools provided herein. 
Example 5. Quality Control. 

The HBR Array technology of the present invention is used to correlate or to determine 
a substantial equivalence' of a . specific batch of an herbal composition (single herb or multiple 
herbs of a formula) to a standardized, or master, batch of a same or substantial similar herbal 

25 composition. The HBR Arrays utilized in this process include the acceptable range of 

quantitative variation for each' of the biological effects (i.e., BioResponse), and possibly a 
global score composed of weighted values assigned to each of the biological effects, which 
may consist of markers from multiple biochemical pathways of a biosystem, 

"Data mining" refers to aprocess used to determine or select which subset of biological 

30 effects is the minimum number of biological effects required in any specific HBR Array. The 
information for data mining results from exposing a biosystem (e.g., a cell line) in a dose 
dependent manner to a standardized herbal composition to establish a standardized HBR 
Array. This standardized HBR Array can then be compared to various HBR Arrays 
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established for test herbal compositions. These test herbal compositions include, but are not 
limited to, different batches prepared at different dates; different batches prepared from raw 
herbs collected at different times; and different batches prepared from raw herbs collected at 
different locations. . / 

Example 6. Improving an Herbal Composition or Identifying New Uses for an 
Herbal Composition. 

HBR Arrays are generated by exposing biosystems to either extracts from individual 
herbs of a formula, or to extracts from the whole formula, and examining the biological effects 
of the extracts. The observed biological effects can be from multiple biochemical pathways of 
a biosystem and/or from multiple tissues of an animal, wherein various markers are evaluated 
for their corresponding qualitative and/or quantitative changes. The resulting HBR Arrays can 
be compared to novel HBR Arrays or to similar HBR Arrays from different herbal 
compositions or herbal compositions prepared by different processes. This procedure is useful 
for selecting a given set of biological effects and the minimum number of markers required to 
predict that a given batch herbal composition has the given set of biological effects. 

In order to construct HBR Arrays, one skilled in the art utilizes various data mining 
tools including, but are not limited to, statistical analyses, artificial intelligence, and database 
research on neural work. The statistical methods of choice include, but are not limited to, basic 
exploratory data analysis (EDA), graphic EDA (such as bushing) and multivariate exploratory 
techniques (e.g., cluster analysis, discrirnmating factor analyses, stepwise linear on non-linear 
regression, classification tree) (see, e.g., STATISTICA™, software packages from StatSoft, 
Tulsa, OK 74104; Tel: 918-749rll l9; Fax: 918-749-2217; www.statsoft.com). 

Data mining tools are used to explore large amounts of HBR Array data in search of 
constructing an HBR Array and consistent pattern within, between or among various HBR 
Arrays. The procedure consists of exploration, construction of an HBR array, and validation. 
This procedure is typically repeated iteratively until a robust HBR Array, or standardized HBR 
Array, is identified. 

Example 7. Establishing a Standardized HBR Array for Ginseng Recipes. 

For the purposes of this example, standard ginseng is chosen to be Panax Ginseng C.A. 
Meyer Gl 15 grown either in Manchuria or in Korea. The climate for growth is between -10 to 
+ 1 0°C with an annual rainfall of 50- 1 00 cm (see Huang in The Pharmacology of Chinese 
Herbs . (1993) pp 21-45, CRC Press, Boca Raton, FL, fully incorporated by reference). Ginseng 
batches will first be characterized by geographic origin, species, plant part (e.g., rhizome, root, 
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leaf skin, seed, bud and flower); growth conditions, processing methods and storage conditions 
both before and after processing. Verification of chemical content for these batches will be 
performed by qualitative HPLC analysis for determination of ginsenoside saponins (e.g., Ro, 
Ral, Ra2, Rbl, Rb2, Rb3, Rc, Rgl, Rg2, Rd, Re, Rf, Rhl, Rh2, NG-R2 and Z-Rl), including 
5 TLC quahtative analysis for lipophilic constituents (see, Elkin et al, Chumg Kuo Yao Li 

Hsueh Pao (1993) 14: 97-100 and Yoshikawa et al, Yakngakn Zasshi (1993) 1 13: 460-467). 
The saponin content of different herbs should be between 2.1 and 20.6% (by weight) 
depending on the species (see Table 1). These data will then be stored, preferably in the 
memory of a computer processor, for further manipulation. 
10 Table 1. Saponin Content of Different Ginseng Herbs.* 



Species 


Total saponins (% by weight) 


Panax ginseng C.A. Meyer . 


2.1-4.4% 


Panax quiquefolius 


4.9% 


Panax notoginseng and Panax japonica 


13.6-20.6% 


Panax japonica var. major ' ' 


9.34% 



*firom Huang in The Pharmacology of Chinese Herbs . (1 993) page 29, CRC Press, Boca Raton, FL. 



Expression biomarkers for standard ginseng (i.e., Gl 15) include the following: IL-8, 
IL-2, GM-CSF, NfkB, ICAM-1, interferon gamma, choline acetyl transferase, trk A, nerve 

15 growth factor (Kim et a/... Planta Med (1998) 64: 110-115; Sonodaet al, 

ImmunopharmacQlog y (1998) 38: 287-294; Baum et al, Eur J Appl Physiol (1997) 76: 165- 
169; Iwangawa et al, Free. Radic Biol Med (1998) 24: 1256-1268; Rhind et al, Eur J Appl 
Phvsiol (1996) 74: 348-360). Alternatively, for a broader batch size, the 400,000 
oligonucleotide group/1.6 cm 2 chip of Affymetrix can be used (U.S. Pat. No.5,556,752). The 

20 expression biomarkers for standard ginseng will be prepared by nucleic acid microarray 

technology using either photolithography, mechanical microspotting or ink jet application (see 
Schena et al, TIBTECH (1998) 16: 301-306). Selected sets of ceils will be contacted with 
standard ginseng for varying periods of time, under varying conditions to generate multiple 
microarray sets. The microarray sets will then be analyzed by hybridization-based expression 

25 monitoring of biochemical extracts via deduction of steady state mRNA levels from 

fluorescence intensity at each position on the microarrays (Schena et al, Science (1995) 270: 
467-470; Schena et al, Proc Natl Acad Sci USA (1996) 93: 10614-10619; Lockhart et al Nat 
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Biotechnol (1996) 14: 1675-1680; -DeRisi et al, Nat Genet (1996) 14: 457-460; Heller et al, 
Proc Natl Acad Sci USA (1997) 94: 2150-2155). The array data sets are then input into 
algorithms to generate statistical expression biomarker values for standard ginseng. 
Biochemical biomarkers for standard ginseng include quantitative analysis for increases in 
5 cycloheximide sensitive -[ 3 H]-leucine incorporation proportional to protein synthesis and [ 3 H]- 
mymidine incorporation reflective of mitosis, (see Yamamoto et al, Arznejmittelforschung 
(1977) 27: 1 169-1 173). For biochemical biomarkers, bone marrow cells will be contacted with 
standard ginseng for varying time periods under varying conditions in the presence of [ 3 H]- 
thynudine (for DNA synthesis) or in the presence and absence of cycloheximide and [ 3 H]- 

1 0 leucine (for protein synthesis) to perform multiple quantitative analysis of biochemical 

biomarkers (i.e., BBM sets).. The BBM sets are then input into algorithms to generate statistical 
biochemical biomarker values, for standard ginseng. Statistical data will then be stored, 
preferably in the memory of a computer processor, for further manipulation. 

Biological response, of a biosystem (i. e. , BioResponses) will be determined using cells 

1 5 and whole animals. For cells*: ginseng batches will be exposed to specific cell types, including, 
but not limited to, fibroblasts, macrophages, monocytes, PMNL, LAK cells, B16-F10 
melanoma cells, THP-1 cells and hippocampal neurons at a concentration of 0.5 mg/ml to 100 
mg/ml. For animal treatments, 0.5-100 mg/kg of ginseng herbal extract will be administered 
orally, by intraperitoneal injection or subcutaneous injection. 

20 To determine a biological response of a biosystem to standardized ginseng, human 

ovarian cancer cells will be inoculated into nude mice, which results in the formation of 
palpable tumors. After tumor formation the mice will be treated by co-administration of cis- 
diammmecichloroplatinum and standard ginseng. Mice will be examined for tumor growth 
inhibition, increase in survival time and lowered adverse side-effects on hematocrit values and 

25 body weight (Nakata et al, Jpn J Cancer Res (1998) 89:733-740). The assay will be repeated 
using various concentrations of standard ginseng to generate measures of central tendency, 
dispersion and variability for each variable. 

The data collected will then be subjected to multidimensional analysis to generate 
multivariant normal distribution sets as a means of determining a baseline correlation between 

30 biological activity and standard ginseng (see Zar, J. H., in Bjostatistical Analysis, 2nd ed. 
(1984), pp 328-360, Prentice -liall, Englewood Cliffs, NJ). A second independent 
determination of a biological response of a biosystem to standard ginseng will be the effect of 
standard ginseng on physical performance during exercise. Rats will be treated for 4 days with 
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standard ginseng at various' concentrations (between 0.5-100 mg/kg/day) and animals will be 
tested for increased plasma free fatty acid level and maintenance of glucose level during 
exercise at approximately 70% V02max (see Wang etal, Planta Med (1998) 64130-133). 
The data generated will becollected and then subjected to multidimensional analysis to 
5 generate multivariant n6rmal distribution sets as a means of .determining a baseline correlation 
between biological activity and standard ginseng (see Zar, J. H., in Biostatistic al Analysis ?.nd 
ed, (1984), pp 328-360, Prentice Hall, Englewood Cliffs, NJ, Herein, fully incorporated by 
reference). The distribution sets for each BioResponse are then put into algorithms to generate 
statistical values for standard ginseng. Statistical data will then be stored, preferably in a 

1 0 memory of a computer processor, for further manipulation. 

Each of these steps (i. e., chemical analysis, generation of biomarker information and 
determination of responses; of a biosystem) is reiterated to generate a large database of 
statistical values. These values are compiled and input into an algorithm to generate 2- and 3- 
dimensional Herbal Response Arrays (HBR Array) for standardized ginseng. Through 

15 reiteration, the resulting arrays (^.'/Standardized Arrays) display the highest correlation 
between composition (mcludihg growth conditions), biomarker information and biological 
response for standardized ginseng.' By determining two or more known associated variables 
for composition and biomarker information values via display on an HBR Array for a test 
batch, the values for biological response variables can be predicted for the test batch by 

20 comparing test values against Standardized HBR Array values for standardized ginseng. The 
resulting prediction will be used to evaluate the quality of a given ginseng batch without 
necessitating the use of an observed biological response of a biosystem (see Example 2). 

Example 8. Evaluation of a Selected Herbal Composition of Ginseng Using a 
Subset of Variables Correlated with a Specific Biological Response. 

25 To evaluate the quality of a test batch herbal composition, data is first collected 

concerning the plant-related parameters for the herbs in the selected herbal composition (e.g., 
plant species, plant parts, geographic origin, growth conditions, processing methods and 
storage conditions). The selected herbal composition is then manipulated such that chemical 
analysis can be performed to determine the chemical content of the herb (see Elkin et al., 

30 Chume Kuo Yao T.i Hsneh Pah (1993) 14: 97-100 and Yoshikawa et al, Yakugaku Zasshi 

(1993) 113: 460-467). Previously obtained ginseng data has demonstrated a strong correlation 
between oxygen consumption during aerobic exercise performance and the presence of a 
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subset of saponin components, especially Rgl and Rbl (Wang et al, Planta Med (1998) 64: 
130-133). "". 

The test batch is then exposed to test cells including, but not limited to, fibroblasts, 
macrophages, monocytes!, PMNL, LAK cells, B16-F10 melanoma cells, THP-1 cells and 
5 hippocampal neurons at a concentration of 0.5 mg/ml to 100 mg/ml to determine expression 

biomarker values. mRNA is isolated from exposed cells which is subsequently manipulated to 
serve as a substrate for hybridization-based expression monitoring of biochemical extracts 
using microarrays comprising IL-8, TL-2 and Interferon gamma cDNA (Schena et al, Science 
(1995) 270: 467-470; Schena et all, Proc Natl Acad Sci USA (1996) 93: 10614-10619; 

10 Lockhart et al, Nat Biotechnol (1996) 14: 1675-1680; DeRisi et al, Nat Genet (1996) 14: 457- 
460; Heller et al, Proc Natl Acad Sci USA (1997^ 94: 2150-2155). Previously obtained 
ginseng data has demonstrated a strong correlation between oxygen consumption during 
aerobic exercise performance . and the induction of the expression biomarkers IL-8, IL-2 and 
Interferon gamma in test cells (Venkatraman et al, Med Sci Sports Exerc (1997) 29: 333-344 

15 and Wang et al, Planta Med (1998) 64: 130-133). For biochemical biomarkers, rat bone 
marrow cells will then be exposed to the test batch and assayed for [ 3 H]-mynudine 
incorporation reflective of mitosis. Previously obtained ginseng data has demonstrated that 
Rbl and Rgl show a strong correlation with DNA synthesis in rat bone marrow cells 
(Yamamoto et al, ArT^eimittelforschung (1978) 28: 2238-2241). 

20 After reiterative analysis, data from each assay will be input into an algorithm to 

generate a test HBR array'for the selected herbal composition based on the enumerated plant- 
related data, including chemical analyses, and data concerning the subset of biomarkers. The 
quality of a test batch will be determined by comparing test HBR and standard ginseng 
Standardized HBR Array variables directed toward analysis of the above observations and 

25 subsets, wherein the demonstration bfthe induction of EL-2, DL-8 and INF gamma mRNA in 
vitro and an increase in [ 5 H].-mymidine incorporation in rat bone marrow cells (including data 
collected on growth conditions, origin, and verification of the saponins Rgl and Rbl) is 
predictive of an equivalent BibResponse effect of the test batch on oxygen consumption as that 
exhibited by standard ginseng. Based on this procedure it can be determined whether or not 

30 the test batch is of a similar or different quality than that of the standard for the given 
biological response or biological response of interest. 
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Example 9. Establishing a Standardized HBR Array for Huang Ling (HL) 
Recipes. 

For the purposes of this example, standard huang ling (HL) is chosen to be Coptis 
chinesis France, from southwest Asia, wherein growth conditions are well known to one 
5 skilled in the art (see Huang in The Pharmacology of Chinese Herbs. (1993), pp 69 and 287- 
288, CRC Press, Boca Raton, FL). Dried rhizomes of Coptis chinesis France will be verified 
for chemical content by quantitative chemical analysis for determination arsenic, berberine, 
caeruleic acid, columbamineVcopsine, cpptine, coptiside-I, coptiside-II, coptisine, coreximine, 
epiberberine, ferulic acid, greehlandicine, isocoptisine, lumicaerulic acid, magnoflorine, 

10 oxybererine, thalifendine, umbellatine, urbenine, worenine, palmatine, jafrorrhizine and 

colubamine (see also Zhu M., Chung Yao Tung Pao (1984) 9: 63-64). Content of the alkaloid 
berberine of different herbs should be between 7-9% (by weight). These data will be stored, 
preferably in the memory, of a computer processor, for further manipulation. 

Expression biomarkers for standard HL include the following: NfkB; bcl-2 analog, Al; 

15 zinc finger protein, A20; ILr2 receptor; cell cycle probes; c-Ki-ras2; growth regulators probes 
and glucocorticoid receptor- dependent apoptosis probes (see Chi et al, Life Sci (1994) 54: 
2099-2107; Yang et al. Naunvn Schmiedebergs Arch Pharmacol (1996) 354: 102-108; Miura 
et al., Biochem Pharmacol / 1997) 53; Chang K.S., J Formos Med Assoc (1991) 90: 10-14). 
Alternatively, for a broader batch size, the 400,000 oligonucleotide group/1.6 cm 2 chip of 

20 Affymetrix can be used (U.S. Pat. No.5,556,752). The expression biomarkers for standard HL 
will be prepared by microarray technology as described in Example 1, including analysis and 
statistical data generation. Bipchemicai biomarkers for standard HL include increase in 
glucocorticoid receptor and inbibitiph of alpha-fetoprotein secretion in HL exposed HepG2 
cells (see Chi et al, Life Sci (1994) 54: 2099-2107). BBM sets are generated and analyzed as 

25 described in Example 1 . Statistical data will then be stored, preferably in the memory of a 
computer processor, for further manipulation. 

Biological response of a biosystem will be determined using cells and whole animals. 
Batches of the selected herbal composition will be exposed to specific cell types, including but 
not limited to, human HepG2 hepatoma cells, human embryonal carcinoma cells and 

30 thymocytes at concentrations from 0.1- lOOmg/ml. For animal treatments 0.1mg-2g/kg of 

Coptic herbal composition (i.e., HL) will be administered orally, by intraperitoneal injection or 
subcutaneous injection. To determine a biological response of a biosystem to standardized HL, 
human embryonal carcinoma 'clone, NT2/D1 is exposed to various concentrations of standard 
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HL and cells will be examined for differentiation into cells with neuronal-like cell morphology 
(Chang K.S„ I Formos Med Assoc (1991) 90: 10-14). The assay will be repeated to generate 
measures and analysis will be performed as described for ginseng in Example 1. A second 
independent determination of a biological response of a biosystem to standard HL will be the 
5 effect of standard HL on diarrhea due to enterotoxigenic Escherichia coli (ETEC). Patients 
with active diarrhea due to ETEC will be treated with various concentrations of HL (e.g., 
2g/kg) and stool volumes will be determined (see, e.g., Rabbani G.H., Dan Med Bull (1996) 
43: 173-185). The assay will be repeated to generate measures and analysis will be performed 
as described for ginseng in Example 1. The distribution sets for each biological system are 

10 then put into algorithms to generate statistical values for standard HL. Statistical data will then 
be stored, preferably in the memory of a computer processor, for further manipulation. 

Lastly, as in Example T, the steps are reiterated to generate HBR arrays for 
standardized HL, wherein the resulting HBR arrays will then be used to predict biological 
activity and evaluate batch quality: Using this method, a Standardized HBR Array can be 

15 generated and updated pefipdically. 

Example 10. Evaluation of a Selected Herbal Composition of Huang Ling Using a 
Subset of Variables Correlated with a Specific Biological Response. 

To evaluate the quality of a selected test batch of an herbal composition of Huang 
Ling, data is first collected concerning the plant-related characteristics (e.g., plant species, 

20 plant parts, geographic origin, growth conditions, processing methods and storage conditions). 
The herbal composition is then manipulated such that chemical analysis can be performed to 
determine the chemical content of the composition (see also Zhu M. T Chung Yao Tung Pao 
(1984) 9: 63-64). ' " . ' /. 

Previously obtained HL data has demonstrated terminal differentiation of human 

25 embryonal carcinoma clones into heUronal-like cells is strongly correlated with the presence of 
berberine (see Chang K.S., JFormos Med Assoc (1991) 90: 10-14). The test batch is then 
exposed to test cells including; human embryonal carcinoma clone, NT2/D1 at a concentration 
starting at a non-toxic' concentration (determination of which is within the skill of the ordinary 
artisan). mRNA is isolated firbm exposed cells which is subsequently manipulated to serve as 

30 substrate for hybridization based expression monitoring of biochemical extracts using 

microarrays comprising TL-2 receptor and NfkB; (see Chi et ah, Life Sci (1994) 54: 2099- 
2107; Yang et al, Naunyn Schmiedebergs Arch Pharmacol (1996) 354: 102-108; Miura et aL, 
Biochem Pharmacol (1997)53; Chang K.S., J Formos Med Assoc (1991) 90: 10-14; U.S. Pat. 
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No.5,556,752), and which can be used to determine down regulation of c-Ki-ras2 gene 
expression in said cells. Previously obtained HL data has demonstrated terrninal 
differentiation of human embryonal carcinoma clones into neuronal-like cells is strongly 
correlated with induction of mitogen probes and down regulation of c-Ki-ras2 gene expression 
5 (see Chang K.S., J Formos Med Assoc (1991) 90: 10-14). 

For biochemical markers, HepG2 cells are exposed to the test composition and cells are 
assayed for increase in glucocorticoid receptor and inhibition of alpha-fetoprotein secretion 
(see- Chi et al, Life Sci (1994) 54: 2099-2107). Previously obtained HL data has demonstrated 
that inhibition of glucocorticoid induced apoptosis is strongly correlated with berberine-type 
10 alkaloids (see Miura et al . ■ Biochem Pharmacol (1997) 53: 13 15-1322). After reiterative 
analysis, data from each assay will be input into an algorithm to generate a test HBR array 
based on the enumerated observational data, chemical data and data concerning the subset of 
biomarkers. 

The quality of a test batch Will be determined by comparing test HBR and standard HL 
15 HBR Array variables directed toward analysis of the above observations and subsets, wherein 
the demonstration of the induction of IL-2 receptor and NfkB, the down regulation of c-Ki- 
ras2 gene expression, an increase in glucocorticoid receptor and inhibition of alpha-fetoprotein 
secretion for HepG2 cells (to including data collected on growth conditions, origin, and 
verification of berberine alkaloid) is predictive of an equivalent BioResponse effect of the test 
20 batch on terminal differentiation of human embryonal carcinoma clones into neuronal-like 
cells and inhibition of dexamethasone induced apoptosis as that exhibited by standard HL. 
Based on this procedure it can be detennined whether or not the test batch is of a similar or 
different quality than that of the . standard. 

Example 11. Evaluation, of Xiao Chai Hu Tang (sho-saiko-to) Using Two 
25 Bioassays. 

To evaluate the . quality of three sources of Xiao Chai Hu Tang, two bioassays were 
used: 1) cell growth inhibition : and 2) hepatitis B virus secretion from infected cells. The Xiao 
Chai Hu Tang composition is made from a mixture of 6-7 herbal plants {Radix Bupeuri, 
Rhizoma Pinelliae, Rhizoma Zingiberis, Radix Scutellariae, Fructus Ziziphi, Codonopsis 
30 Pilosula, Radix Ginseng and Radix Glycyrrhizae, see Table 2 for relative amounts, by weight). 
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Table 2. Composition of Xiao Chai Hu Tang. 





Plant Species 




Radix 
buplettri 


Rliizoma 
pinelliae . 


Rliizoma 
zingiberis 


Radix 
Scutellariae 


Fructus 


Codonopsis 
pilosula 


Radix 
ginseng 


Radix 

glycyrrhizae 


Relative Amount by Weight 


Singapore 


1 




\ 0.375- 


; 0.375 


0.375 




0.375 


0.375 




1 


0.717 


' .b.sfi 


0.492 




0.429 




0.288 


Taiwan 


1 


0.25 . 


0.375 


0.375 


0.25 




0.375 


0.375 



The three 'Recipes" 'originate in either Singapore, Korea or Taiwan. Batches were 
evaluated for toxicity and : for the ability to inhibit hepatitis B virus as detected by DNA 
5 quantitation or detection of hepatitis B surface antigen (HbsAg) (see Dong et al., Proc Natl 
Acad Sci USA (1991) 88: -8495-8499). • 

Briefly, one gram of preparation was added with 10 ml of water. The mixture was 
boiled for 30 minutes. The supernatant was collected after centrifugation and filtered through 
a 0.22 j_tm filter. Two cell types were used: a) 2.2.15 cells which secrete hepatitis B virons 

10 (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Ac ad Sci USA (1987) 84: 1005- 
1009) and b) HepG2 cells (ATCC cat # HB-8065). One to fifty dilutions were used for each 
assay. The cell growth inhibition assay was performed for 72 hours. All other procedures 
were performed as described by Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499. 
The results of the assays using the three batches is displayed in Table 3. Based on these data, 

15 the Taiwan source would be selected as a standard herbal composition because of its low 

toxicity combined with its effectiveness in reducing secretion HbsAG (which is proportional to 
viral release) by more than haif. 
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Table 3. Bioassay of Xiao Chai Hu Tang (Sho-saiko-to). 



Source 


Cell Growth Inhibition (%) 


Hepatitis B Virus (secreted) 
% Inhibition 




HepG2 


2.2.15 


DNA 


HbsAG 


Singapore 


73 


100 


65 


38 


Korea 


13 " 


60 


20 


42 


Taiwan 


o ■ 


42 


0 


47 



The data presented in Tables 2 and 3 for the Taiwan herbal composition constitute the 
5 initial data for the standardized HBR Array for this herbal composition. Therefore, this data 
set would initially include the source of the herbal composition, the plant species and relative 
amounts of each the herbal composition, and two BioResponses {i.e., cell growth inhibition 
and hepatitis B virus secretion from infected cells). 

Using the procedures set forth in the schematic of Figure 1 and in Examples 1 and 3, 
10 additional data can be collected on plant-related data, markers and BioResponses for the 

standard herbal composition. This additional data is added to the initial standardized HBR 
Array to generate an expanded standardized HBR Array. Appropriate analyses of the resulting 
database can be conducted as set forth in the detailed description and the examples in order to 
ascertain the subset of variables which is most highly correlated or associated with the 
15 BioResponse of interest. Batch HBR Arrays may be determined using the methods depicted in 
Figure 2 and in the procedures of Examples 2 and 4. 

The resultant batch HBR Array can be compared to the standardized HBR Array so as 
to predict the BioResponse of the'batch herbal compositions. 
Example 12. Herbal Preparation 
20 The Standardized protocol : fdr the herbal extract preparation was as follows: The 

ingredients of herbal raw materials- with proper ratios were placed in a jacketed reactor and 
extracted with water at an elevated constant temperature with mixing. The solid was separated 
from the liquid with a 120-mesh screen. The resultant filtrate was collected and then 
concentrated by evaporating the water under reduced pressure. The concentrated liquor was 
25 spray dried at elevated temperature to yield granulated powder. This bulk substance was then 
formulated into the desired dosage form. 
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Example 13. Evaluation of Huang Qing Tang 

Huang Qing Tang (HQT) is an ancient Chinese botanical formula composed of four distinct 
herbs: Scutellariae (scute),. Glycyrrhizae (licorice), Paeonie lactiflora pallus (white peony 
root), and Fructus zizipho (date): (Table 4). This herbal formula has been long used in Asia to 
5 treat a variety of gastrointestinal ailments since 300 AD. 



Table 4. Herbal Ingredients of TCM Formula HQT 



Scientific Name 


Common Name 


Traditional Use 


Scutellariae Radix 


Scute Baical 
Skullcap root 


Used to reduce capillary permeability: to reduce 
inflammation: to treat enteritis and dysentery: increase the 
secretion of bile to treat jaundice: to relieve muscle 
spasms to treat coughing: to expel parasites. 


Glycyrrhizae 

Radix 

(Gancao) 


Licorice Ro'ot 


Used to moisten the lungs and stop coughs: to relax spasm 
and stop pain: to moderate the action of herbs; to reduce 
fire and release toxins. 


Fructus Ziziphi 


Date 


Has diuretic and strengthening effects. 


Paeonie lactiflora 
pallus radix 


White Peony 
Root . . • . . 


Used to suppress and soothe pain; to soothe ligaments and 
purify the blood. 



Biological and Enzyme Assays . ; 
10 Table 5. Batch Properties (HQT) 



Property 


Batch A 


Batch B 


Batch C 


Origin 


Taiwan, Sun-Ten 


Taiwan, Sun-Ten 


Taiwan, Sun-Ten 


Preparation method 


Standard.". 


Standard 


Boiled 30 min. 


Plant part 


Root ; 


.Root 





Briefly, one. gram of each batch of Huang Qing Tang (HQT) was added with 10 ml of 
water (1 mg/ml). The mtxtine was treated as outlined in Table 5. The supernatant was 
collected after centrifugation and filtered through a 0.22 um filter. Two cell types were used to 
15 test for biological effects of each batch of HQT: a) Jurkat T cells (ATCC cat #TTB-152) and b) 
HepG2 cells (ATCC cat # HB-8065). One to fifty dilutions were used for each assay. Frozen 
cells (10 7 /ml) were quickly thawed in a water bath at 37 °C. The cells were then diluted in 10 
ml of pre-warmed media (see Life Technologies, Inc., Catalogue and Reference Guide, 1998- 
1999, Cell Culture section) followed by centrifugation at 1500 rpm for 5 min. The supernatant 
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was then discarded and the cells were cultured in 100 ml media at 37 °C, 5% C0 2 . After 2 

days, the cells were counted (approximately 8 x 10 5 /ml, total 100 ml). 

Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by 

DNA quantitation (see Dong etal, Ptoc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, 
5 one gram of preparation was added with 10 ml of water. The mixture was boiled for 30 

minutes. The supernatant was collected after centrifugation and filtered through a 0.22 |j.m 

filter. HepG2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. 

Ace- see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. 

One to fifty dilutions were used for each assay. The cell growth inhibition assay was 
1 0 performed for 72 hours. All other procedures were performed as described by Dong et al, 

Proc Natl Acad Sci USA (1991) 88: 8495-8499. 

P-glucuronidase was assayed as HQT is known for its anti-diarrhea properties. 

Different HQT extracts were added to triplicate wells of a 96-well plate which contained 

O.lmM phenolphthalein glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of dialyzed p- 
15 glucuronidase (from E. ColU purchased from Sigma™) to a final volume of 80 ul. After 2 hr 

incubation at 37°C, the reactions were terminated with 200 ul of stopping solution which 

contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was monitored with a kinetic 

microplate reader at 540 nm. , 

The results of the assays using the three batches are displayed in Table 6. Based on 
20 these data, HQT sources A and B have relatively low toxicities combined with higher 

inhibitory activity relative to batch HQT C {i.e., approximately 5 fold greater toxicity toward 

HepG2 cells and 3.3 fold less inhibitory activity against P-glucuronidase than either HQT A or 

B, see Table 6). 

25 



30 
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Table 6. Biological Assay of Three Preparations of HOT* 





• ' E. Coli 
P^GIucuronidase 


HepG2 


Jurkat 


HBVJ 
DNA 


HQTA 


' •: 0.6 : 


1.50 


0.76 


None 


HQTB 


• -i ! 0.7 


1.6 


0.81 


ND 


HQTC 


• . 2.2 • '. 


0.32 


ND 


ND 


♦Values 

represen 

values. 

ND,not 

detennir 


t, % of Control 

ted. 









Evaluation o f HOT Effects on Protein Expression 

HepG2 cells (1 x 10 6 ) were seeded in 25 cm 2 flasks in 3.0 ml of RPMI-1640 medium 
(see Life Technologies,. Iric:,-Gatalqgue and Reference Guide, 1998-1999, Cell Culture section) 
24 hr before the drug addition. The cells were treated with or without herbal medicine, where 
the former is added at two final concentrations of 0.2 mg/ml or 4 mg/ml, respectively, and 
incubated at 37°C for 24 hours: The medium was removed and the cells were washed twice 
with cold PBS. The cells were harvested into 1 ml of PBS and centrifuged at 10,000 rpm for 2 
minutes, extracted on ice with a buffer containing 50 mM Tris-Cl (pH 7.5), 0.2 mM PMSF and 
10% glycerol, followed by three freeze-thaw cycles. Potassium chloride was added to the cell 
lysate at a final concentration of 0.15 M prior to centrifugation. The protein concentration was 
determined and the cell extract was electrophoresed according to the method of Laemmli 
(Nature (1970) 227:680-685). Western blots were performed by standard techniques known in 
the art, see for example Sambrook, et dl (1989). The antibodies used were directed to the 
following proteins: Topo l;: Stat' (20707); Cyclin Bl; MAPK (Ab2) and Nm 23 HI. 

Figure 4 demonstrates that the higher concentrations of HQT A or HQT B differentially 
effects the expression of cyclin B 1 . polypeptide. 

HPLC Analy sis f ; , \ \.';Z y 

The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere™ 
column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer 
(Perkin Elmer). The wavelengths for UV' detection were monitored at 280 nm and 340 nm. The 
mobile phase was pumped ; at 1 nil/rnin and consisted of Solvent A: H 2 0 and Solvent B: 20% 
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MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 
2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 
minutes; and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 
minutes. This was followed by the addition of 100 % solvent A for 5 minutes. The HPLC 
5 markers are baicalin and-baicaleih. -■ 
Mass Spectrometry . 

The herbal extract was analyzed by Mariner™ ESI-TOF Mass Spectrometry (MS) from 
PE Biosystems. Control tracings were generated using baicalein and baicalin, two known 
active ingredients in HQT. ' : 
1 0 HQT samples in water and acid treated batches were been analyzed by HPLC and Mass 

Spectrometry. While water treated HQT batches A and B had distinct HPLC and MS tracings, 
acid treated batches gave almost identical patterns (data not shown). 

Algorithm " ' ; ' V 

The data collected form part of the multidimensional analysis used to generate 
1 5 multivariant normal distribution sets as a means of determining a baseline correlation between 
biological activity and standard HQT chemical (HPLC and Mass Spec), and origin/growth 
characteristics. 

Example 14. Individual components 
A. Licorice . 

20 Evaluation of Glycyrrhizae Radix (licorice) 

Licorice is useful for moistening the lungs and reducing coughs, helps to relax spasm 
and pain. The properties of the licorice batches used in this example are presented in Table 7. 



Table 7. Batch Properties (licorice) 



Property 


Batch A 


Batch B 


Batch C 


Batch D 


Plant Name 


Glycyrrhizae 
Radix . 


Glycyrrhizae 
Radix ; 


Glycyrrhizae 
Radix 


Glycyrrhizae 
Radix 


Origin 


Inner Mongolia 


Inner Mongolia 


U.S., Kin Man 
Herb Center 


U.S., Kin Man 
Herb Center 


Preparation 
method 


Standard "'" 


Standard 


Boiled 30 min. 


WarmH,0, 30 
min. 


Plant part 


Root. • : 


Root 







Biologica l and Enzyme Assays 
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To assay the quality, of herbal sources, each herbal extract supernatant was assayed and 
the analysis was repeated three times. For a given sample to be assayed, 1 gram of herbal 
powder was dissolved in 10 ml of 80° C deionized water (neutral pH) in a polypropylene tube. 
The tube was then incubated as outlined in Table 7, then centrifuged to obtain the supernatant. 
5 Batches of licorice were tested against either HepG2 cells (ATCC cat # HB-8065) or Jurkat T 
cells (ATCC cat #TIB-152) or both. Cells were cultured for 24 hours as described above. 

Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by 
DNA quantitation (see Dong et al. Proc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, 
one gram of preparation was added with 10 ml of water. The mixture was boiled for 30 

10 minutes. The supernatant was collected after centrifugation and filtered through a 0.22 um 
filter.2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see 
Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. One to fifty 
dilutions were used for each' assay. The cell growth inhibition assay was performed for 72 
hours. All other procedures were performed as described by Dong et al., Proc Natl Acad Sci 

15 USA (1991) 88: 8495-8499): ; ' 

Again, p-glucuronidase was assayed. Different licorice extracts were added to triplicate 
wells of a 96- well plate which contained 0.1 mM phenolphthalein glucuronidate, 70 mM Tris- 
HC1 (pH 6.8) and 0.8 ng of dialyzed beta-glucuronidase (from E.Coli, purchased from Sigma) 
to a final volume of 80 ul and assayed as above. 

20 The results of the assays using the two batches is displayed in Table 8. Based on these 

data, licorice batch A was much more toxic to Jurkat cells than batches B (approximately 9 
fold) and a more effective inhibitor of p-glucuronidase (see Table 8). 



Table 8. Biological Assay of Four Preparations of Licorice* 




• :. "E. Cpli ' .. 
P-Glucuronidase 


HepG2 


Jurkat 


HBVJ 
DNA 


Licorice A 


■>■ 11. 


1.07 


0.41 


None 


Licorice B 


•/ :ND 


ND 


3.6 


ND 


Licorice C 


• 2.1 • 


ND 


ND 


ND 


Licorice D 


ND 


ND 


>2.0 


53.8 


♦Values re 
values. 
ND, not de 


jresent IC50 . %, % of Control 
tennined. * ,v. ( - • 
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Expression Assay 

In order to assay gene expression, Jurkat T cells were treated with herbal extract as 
follows: Jurkat cells (10 7 /ml)were quickly thawed in a water bath at 37 °C. The cells were 
then diluted in 10 ml of pre-warmed media (see Life Technologies, Inc., Catalogue and 
Reference Guide, 1998-1999, Cell Culture section) followed by centrifugation at 1500 rpm for 
5 mm. The supernatant was then discarded and the cells were cultured in 100 ml media at 37 
°C, 5% C0 2 - After 2 days ? the cells were counted (approximately 8 x 10 5 /ml, total 100 ml). 

The herbal extract solution was prepared as outlined above (e.g., 2 g of an herbal 
powder to obtain 20 ml of sterile solution (0. 1 g/ml). The cells were divided into 3 flasks at a 
density of 2.5 x 10 5 /mL 1-00 ml/each flask. Assays were carried out with control (no extract), 
and 10 ml of extract at 10 mg/ml, and 1 mg/ml. Again, toxicity results were used to deteimine 
the "high" and "low" concentrations for any given extract. After extract addition, cell cultures 
were incubated for 24 hours' under conditions as outlined above. The cells were counted and 
subsequently collected in 50 ml centrifuge tubes. The resulting cell pellet was treated with an 
RNA isolation means to extract mRNA (see, for example, Sambrook et al., 1989 at pages 7.3- 
7.39). ' . ' •'• " , 

Microarray '■■}; :'■ ' 

Microarray printing was carried out as follows: 

Human gene clones' were obtained from the IMAGE Consortium libraries through its 
distributors and comprise genes from various tissues. Most clones have been partially 
sequenced and the sequences were : available as expressed sequence tags in the dbEST database 
of GenBank. Clones were cultured and amplified using commercially available primers prior to 
application on nylon membranes (Chen et al., Genomics (1998) 51:313-324). Approximately 
10 ng of each amplified target was applied on a positively charged nylon membrane using a PC 
(personal computer) controlled arraying system. The arraying system allows high density 
spotting and is capable of depositing 31,000 spots on a piece of nylon membrane measuring 18 
by 27 mm using a 24-pin arraying tool. 

cDNA probe and Membrane Hybridization 

TwO microgram of each mRNA sample (mRNA was isolated as outlined above) was 
labeled with biotin and/or digdxigenin using random primed reverse transcription. The labeled 
samples were treated with alkali and the resulting labeled nucleic acids were precipitated prior 
to use in hybridization. Membrane hybridization and washing were carried out using the 
labeled probes as disclosed in Chen et al. (1998). To detect the spots on the membrane in dual 
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color mode (i.e., both biotin and chgoxigenin), P-galactosidase-conjugated streptavidin (Strept- 
Gal) and alkaline phosphatase-conjugated digoxigenin antibody (anti-Dig- AP) were employed. 
After color development, image digitization using an imaging means was employed (e.g., a 
flatbed scanner or digital .camera). Quantitative measurements were determined by computer 
5 analysis which uses a program that measures the integrated density of the primary color 
components of each spot, performs regression analysis of the integrated density data and 
locates statistical outliers as differentially expressed genes. 

- Gene expression data for samples L 2 and licorice (ST1 17) 

Extract 1, 2 and 6 corresponding to extract of Cordyceps sinensis, Poria cocos (ST 

10 027) and licorice, respectively, were assayed by the following method: Batches were evaluated 
for toxicity using Jurkat T cells. 

The extracts were prepared as outlined in Example 6. The cells were divided into 24 
well culture plates by adding 1 ml of Jurkat cells at a density of 5 x 10 5 /ml. Assays were 
carried out with control (no extract), and 5 concentrations of extracts as described (see Table 

15 9). The high and low concentrations for the cell culture assays were varied between 10 mg/ml 
and 0.05 mg/ml (i.e., mg dry weight of herbal extract per ml) depending on the toxicity of the 
extract to cells. For certain samples the toxicities at 10 mg/ml were such that "high" and "low" 
concentrations were adjusted downward, nevertheless, at least one order of magnitude between 
extremes was maintained. For example, for licorice (ST1 17) the "high" was 0.5 mg/ml and the 

20 "low" was 0.05 mg/ml (see. Table 9). After extract addition, cell cultures were incubated for 24 
hours under conditions as outlined in Example 6. The cells were counted and the resulting data 
tabulated to demonstrate extract toxicity. The resulting data is shown in Table 9. 

25 
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Table 9. Survival Cell Number at Different Concentration of Herbal Extract Solution 



nojcltfnil 

1 Cordyceps sinei 
mycelium 

2 ST027 

3 STO-44 

4 ST051 

5 ST093 

6 STJ17 

7 ST123 

8 STI28 

9 ST134 
. 10 ST237 



Note: original cell number is 5xl0Vml and the number to 10xl6 3 /ml after 24h incubation, 
describes all dead cells. 



Protocol: 

1 . Add 1 ml of 5xl0 5 /ml Jurkat cells into 24 well culture plates. 

2. Prepare 12 kinds of herbal extract solutions and sterilize. 

3. Test 5 concentration per sample. 10 mg/ml, 5 mg/ml, 2 mg/ml, 1 mg/ml, 0.5 mg/ml 

4. Culture the cells for 24 h in 37C with 5% C0 2 incubator. 



20 



In each analysis, 144X96 genes (i.e., 13,824 genes) were analyzed (data not shown) 
and about 100 genes showed significant differences in comparison with that of control (Table 
10). Some of the.genes were? up-regulated and others were down-regulated. The magnitude of 
the difference with the control sometimes varied depending on the relative amount of the 
herbal composition to which the particular cells were exposed. Numbers under CI (control 
treatment) and H or L (herbs) represent intensities of mRNA expressed after subtraction of 
background (Table 10). The gene designation is encoded in Array AD, which can be traced to 
a specific GenBank clone.' The level of expression was determined by H or L divided by C. 
Only a fraction of 13,824 genes in each herb treated samples showed significant changes, 
namely, up, down or unchanged (see 1 Table 10). 
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In this manner, we are able to correlate specific gene expression with the exposure of a 
cell to no, low (L) or high (H) amounts of an herbal composition. Many of the genes identified 
in this way code for proteins important in known metabolic or biochemical pathways. Many of 
5 these proteins have direct and indirect effects on certain physiological, morphological and 
psychological parameters. Thus, this method permits the association of a particular genetic 
fingerprint of an herbal composition with its array biological effects. Such associations can be 
used to profile or characterize an herbal composition for the purposes of Quality Control and 
Quality Assurance and evaluating pharmacological or toxicological properties. The role of 
1 0 primary and secondary herbs in an herbal formula can also be assessed by this approach. 

HPLC Analysis 

The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere™ 
column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer 
(Perkin Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The 

15 mobile phase was pumped at 1 ml/min and consisted of Solvent A: H 2 0 and Solvent B: 20% 
MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 
2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 
minutes; and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 
minutes. This was followed by the addition of 100 % solvent A for 5 minutes. The HPLC 

20 marker is glycyrrhizin. ' ' ' 1 

Al gorithm 

The data collected form part of the multidimensional analysis used to generate 
multivariant normal distribution sets as a means of detennining a baseline correlation between 
biological activity and standard licorice molecular, chemical (HPLC and Mass Spec), and 
25 origin/growth characteristics. 

B. Scute 

Evaluation of Radix Scutellariae (Scute) 

Scute has been found to be useful in reducing capillary permeability and inflammation 
It can also be used treat enteritis and dysentery, increases the secretion of bile to treat jaundice; 
30 to relieve muscle spasms; to treat coughing and to expel parasites. The properties of the scute 
batches used in this example are presented in Table 11. 
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Table 11. Batch Properties (Scute) 



Property 


Batch A 


Batch B 


Batch C 


Batch D 


Plant Name 


Scutellariae 
radix 


Scutellariae radix 


Scutellariae 


Scutellariae radix 


Origin 


Sanxi Province. 


U.S., Kin Man Herb 
Center 


U.S., Kin Man 
Herb Center 


U.S., Boston 


Preparation 
method 


Standard . 


Boiled, 30 min 


Warm H 2 0, 30 


Boiled , 2 hours 


Plant part 


Root 









Biological and Enzvme Assays 

Briefly, one gram of each preparation of scute extract was added with 10 ml of water (1 
5 mg/ml). The mixture was treated as Outlined in Table 11. The supernatant was collected after 
centrifugation and filtered through a 0.22 jam filter. Batches of scute were tested against either 
HepG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #TIB-152) or both. One to 
fifty dilutions were used for each assay. Cells were cultured for 24 hours as described above. 
Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by 

10 DNA quantitation (see Dong et al; Proc Natl Acad Sci ITS A (1991) 88: 8495-8499). Briefly, 
one gram of prep aration was added with 10 ml of water. The mixture was treated as outlined 
in Table 11. The supernatant was collected after centrifugation and filtered through a 0.22 um 
filter.2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see 
Ace et al. Proc Natl Acad Sci USA " (1987) 84: 1005-1009) were used in this assay. One to fifty 

15 dilutions were used for each assay. The cell growth inhibition assay was performed for 72 

hours. All other procedures were performed as described by Dong et al., Proc Natl Acad Sci 
USA (1991) 88: 8495-8499: : 

For P-glucuronidase, different scute extracts were added to triplicate wells of a 96-well 
plate which contained O.lmM phenolphthalein glucuronidate, 70 mM Tris-HCl (pH 6.8) and 

20 0.8 ng of dialyzed P-glucuronidase (from E. Coli, purchased from Sigma) to a final volume of 
80 pi. After 2 hr incubation at 37°C, the reactions were terminated with 200 pi of stopping 
solution which contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was 
monitored with a kinetic microplate reader at 540 mn. 
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The results of the assays using the three batches is displayed in Table 12. 
Table 12. Biological Assay of Four Preparations of Scute* 

"E. Coli HepG2 Jurkat HBVJ DNA 
P-GIucuronidase 

Scute A - .1.5. . 0.33 0.45 None 

Scute B ' . r.8 . : ' : ND ND ND 

Scute C . . 0.3 ND ND ND 

ScuteD ND ' 0.65 ND 27.5 

""Values represent J, % of Control 

IC50 values. 

ND, not determined. 

Evaluation of Scute Effects oh Protein Expression 

HepG2 cells (1 x 10*) were seeded in 25 cm 2 flasks in 3.0 ml of RPMI-1 640 medium 
(see Life Technologies^ Inc.j Catalogue and Reference Guide, 1998-1999, Cell Culture section) 
24 hr before the extract addition. The cells were treated with or without herbal medicine, 
where the former is added at two final concentrations of 0.2 mg/ml or 4 mg/ml, respectively, 
and incubated at 37°C for 24 hours. The medium was removed and the cells were washed 
twice with cold PBS. The cells were harvested into 1 ml of PBS and centrifuged at 10,000 rpm 
for 2 minutes, extracted on ice with a buffer containing 50 mM Tris-Cl (pH 7.5), 0.2 mM 
PMSF and 10% glycerol, followed by three freeze-thaw cycles. Potassium chloride was added 
to the cell lysate at a final concentration of 0.15 M prior to centrifugation. The protein 
concentration was detennuied and the cell extract was electrophoresed according to the method 
of Laemmli U.K. (Nature (1970) 227:680-685). Western blots were performed by standard 
techniques known in the art, see for example Sambrook, et al (1989). The antibodies used 
were directed to the following proteins: Topo I; Stat (20707); Cyclin Bl; MAPK (Ab2) and 
Nm23 HI. 

Figure 4 demonstrates that scute batches A and B do not differentially affect the 
expression of the polypeptides resolved on Western blots. 
HPLC Analysis r ■ •■• 

The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere™ 
column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer 
(Perkin Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The 
mobile phase was pumped at'l rhl/min and consisted of Solvent A: H 2 Q and Solvent B: 20% 
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MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 
2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 
minutes; and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 
minutes. This was followed by the addition of 100 % solvent A for 5 minutes. The HPLC 
markers are baicalin and baicalein. 

Scute batches in water and acid treated samples were analyzed by HPLC. Water and 
acid treated batches were virtually indistinguishable. 
- Algorithm 

The data collected form part of the multidimensional analysis used to generate 
multivariant normal distribution sets as a means of detenrnning a baseline correlation between 
biological activity and standard scute chemical (HPLC), and origin/growth characteristics. 

C. White Peony Root " , 

Evaluation of Paeonie lactiflora nallus radix fPeonyl 

Peony is used to suppress and soothe pain. It is also known to soothe ligaments and 
purify the blood. The properties of the peony batches used in this example are presented in 
Table 13. 



Table 13. Batch Properties (Peony) 



Property 


Batch A 


Batch B 


Plant Name 


Paeonie lactiflora pallus 


Paeonie lactiflora pallus 


Origin 


Anwey Province 


U.S., Boston 


Preparation metnod 


Standard • 


Boiled 2 hours. 


Plant part 


Root • 


Root 



Biological and Enzvme Assays 

Briefly, one gram of each preparation of scute extract was added with 10 ml of water (1 
mg/ro.l). The mixture was treated as outlined in Table 13. The supernatant was collected after 
centrimgation and filtered through a 0.22 um filter. Batches of peony were tested against either 
HepG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #TIB-152) or both. One to 
fifty dilutions were used for each assay. Cells were cultured for 24 hours as described above. 

Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by 
DNA quantitation (see Dong vt al. Proc Natl Acad Sci USA (1991) 88: 8495-8499). 
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Briefly, one gram of preparation was added with 10 ml of water. The mixture was 
treated as outlined in Table 13. The supernatant was collected after centrifugation and filtered 
through a 0.22 pm filter: 2.2.15 cells which secrete hepatitis B virons (kindly provided by 
Professor G. Ace; see Ace et al. Pror^tl Acad Sci USA (1987) 84: 1005-1009) were used in 
this assay. One to fifty dilutions were used for each assay. The cell growth inhibition assay 
was performed for 72 hours; All other procedures were performed as described by Dong et al, 
Pmr. Natl Acad Sci USA (1991) 88: 8495-8499. 

' Different peony extracts were added to triplicate wells of a 96-well plate which 
contained O.lmM phenolphthalein glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of 
dialyzedbeta-glucuronidase (from E.coli, purchased from Sigma) to a final volume of 80 pi. 
After 2 hr incubation at 37°C,-fhe reactions were terminated with 200 pi of stopping solution 
which contained 0.2 M Glycine iand 0.2 M NaCl (pH 10.4), and the OD was monitored with a 
kinetic microplate reader at; 540 nin. Results are shown in Table 14. 

Table 14. Biological Assay of Two Preparations of Peony* 





:E.CoU 

P^Giucuronidase 


. HepG2 


Jurkat 


HBVf 


Peony A 


2.8 


>1.5 


1.1 


None 


Peony B 


>2.5 


ND 


ND 


ND 


*Values represent IC 5 o 
values. 

ND, not determined. 


%,% of Control 









The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere column 
(5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin 
Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The mobile 
phase was pumped at 1 ml/min arid consisted of Solvent A: H 2 0 and Solvent B: 20% MeOH 
with the following gradient: i) the solvent was 100% solvent A for the first 5 minutes; 2) the 
solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; 
and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 minutes. 
This was followed by the addition of 100 % solvent A for 5 minutes. HPLC marker is 
paeoni florin. 

Peony batches were analyzed' by HPLC as shown in Figure 5. 
Algorithm 
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The data collected form part of the multidimensional analysis used to generate 
multivariant normal distribution sets as a means of deteraiining a baseline correlation between 
biological activity and standard peony chemical (HPLC), and origin/growth characteristics. 

D. Date 

5 Rvaluation of Ziziphi Fructus ( Date) 

Date has been used for diuretic properties and strengthening effects. The properties of 
the date batches used in this example are presented in Table 15. 



10 Table 15. Batch Properties (Date) 



Property 


Batch A 


Batch B 


Batch C 


Plant Name 


Ziziphi Fructus 


Ziziphi Fructus 


Ziziphi Fructus 


Origin 


Hebei Province. 


U.S., Kin Man Herb 
Center 


U.S., Kin Man Herb 
Center 


Preparation method 


Standard' .. . 


Boiled, 30 min 


Warm H 2 0, 30 min. 


Plant part 


Fruit ," . .. 







Biological and Enzyme Assays 

Briefly, one gram of each batch of scute extract was added with 10 ml of water (1 

mg/ml). The mixture was treated as. outlined in Table 15. The supernatant was collected after 
15 centrifugation and filtered through a. 0.22 pm filter. Batches of date were tested against either 

HepG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #TIB-152) or both. One to 

fifty dilutions were used for each assay. Cells were cultured for 24 hours as described above. 
Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by 

DNA quantitation (see Dong ei al, Proc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, 
20 one gram of preparation was added with 10 ml of water. The mixture was treated as outlined in 

Table 15. The supernatant was collected after centrifugation and filtered through a 0.22 pm 

filter. HepG2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. 

Ace; see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. 

One to fifty dilutions were us'ecl for each assay. The cell growth inhibition assay was 
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perfonned for 72 hours. 1 All other procedures were performed as described by Dong et ah, 
Pror. Natl Acad Sci USA (1991) 88: 8495-8499. 

Different peony extracts were added to triplicate wells of a 96-well plate which 
contained O.lmM phenolphthalein.glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of 
5 dialyzed beta-glucuronidase (from E. Coli, purchased from Sigma) to a final volume of 80 ul. 
After 2 hr incubation at 37°G, the reactions were terrninated with 200 ul of stopping solution 
which contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was monitored with a 
kinetic microplate reader at 540 nm. Results are shown in Table 16. 



Table 16. Biological Assay of Three Preparations of Date* 





KColi 
P-Glueuronidase ' 


HepG2 


Jurkat 


HBVJ DNA 


Date A 


' 1,2 ■ 


1.5 -. 


5.1 


None 


DateB 


•ND 


>2.0 


ND 


52.3 


DateC 


2.5 


ND 


ND 


ND 


♦Values represent 

IC50 values. 

ND, not detennined. 


%, % of Control 









10 HPLC Analysis • - .. 

The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere column 

(5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin 

Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The mobile 

phase was pumped at 1 rnl/min and consist of Solvent A: H 2 0 and Solvent B: 20% MeOH with 
15 the following gradient: 1) the solvent was 1 00% solvent A for the first 5 minutes; 2) the 

solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; 

and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 minutes. 

This was followed by the addition of 100 % solvent A for 5 minutes. HPLC markers for date 

are chelidonic acid and cAMP. 
20 Date batches samples were analyzed by HPLC as shown in Figure 6. 

Algorithm 

The data collected form part of the multidimensional analysis used to generate 
multivariant normal distribution sets as a means of determining a baseline correlation between 
biological activity and standard peony chemical (HPLC), and origin/growth characteristics. 
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Example 15. Characterizing Herbal Medicines bv Nucleic Acid Microarray 
Analysis. 
Introduction. 

The rapid development of nucleic acid microarray technology has led to an explosion 
5 of gene expression data (Lander, 1999, Duggan et al., 1999). Four characteristics of the gene 
expression account for the great value of using nucleic acid microarrays to study the gene 
expression profiles, (i) Nucleic acid microarray makes it easier to measure the transcripts of 
thousands of genes at once, (ii) Close association between the function of a gene product and 
its expression pattern makes gene function predictable, (iii) Cells respond to the micro- 

10 environmental changes by changing the expression level of specific genes, (iv) The sets of 
genes expressed in a cell determine, what the cell is derived of, what biochemical and 
regulatory systems are involved, and so on (Brown and Botstein, 1999). By using a microarray 
system, the above features can be studied in an ensemble manner. 

The expression of any. desired number of genes can be detected using the nucleic acid 

1 5 microarray technology. Tor; example, up to about 20,000 genes may be placed on a single 
array. We have developed a nucleic acid microarray with colorimetric detection system 
(Microarray/CD) (Chen ei al., 1998). Gene expression profiles of different cell lines were 
studied using microarray filter membranes (2.7 cm x 1.8 cm) with about 10,000 cDNA 
representing approximately 10,000 distinct human transcripts. The sensitivity and detection 

20 limits of the rmcroarray/CD ; system- have been characterized and are comparable to the system 
with radioactive detection or the system with laser induced fluorescence detection (Bertucci et 
al., 1999). "».;'■■■ 

As previously described; cellular gene expression profiles portray the origin, the 
present differentiation of the cell, and the cellular responses to external stimulants. In other 

25 words, the gene expression profiles reveal the state of the cell and microarray is a perfect tool 
for the rendering purpose. In the present studies, we apply the microarray/CD system to 
characterize cellular responses to external stimulants, in this case, the Chinese herbal 
medicines. Conversely, we also based on the stimulated gene expression profiles to classify 
different herbal medicines: ' 

30 Figure 7 is a flowchart depicting a general method that may be used for establishing an 

expression response data set for cells treated with an herbal composition. The method 
comprises the steps of: 
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(a) Determine the IC 5 6 concentration of ian herbal composition by incubating various 
concentrations of the herbal medicine in mammalian cell cultures and identify the 
concentration that leaves 50% of survival cells after a predetermined time. 

(b) Incubate the mammalian cell cultures with herbal extracts of various fractions of 
IC50 concentrations. 

(c) Harvest and count the cultured cells after a predetermined culture time. 

(d) Immediately lyse the cells after they are removed from the incubator and extract 
mRNA from cell lysate. 

(e) Label the mRNA by reverse transcription reaction to turn mRNA into labeled cDNA. 

(f) Mix the labeled cDNA with control cDNA of plant origin and perform hybridization 
to a microarray of mammalian gene probes. 

(g) Measure expression level of genes by analyzing digitized images of the microarray 
hybridization results. 

(h) Perform data pre-processing to select data for statistical analysis. 

(i) Acquired expression data generated by microarray experiments of an herbal 
composition with various concentrations. 

(j) Data pre-processing to select the genes with statistical significance in cells treated 

with different concentrations of the herbal medicine, 
(k) Categorize expression profiles into clusters by statistical methods such as the self- 

organizing-map algorithm. 
(1) Deduce the characteristic expression profiles for the herbal medicine based on the 
expression profile clusters. 
Figure 8 is a flowchart demorWtrating how data sets of expression data for various batches of 
the herbal composition are integrated to make an expression profile database for the particular 
herbal composition. The expression profile database then becomes part of the HBR Array. 

HBR Arrays containing expression profiles may also be used to identify an unknown 
herbal composition. Figure 9 is a flowchart depicting a general method for identifying an 
unknown herbal composition, the method comprising the steps of: 

(a) Construct an HBR Array c6ntaining characteristic expression profiles for an herbal 
medicine or a collection of expression profiles of various herbal medicines by the 
aforementioned steps. 

(b) Obtain the characteristic expression profile data set of the unknown herbal 
composition. 
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(c) Compare the HBR Array containing the characteristic expression profile induced by 
the said unknown herbal composition with a standardized HBR Array containing 
expression data by algorithms such as the Hannning distance algorithm. 

(d) Score possible ahgriments to identify the most probable herbal composition whose 
5 characteristic expression profiles are archived in the said HBR Array. 

Scoring possible alignments of HBR Arrays containing expression profiles may be performed 
using hierarchical cluster analysis of the Hamming distance matrix. Use of hierarchical cluster 
analysis for the Hamrning distance matrix is well known in the art. 

The gene expression profiles may also be incorporated into the standardized HBR Array. 

10 As has been already discussed, the standardized HBR Array containing such gene expression 
profiles induced by an herbal composition can be used for studying the pharmacological 
mechanisms of the herbal composition, for discovering new application of the herbal 
composition, and for designing optimized formulation of a complex herbal preparation. As can 
be seen from the flowchart of Figure 10, the method may be generally outlined as comprising 

15 the steps of: 

(a) Construct a data set containing the characteristic gene expression profiles for an 
herbal composition. • . ' 

(b) Score each gene by the consistency of its expression profiles in the data set using 
known statistical parameters, such as the coefficient of variation. 

20 (c) Based on the statistical scoring, gene expression profiles for an herbal medicine are 

selected to be incorporated into the standardized HBR Array. 
HBR Arrays containing gene expression profiles may also be used to identify signature 
gene expression profiles induced 'by. individual chemical constituents in an herbal composition 
consisting of complex chemical constituents, as outlined in the flowchart of Figure 11. The 
25 method comprises the steps of: 

(a) Construct a HBR Array containing characteristic gene expression profiles for an 
herbal composition by the aforementioned steps. 

(b) Determine me composition of chemical constituents in an herbal medicine by high 
performance liquid chromatography (HPLC) or liquid chromatography mass 

30 spectrometry (LC-lvlASS). 

(c) Repeat the step (b) for various batches of herbal medicine preparations. 

(d) Score the correlationcoerficients between the expression levels of each gene with the 
amount of individual chemical constituent in an herbal preparation. 
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(e) The signature gene expression profiles for individual chemical constituent are selected 
with a Pearson correlation coefficient exceeding 0.99 or smaller than -0.99. 

Any herbal composition can then be characterized through the use of gene expression 
profiles generated through the usfe of nucleic acid microarrays. Moreover, one can choose any 
number of genes that are differentially expressed to be included in the data set represented the 
gene expression profiles. For example, one may choose about 10 genes, about 100 genes, 
about 500 genes, about 1000 genes, about 1500 genes, abut 2000 genes, about 2500 genes or 
more, or any number in between. 

The prescription of the Chinese herbal medicine Scute and Licorice combination 
(Huang Chin Tang) stops diarrhea, relieves spasms and clears fever. The ingredients of Huang 
Chin Tang are Scute, Peony, Licorice and Jujube. This recipe has been used for more than 
1000 years but the chemical and biomedical studies on the prescription have not been carried 
out until recent decades. In tins study we used the nucleic acid microarray technology to study 
the gene expression profiles' of: herbal medicines treated cells. Our aims are to demonstrate the 
feasibility of using the microarray/CD system for classification of different herbal 
compositions or different preparations and to find the predictor genes (marker genes) for the 
Huang Chin Tang prescription. The long-term goals are to find the correlation of the 
biochemical ingredients in each herbal composition with the gene expression profiles of 
various treated cells and to decipher the molecular pharmacological mechanisms of the 
Chinese herbal medicines ih a rational fashion. ' 

Materials and Methods. 

1. Development of a cell banking system. 
Purpose: Microarray System is a sensitive detection method to monitor gene expression 

patterns of cells: It is necessary to build a Cell Banking System with a Master 

Cell Bank' (MCB) and a Working Cell Bank (WCB) to minimize ceU variability 

for herbal medicine testing. 
Scope: The Cell Bank System is used for all types of cells in microarray studies. 

Apparatus: CO z Air-Jacketed Incubator (NUAIRE™ DH autoflow) 

Centrifuge (KUBOTA 2100) 

Freezing vial (Corning Costar, Cat. #430659) 

Tissue culture flask 750 ml (Falcon, Cat. #3045) 

Tissue culture dish 1 50x25 mm (Falcon, Cat. #3025) 
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Cell: Jurkat T cell from Dr. Alexandra Ho 

Reagents: RPMI Medium 1640 (GIBCO BRL, Cat. #31800-014) 
Dimethyl Sulphoxide (DMSO) (Sigma, Cat. #D-2650) 
Fetal Bovine Serum (HyClone, Cat. #SH3 0070.03, Lot#AGL7258) 
2-mercaptoefhanol (GIBCO BRL, Cat. #21985-023, 5xl0' 2 M) 

Media: I. Culture medium: 90% RPMI + 10% fetal bovine serum + 2- 

mercaptoethanol (5x1 0" 5 M) 
H. Freezing medium: 90% RPMI 1640+1 0% DMSO 

Procedure: 

A. Master Cell Bank 

1 . Follow the standard sterile procedure of cell culture. 

2. Seed Jurkat T cells in culture medium in flask at 37°C with 5% C0 2 incubator. 

3. After incubation for two days, count the cell number and spread the cells to two flasks. 
Note. Cell density is kept about 5xl0 4 - 2xl0 6 per ml. 

4. Culture and count the cell number until the number reaches 2xl0 7 . 

5. Collect the cells in 50 ml centrifuge tubes and spin at 1300 rpm (300xg) for 5 min. 

6. Discard the supernatant and re-suspend the cell pellet with chilled freezing media. The 
cell number in each vial is about lxl 0 6 per ml. 

7. Slowly freeze the cells by the following temperature profile: -20°C for 2 hr, 

-80°C for 24 hr and then place the cells in liquid N 2 storage. Store a total of 20 frozen 
vials in MCB. 

B. Working Cell Bank ".' 

1. Retrieve one vial of cell from MCB in liquid N 2 tank and quickly thaw at 37°C water 
bath. 

2. Transfer the cells into , 10 ml of warm culture medium. 

3. Spin down the cells at 1300 rpm (300*g) for 5 min. Discard the supernatant. Culture 
the cells with 20 ml medium in a flask. " 

4. Sub-culture the cells to 2 flasks. 

5. Seed 5xl0 7 cells with 500 ml culture medium in each flask with stirring for a total of 2 
flasks. Culture the 1 cells for 2 days. 

6. Culture until the cell density reaches lxl0 6 /ml and a total volume of 1 L. 

7. Prepare freezing media by adding 100 ml of fetal bovine serum and 10 ml of DMSO. 
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Centrifuge, discard the supernatant and re-suspend the cells to 1 10 ml of freezing 
medium. 

Dispense 1 ml to every freezing vial (10 million cells per vial) for a total of 100 vials. 
Slowly freeze the ceils by the temperature profile described above. 



2. Determination of growth inhibition concentration of herbal extract in cell 
cultures. 

Most drugs are toxic to cells. This experiment is designed to examine the 
toxicity of herbal extracts in Jurkat T cells and to deterrnine the growth 
inhibition concentration of herbal extracts that keeps the cells alive. 
This assay can be used in all kinds of herbal extracts to examine the toxicity. 
C0 2 Air-Jacketed Incubator (NUAIRE™ DH autoflow) 
Counting chamber (Hemacytometer, Reichert, USA) 
Microscope (Zeiss, Axiovert 1 00) 
Jurkat T cell 

RPMI Medium 1640 (GIBCO BRL, Cat. #31800-014) 
Fetal Bovine Serum(Hy Clone, Cat. #SH30070.03, Lot #AGL7258) 
2-mercaptoethanol (GIBCO BRL, Cat. #21985-023, 5xl0" 2 M) 
Culture media: 90% RPMI + 10% fetal bovine serum + 2-mercaptoethanol 
(5xlO" 5 M) 

Disposable sterile syringe filters (0.2 m, Corning, Cat. #21052-25) 
Herbal extracts: 



Purpose: 



Scope: 
Apparatus: 



Cell: 

Reagents: 



1. 


Cordyceps Sinensis Mycehum 


2. 


ST 024 




3. 


ST 044 




4. 


ST 051 




5. 


ST 093 




6. 


ST 117 




7. 


ST 123 




8. 


ST 128 




9. 


ST 134 




10. 


ST 237 




11. 


PHY906-303503: Complex mix composed of 4, 6,7, 10 
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12. PHY906-284003: Complex mix composed of 4, 6, 7, 10 
Procedure: 

A. Herbal Extract Preparation 

1. Dissolve 1 gram of herbal powder in 10 ml of 80 °C deionized water (neutral pH) in a 
polypropylene tube. 

Incubate the tube at : 80 °C water bath for 30 minutes with gentle shaking then 
centrifuge at 4000 rpm (1500*g) for 5 min to obtain the supernatant. 
Centrifuge at 11000 rpm (14000xg) for 10 min to collect the supernatant. 
Using the disposable sterile syringe filter to filter the supernatant. 



2. 



7. 



B. Cell Survival Test 

1 . Culture Jurkat T cells as described above. 

2. Dispense 1 ml of 5xi0 5 /ml cells per well to 24-well culture plates. 

3. Prepare 12 kinds of herbal extract solutions. The extract solutions must be freshly 
prepared and used immediately. 

4. Add 100, 50, 20, 10, 5 pi of each herbal extract solution into the 24 well culture plates 
to get the five different concentrations: 10, 5, 2, 1, 0.5 mg/ml. 

5. Culture the cells for 24 h at 37°C in an incubator filled with 5% C0 2 . 

6. Count the number of cells per well. Mix 10 pi of cell solution with 10 pi of Trypan 
blue dye and load into cell counting chamber. 
Count the four major square areas to calculate the cell number. 

(number of cells in 4 areas)/4 x 10 4 x dilution factor = number of cells per ml 
3. Profiling gene expression patterns of Jurkat T cells treated with herbal 
extracts. 

Profile the gene expression patterns of Jurkat T cells treated with herbal 
extracts. A high-density nucleic acid microarray with colorimetric detection 
system is used. 

Heat block (Boekel, Model 1 10002) 
Spectrophotometer (Beckman, DV640) 
Centrh%ge(KUBOTA 1910) 
Water bath (SLM AMINCO, Model 800) 
Hybridization incubator (YIH DER OH-800) 
Heat sealer(TISH-300, TEW) 
Reagents: RNAzol™. B (Tel-Test, Cat. #CS-104) 



Purpose: 



Apparatus: 
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Oligotex iriRNA Midi Kit (Qiagen, Hilden, Germany) 
Hybridization Bags (GIBCO BRL, Cat. #18278-010) 
EasiSeal (Hybaid, Cat. #HBOSSSEZlE) 
Glass slides (Matsunami, S2214, Japan) 

Aerosol Resistant Tips (ART tip) (molecular BlO-Products, Cat. #2139) 

Random hexartier primer (GIBCO BRL, Cat. #48190-01 1) 

Reverse transcriptase and 5* buffer (GIBCO BRL, Cat. #18064-014) 

RNase inhibitor (GIBCO BRL, Cat. #10777-019) 

Biotin-16-dUTP (Boehringer Mannheim, Cat. #1093070) 

Dig-il-dUTP (Boehringer Mannheim, Cat. #1558706,) 

Blocking powder for hybridization (Boehringer Mannheim, Cat. #1 096 1 76) 

Bovine serum albumin (Sigma, Cat. #A2153) 

20X SSC (Amresco, Cat.#0918S-2-20XPTM5L) 

SDS (Merck, Cai.#113760) 

Dextran Sulfate (Sigma, Cat. #D6001) 

Streptavidin-p-galactosidase (GIBCO BRL, Cat. #19536-010) 
Anu-digoxigenih-AP Fab fragments (Boehringer Mannheim, Cat. #1093274,) 
X-gal (GIBCO BRL, Cat: #15520-018) 
Maleic acid (Sigma, Cat. #M1125) 
N-lauroylsarcosih (Sigma, Cat. #L5777) 
Fast red TR/AS-MX substrate kit (PIERCE, Cat. #34034) 
Polyethylene glycol (Sigma, Cat #P2 1 39) 
Reagent Preparation: 

lx hybridization buffer (4X SSC, 0.1% N-lauroylsarcosine, 0.02% SDS, 1% BM 



blocking reagent) 1 .. . ... .. 

20x SSC 16 m 

1% N-lauroylsarcosine 8 ml 

10% SDS 160 \x\ 

BM blocking powder .' 0.8 g 

H z O . " ... - 51" ml 



total 80 ml 
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Heat to 65°C to dissolve the powder then store at -20°C. 

50% PEG-8000 (Polyethylene Glycol) 
PEG-8000 : 10 g 

5 H 2 0 up to 20 ml 

Heat to 65°C to dissolve then autoclave. Aliquot and store at -20°C. 

- 10* TBS (100 mM Tris, 1.5M NaCl, pH 7.4) 
Trisbase ; 12.1 g : 

10 NaGl 87.6 g 

H 2 0 up to 1000 ml • 

120 mM X-gal 

X-gal 100 mg ' ; v 

15 DMF 2 mi ' : 

Store at -20°C. 

X-gal Substrate Buffer (ImM MgCl 2 , 3 mM K 3 Fe(CN) 6 , 3 mM K4Fe(CN) 6 in IX TBS 
buffer) 500 ml 

20 10XTBS/pH7.4 50 ml 

Potassium Ferrocyamde 633.5 mg 
Potassium Ferri cyanide 493.9 mg 
MgCl 2 101.6 mg 

25 Filter and store at -20°C. 

BM Blocking Dilution Buffer/pH7.5 (0.1 M maleic acid, 0.15 M NaCl) 

IMMaleicAcid ." 100 ml 

5M NaCl ' ' 30 ml 

30 Solid NaOH 7.5 g 

H 2 0 up to 1000 ml 

10% Blocking Reagent 100 ml 
Blocking Powder 10 g 
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Blockang Dilution Buffer (no tween 20) 
Heat to 70°C then autoclave: Store at 4°C. 

20% Dextran Sulfate 

Dextran Sulfate "2g 
H 2 0 Up to 8 ml 

Autoclave then store at -20°C. 

DEPC-treated water 800 ml 

Diethyl Pyrocarbonate (store at 4°C) 400 ul 

H 2 0 800ml 
Put in 37°C shaking water bath for 4 hrs and then in 37°C warm room for overnight. Autoclave 
the solution for 45 minutes or 25 minutes each for two times. 
Procedure: 

A. Preparation of Jurkai'T Cells 

1 . Thaw one vial of Jurkat T cells. Transfer to 1 0 ml growth medium. The number of vials 
to thaw depends oh the number of test to perform. In general, one vial of cells is needed 
for 2 herbal extract tests. 

2. Re-suspend the cells in 50 ml of culture medium. 

3. Incubate the cells for one day. Add 150 ml of culture medium and divide into two 
flasks with 100 ml each. 

4. Culture the cells for 3 days. 

5. Change the medium and distribute 4 flasks with 100 ml of medium each. 

6. Culture the cells for 2 days. 

7. Count the cell number. Collect the cells and spin down. Re-suspend the cell pellet with 
culture medium to 5xl.0 ? cells per ml and 100 ml per flask. 

8 . Culture the cells for 3 hrs before adding the herbal extract. 
B. Herbal extract treatment ' 

1 . Prepare the herbal extract (see herbal extract preparation). 

2. According to the growth inhibition concentration of herbal extract determined by the 
cell survival experiments, calculate the "50% growth inhibition concentration of each 
herbal extract. : 
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3. Define the 50% growth inhibition concentration as H. Treat the cell with serial dilution 
of herbal extracts with the following concentrations: H, H/2.5, H/5, H/10, H/20. 

4. Culture the cell for 24. hrs. 

5. Collect cells and count the number of cells. Centrifuge to obtain cell pellet. 
5 6. Wash the cell pellet with 1 *PBS once. 

7. Discard the supernatant. The cell pellet is ready for total RNA isolation. 

C. Isolation of total RNA 

1. - Add 1 ml of RNAzol™ B per 10 7 cells. Homogenize the cell pellet but do not vortex. 

2. Add 0.1 ml chloroform per ml of homogenate, cover the samples tightly, shake 
10 vigorously for 1 min .(do not vortex). Place on ice for 15 min. 

3 . Centrifuge at 1 2000. rpm (1 3500xg) at 4°C for 1 5 min. 

4. After centrifugation, uie homogenate develops two phases: a lower blue phenol- 
chloroform phase and a colorless upper aqueous phase. DNA and proteins are in the 
interphase and the organic phase. Transfer the aqueous phase to a new tube, add an 

15 equal volume of isopropanoiand store the samples at -80°C. Note. The range of 

isopropanol addition is 1 from 0.7 to 1 volume of the aqueous phase solution. 

5. Keep the samples at -80°C until use. Let the sample completely thaw before 
centrifuging and mix 2 to 3 times by inverting the tube. Centrifuge samples for 15 min 
at 13000 rpm(15000x < £r). 

20 6. Remove the supernatant and wash the RNA pellet once with 1 ml of 75% ethanol. 
Centrifuge for 3 niih at- 13000 rpm (15000x^) and at 4°C. 

7. Discard me -supernatant. Dry the pellet under vacuum for 1 min. Note. Do not let the 
RNA pellet dry completely. Tt will greatly decrease its solubility. 

8 . Dissolve the RNA pellet in' 50-1 00 ul of diethylpyrocarbonate (DEPC) - treated water 
25 by pipetting. Note. If the pellet is hard to dissolve, incubating the pellet for 10 - 15 min 

at 60°C mayhefpY '■'■'< ' * 

9. Measure absorbance af 260 nm (A260) and 280nm (A280) with a spectrophotometer. 
Concentration analysis: OD 2 6o X 40 ng/ul x dilution factor = total RNA (ng/p.1). 

D. Isolation of polv-A+ rtiRNA from total RNA 

30 1 . Determine the amount of starting RNA and the appropriate volume of Buffer OBB and 
Oligotex Suspension solution to be added in the RNA solution according to the Table 
17. 
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Table 17. Buffer amounts for Oligotex mRNA Spin-Column Protocol. 



Total RNA 


Add RNase-free 


Buffer OBB 


Oligotex 


Prep size 




water to: (pi) 


0*1) 


Suspension (pi) 




20 ug 


100 


100 


6 


Mini 


0.25 mg 


250 - 


250 


15 


Mini 


0.25-0.5 mg 


500 


500 


30 


Midi 


0.5-0.75 mg 


500 


500 


45 


Midi 


0.75-1 mg 


500 


500 


55 


Midi 



The following procedures are based on using 500 ug total RNA as an example. 

2. ' Add 500 pi of 2 x Binding Buffer and 30 pi Oligotex Suspension to the total RNA 
5 sample. Mix the contents thoroughly by flipping the tube. 

3. Incubate the sample for 10 min at 70°C. 

4. Incubate for 20 min at room temperature. 

5. Centrifuge for 2 min at maximum speed (14000 to 18000*g) and aspirate the 
supernatant. 

10 6. Re-suspend the pellet in 400 pi of Wash Buffer O W2 and transfer onto spin column 
and centrifuge the spin column for 1 min. 

7. Wash with 400 ul of OW2 and centrifuge as above. 

8. Add 20 ul of preheated (70°C) Elution Buffer onto the column and re-suspend the 
resin. Close the microcentrifuge tube. 

15 9. Put the spin column with 1 .5 ml microcentrifuge tube at 70°C for 3 min. 

10. Centrifuge the column at maximum speed for 2 min at room temperature. 

11. Elute again, (repeat step 8 to 10 to get better yield) 
E. cDNA Labeling ' 

1. Mix 2pg mRNA, 1 pi of control plants mRNA for single color label (Hat22: lxlO 9 , 
20 Rbcl: 5xl0 8 , Ga4: lxlO 8 , Rca: 5*10 7 , Asal: lxlO 7 , Atps: 5*10 6 molecule/pl), 6 pi of 

50 mM random hexamer and DEPC-H 2 0 to 28.88 pi final volume. For dual-color 
mode, use 2 pg of mRNA each in Biotin or Dig labeling and individual addition of 
control plants mRNA: 1. biotin labeling: Hat22: 1x10 s , Rbcl: 5xl0 7 , Ga4: 2xl0 7 , Rca: 
lxlO 7 , Asal: lxlO 7 , Atps: TxlO 7 , Hat4: lxl0 7 /pl. 2. Dig labeling: Hat22: lxlO 7 , Rbcl: 
25 lxlO 7 , Ga4: lxlO 7 , Rca: lx.10 7 , Asal: lxlO 8 , Atps: 5xl0 7 , Hat4: 2xl0 7 /pl. 

2. Denature for 10 min at 70°C, then chill quickly in ice for 5 min. 

3. Add 10 pi of 5x first strand buffer, 5 pi of 0.1 M DTT, 1 pi of 25 mM dATP, dCTP, 
dGTP mixture, 1 pi of 2 mM dTTP , 2 pi of 1 mM biotin- 16-dUTP, or Dig-1 1-dUTP (1 
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mM), 0.63 ul of 40 U/ul RNAsin and 1.5 ul of Superscript II (reverse transcriptase, 
GIBCO BRL)(200 U/ul). 

4. Mix well and incubate for 1 0 min at 25°C, then for 90 min at 42°C. 

5 . Stop the reaction for 5 min at 94°C. 

5 6. Add 5.5 ul of 3 M NaOH for 30 min at 50°C. 

7. Add 5.5 ul of 3 M. CH 3 COOH for 30 min at 50°C. 

8. Precipitate the labeled cDNA by adding 34 ul of water, 50 ul of 7.5 M ammonia 
acetate, 10 fig of linear polyacrylamide as carrier and 380 ul of absolute alcohol. 

9. Incubate the sample for 30 min at -80°C. Centrifuge at 13000 rpm for 15 min. 
10 10. Wash the pellet with 1 ml of 70% ethanol and centrifuge at 1 3 OOOrpm for 5 min. 

11. Dissolve the pellet in 36 ul of autoclayed H 2 0. For dual color, combine two labeled 

cDNA together. 
F. Array Hybridization 

1. The filter membrane carrying the 9600 EST PCR products is pre-hybridized in 5 ml of 
15 lx hybridization buffer (4X SSC, 0.1% N-lauroylsarcosine, 0.02% SDS, 1% BM 

blocking reagent (Boehringer Mannheim)), and 50 fxg/ml salmon sperm DNA (GIBCO 
BRL) at 63°C for 1.5 hours. Note. You can prepare 80 ml of lx hybridization buffer 
and store it at -20°C- Thaw the buffer at 60°C before use. 

2. Stick one side of adhesive EasiSeal® square to a clean glass slide and place the pre- 
20 hybridized membrane in the center of the square with the spots facing up. 

3. Mix the probe with 2 ul of poly-d(A)10 (10 ug/ul) and 2 ul of human Cot-1 DNA (10 
ug/ul) (GIBCO BRL) and 40 ul of 2* hybridization buffer to 80 ul final volume. 

4. Denature the probe mixture at 95 °C for 5 min and then cool on ice. 

5. Seal the filter membrane with the probe solution in the hybridization bag. 
25 6. Incubate at 95°C for 5 min and then at 63°C for 12-16 h (overnight). 

7. Wash the filter membrane twice with 5 ml of 2x SSC, 0.1% SDS for 5 min at room 
temperature. ; 

8. Wash three times for 15 min each with 5 ml of O.lx SSC, 0.1% SDS at 63°C. 

9. Block the filter membrane with 5 ml of 1% BM blocking reagent containing 2% 
30 dextran sulfate at room temperature for 1 h. 

10. Incubate with 5 ml . mixture containing 700x diluted Streptavidin-p-galactosidase 
(1.38U/ml, enzyme activity)(GIBCO BRL), lOOOOx diluted anti-Digoxigenin-alkaline 



WO 01/66803 PCT/US01/07608 

-80- 

phosphatase (0.075U/mL enzyme activity)(Boehringer Mannheim), 4% polyethylene 
glycol 8000 (Sigma), and 03% BSA in lx TBS buffer at room temperature for 2 hours. 
Note. This formula is for dual-color mode. For single color mode, anti-Dig- AP is not 
needed and the incubation time can be reduced to 1 hour. 
5 11. Wash with 1 * TBS buffer three times for 5 minutes each. 

12. Freshly prepare X-gal substrate solution (1.2 mM X-gaL ImM MgCl 2 , 3 mM 
K3Fe(CN) 6 , 3 mM ^e(CN) 6 in l x TBS buffer) by mixing 50 ul of 120 mM X-Gal 
" . and 5 ml of X-Gal substrate buffer. Immerse the filter membrane in the X-gal substrate 
solution for 45 min at 37°C with gentle shaking. 
10 13. Wash with lx TBS. 

14. Dual color development: stain the membrane with 5 ml of Fast red TR/naphthol AS- 
MX substrate (Pierce- Rockford, EL) at room temperature for 30 minutes with gentle 
shaking. : « . :; \ ■:. 

15. Wash with deionized water. Stop the reaction with lx PBS containing 20 mM EDTA 
15 for 20 min. 

16. Air dry the filter membrane. 
Results. 

1. Determine the growth inhibition concentrations of herbal extracts in cell 
cultures. 

20 Each herbal extract has different cellular toxicity, thus it is necessary to determine the 

growth inhibition concentration of every herbal medicine before treating the cells. Five serial 
dilutions of herbal extract' (lO/S, 2,1, 0.5 mg/ml) were added to 5xl0 5 /ml cultured cells and 
incubated for 24 hours in an incubator at 37°C with 5% C0 2 . The numbers of survival cells at 
different concentrations of herbal extracts are shown in Table. 18. 

25 



30 



35 
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Table 18. Survival Cell Number at Different Concentrations of Herbal Extracts 















experimental 
concentration 




no.xlu /ml 


mg/ml 


5 

mg/ml 


2 

mg/ml 


1 

mg/ml 


0.5 
mg/ml 


No 
drug 


High cone 
(mg/ml) 


Low cone 
(mg/ml) 




Cordyceps sinensis 
mycelium 


8.4' 


11.9 


11.5 


9.2 


9.0 


12.4 


10 


1 


2 


ST024 


4-2 












10 


1 


3 


ST044 




5.9 


8.4 


9.9 


9.4 




5 


0.5 


4 


ST051 








1.7 


5.4 




0.5 


0.05 


5 


ST093 






1.9 


3.8 


4.4 




0.5 


0.05 


6 


ST117 




1.6 


3.6 


4.6 


5.8 




0.5 


0.05 


7 


ST123 


3,4 


6.4 


8 


9.3 


7.8 




, 5 


0.5 


8 


ST128 


■ 3.5 


7.7 


7.9 


7.7 


8.3 




5 


0.5 


9 


ST134 


2.9 


6.1 


11.2 


9.6 


9.8 




5 


0.5 


10 


ST237 






2.5 


6.6 


8.7 




1 


0.1 


11 


PHY906-3035O3 






4.5 


4.9 


8.2 




1 


o.i ! 


12 


PHY906-284003 






3.8 


5.7 


7.6 




1 


0.1 



Note: The original cell number was SxlO^/ml. The number increased to lOxlO^/ml after 24h 
incubation. "-" indicates all.dead cells. . 



The cell number without herbal extract addition doubled after 24 hours incubation. On 
the other hand, the number of survival cells varies with different herbal medicine treatments. 
We chose the 50% growth inhibition concentration (IC50) as the high concentration and one- 
tenth of it as the low concentration. In order to maintain consistence of the Jurkat T cell line, a 
cell banking system was established. In the cell bank, a total of 100 vials of cells (10 million 
cells per vial) were frozen in a.-150°C freezer. 

2. Molecular classification of herbal medicines by nucleic acid microarray 
analysis. 

Analysis of 3 single-element herbal medicines. Three single-element herbal medicines, 
Cordyceps Sinensis Mycelium (CSM), ST024, and ST1 17 were used to treat the cell cultures 
as described in the methods section. Gene expression measurements were performed by using 
microarrays of 13824 cDNA fragments each representing a distinct human transcript. For the 
data analyses, gene spots of high data quality were selected. The selection was based on signal 
to background ratio greater than 2.5 or the coefficient of variation (CV) of spot area smaller 
than 10%. All the data sets were normalized with the control cells, which received no herbal 
treatments. The spot intensity was rounded up to 10 for those intensities that were less than 10. 
Based on the selection criteria; a total of 492 genes with differential expression ratio greater 
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than 1.5 were selected for cluster analysis. These data pre-processing procedures were 
performed by the program >!DataExtract" and "Rat^" developed in-house. 

These 492 genes were cluster analyzed by the average-linkage method. The distance 
between genes is used asthe linear correlation coefficient or resemblance coefficient. The 
5 cluster analysis programs, Cluster and t reeView, were based on hierarchical clustering method 
and were written by Dr. Michaei Eisen of Stanford University (Eisen, 1999, Eisen et al., 1998). 
The results are shown in Figure 12. From Figure 12C, one can clearly identify that three 
different herbal medicines of high and low concentrations are each clustered together. For 
instance, CSM-L is closer to :CSM-H and less similar to ST024 or ST1 17 in the clustering tree. 

10 A different clustering algorithm, the self-organizing map algorithm, which is based on non- 
hierarchical method, yields the same results (data not shown). From the clustering results 
shown in Figures 12A &' 13 A; several features are noted. (1) 4 genes were up regulated by 
ST1 17 treatment but down regulated by other herbal treatments (Figure 12B). (2) 34 genes 
were down regulated by CSM, treatment but up regulated by others (Figure 13B). (3) 2 genes 

15 were up regulated by all the three herbal treatments, one is Malic enzyme 2 and the other one is 
an anonymous gene (clone 0: 328351) (Figure 13C). (4) 12 genes were highly induced by the 
high concentration treatment and less induced by the low concentration treatment in all the 
three herbal medicines (Fig. 13D). 

Analysis of 2 preparations of multi-element herbal medicines. Two batches of Huang 

20 Chin Tang, PHY906-303503 (#1 1) & PHY906-284003 (#12), each with low and high 

concentrations were used to/treat the cell cultures in three independent experiments. The gene 
expression profiles Were acquired with microarrays of 9600 non-redundant cDNA elements. 
After the data pre-processing procedures as described above about 5000 genes were selected 
for the subsequent data analysis'. There were 3 repeats for each herbal treatment. For data 

25 analysis, we use a modified method based on the one reported by Slonim et al. (Slonim, 1999). 
The following algorithm is designed to search for the candidate marker genes that have high 
differential expression ratios but low deviation in the three repeats. We designate a P(i) value 
to account for the gene i with the aforementioned features. 

30 P(zO = square root of <E(u^-u^) 2 ')/ (<jc+2o-m) 

H : Mean expression levels in three repeat experiments for herbal treated cells 00 or 
untreated control cells (jic): ■ 
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a : Standard deviation of the expression levels in three repeat experiments for herbal 
treated cells (a m ) or untreated control cells (ct c ). 

We calculated the P(z") value for each gene and selected 500 genes with the highest 
5 scores as candidate genes : for cluster analysis (Figure 14). The values of each gene were 

averaged over the 3 repeats.- As shown in Figure 14B, two different concentrations of #12 are 
clustered together (12-H & 12-L). The higher concentration of the #1 1 preparation is closer to 
the "#12 preparation cluster than fhe'tower concentration of #11 preparation. However, all these 
clusters have similar resemblance coefficient (distance between clusters) compared with the 

10 tree shown in Figure 12. These results suggest that the gene expression profiles of #1 1 and #12 
preparations of Huang Chin' Tang are similar. The results are justified based on the fact that 
these two preparations are" based on the same herbal medicine mix. 

Several features are noted in the expression profiles illustrated in Figures 14A & 15. 
The averaged gene expression levels are shown in Figure 15A The Boxl encloses genes that 

15 were down regulated in #1 i-L treated cells but up-regulated in others. These genes include 2 
tRNA synthetase (isoleucirie' and methion), RNA polymerase II polypeptide B (Clone ID 
42020), KIAA0212 gene (Clone ID 310497, containing ATP/GTP-binding site motif A), and 
KIAA0577 (Clone lD 29263, ATP-dependent RNA helicase). It is interesting to note mat 3 out 
of the 6 genes were involved in the RNA replication. Box2 encloses the genes that were up 

20 regulated by all the #11 arid '#12 treatments. Box3 encloses the genes which showed no 

response by #1 1-L treatment but were down regulated by the others. Box4 encloses the genes 
that were highly repressed'by low concentration herbal treatment but were less repressed by 
high concentration herbal tieatment. Finally, in Boxl and Box3, the expression profiles of #1 1 
treated cells are different from the profiles generated by the other 3 treatments. This result is 

25 consistent with the finding depicted in Figure 14B. 

Combining the data sets of the gene expression profiles of the 3 single-elements and the 
2 preparations of multi-element herbal medicines together, a couple features are noted as an 
illustration. The K1AA021 2 gene (Clone ID 310497, containing ATP/GTP-binding site motif 
A) was highly induced by all the^high concentration herbal treatments except that it was only 

30 mildly induced by the #H-t and the CSM treatments. Two genes, an anonymous gene, (Clone 
ID 510908) and Proteasome chain 7 precursor (Clone ID 70088) were highly up-regulated on 
all the treatments except down-regulated by the CSM treatment 
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We next worked on the crux to cluster analyze the gene expression profiles of the 5 
different types of herbal medicine treatments. The data pre-processing procedures were 
performed as aforementioned and 500 genes were selected for cluster analysis. A hierarchical 
clustering was performed by the program "Cluster" described above. The hierarchical tree was 
cut at the position where the range of distance between clusters is the greatest (Romesburg, 
1989) and the result is shown in Figure 16 A. The 3 single-element herbal medicines, CSM, 
ST024, and ST1 17 are clustered together. Out of the 2 different batches of the multi-element 
herbal medicine, #12-H, #12-L and"#l 1H are clustered together and the #1 1-L stands by itself. 
The result suggests that higher similarity exists between the #1 1 and #12 as compared with that 
of CSM, ST024, and ST117, In order to better classify the different herbal medicines, the data 
analysis algorithm was imprpved-by standardizing all the data sets so that the expression level 
of each gene across me different data sets has zero-mean and unit-variance (Tavazoie et al., 
1999; Chen et al., 1999). This yields the transformed variables: 

Xi = (xi-p x )/ ct x 

p x : mean expression levels in the data set 

cr x : standard , deviation of the expression levels in the data set 

Xj: un-transformed gene expression level 

Xi : transformed gene expression level 

After standardizing the data set #1 1 and #12 are clustered together as shown in Figure 
16B. The ST024 and ST117 are clustered together and the CSM is in an independent cluster. 
Furthermore, the clustering also suggests that CSM is more similar to #1 1 and #12 than to the 
ST024 and ST1 17. Another clustering algorithm, self-orga n iz i ng maps, was performed with 
the same standardized data sets and yields the same result as the hierarchical clustering (Figure 
16C). 

Class predictors for discriminat ing #1 1 and #1 2 herbal treat e d expression profiles . The 
above cluster analyses for the #1 1 and #12 show that they are similar and further classification 
is difficult by the hierarchical^ clustering or self-organizing maps methods with the data set 
containing the 500 genes of the highest P(i) values. We then modified the algorithm to select 
genes with larger expression ritio difference between #1 1 and #12 herbal treated cells, but 
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smaller variation in the two herbal treated cells. The T(0 value is defined to score this feature 
as following: ' 

T(z) = logGni)-log(ni2) /(CTn+0-12) 

Mean expression. ratios, in three time experiments for #11 treated cells (ji n ) or #12 
treated cells (fin) ■ . 

Standard deviation of the expression ratios for #11 treated cells (<Tn) or #12 treated 
cells (cth) 

We calculated the T(z) value for each gene and selected 50 genes with the highest 
as class predictors.(Figure 17). 1 8 genes were up-regulated by the #1 1 treatment and 
down-regulated by.the #12.treatment. The rest of the genes were up-regulated by the #12 
treatment and down-regulated by the #11 treatment. We then used these class predictors to 
15 classify two test herbal preparations based on a modified method described by Golub et ah, 
1999. ■* 

Two different batches of Huang Chin Tang preparations, PHY010401 (#16) & 
PHY010402 (#17) were obtained from Sun Ten Pharmaceutical Co. and were used for the 
class prediction test. The gene expression profiles of #16 and #17 preparations were 
20 normalized with the expression profile of the untreated control cells and standardized with the 
class predictors. Each predictor gj votes for either #1 1 or #12 herbal preparation depending on 
whether its expression level x$ is closer to #1 1 or #12. The vote for each gene is given by v* = 
|xi - (utn+Hc)^!, where : - : ' • 

25 (J.: Mean expression ratio, in three repeat experiments for #11 (nn) or #12 herbal treated 
cells (n.12) " : V ;*; 

The average votes Vn and V12 were collected from the predictor genes correlated with 
the predictor on #1 1 and #12, respectively. The prediction strength (PS) reflects the margin of 
30 victory and it was defined as PS = (V n -V 12 )/ (Vn+Vi 2 ). If the PS was greater than 0, it 

indicated that the herbal preparation was more similar to #1 1 and less similar to #12. The 
results obtained from the analyses on,#16 and #17 indicated that #16-H was similar to #1 1 
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(PS=0.1) and #16-L, #17-H, arid #17-L were similar to #12 (PS=-0.29, -0.21 and -0.2, 
respectively). Based on the information of #16 and #17 preparations, this test failed to 
correctly identify #16-L as more similar to #1 1. 
Discussion. ' 

5 Characteristic gene- exp r ession profiles for herbal medicines treated cells . The predictor 

genes are selected from the differentially expressed genes in two herbal treated cells. These 
genes represent the cellular responses to the herbal medicine treatments. In this study, we have 
identified some interesting genes based on their responses across different herbal treatments 
(Figures 13, 15, and 17). These genes are valuable assets in studying the signaling pathways of 

1 0 cells in response to the herbal stimulation and in deciphering the molecular pharmacological 
mechanism of herbal medicines. 

Classification of herbal' medicines bv nucleic acid microarrav analysis. In summary, a 
two-step classification procedure is proposed. An initial classification procedure based on the 
standardized data sets and the clustering algorithms is performed and followed by a final 

15 classification procedure with' the class predictors. All the procedures can be integrated in a 
computer program. In these pfemriinary studies, all the genes have the same contribution for 
classification. When me data sets are large enough, the weight for each gene (or predictor) can 
be acquired from the linear correlation coefficient (Golub et al., 1999, Chen et al., 1999). 
The #1 1-L failed to be clustered with #1 1-H by significant association. A similar 

20 preparation of #1 1, the #16 preparation yielded the same results. It was interesting to discover 
that no matter what clustering algorithms were applied, the #1 1 and the #16 preparations did 
not yield the expected results/Even with independent experiments, the results remained the 
same. The reason behind the failure will be investigated with detailed information of #1 1 and 
#16 preparations. v 

25 Quality control and analysis in microarrav system . The quality of acquired microarray 

data and the choice of the statistical analysis methods are both important factors for achieving 
meaningful results. We have recognized that the variation among arrays contributes to errors in 
measuring gene expression levels. Based on the data in this report, we have found that for 
every herbal preparation, the high concentration treated expression profiles always cluster with 

30 its lower Concentration counterpart (Figures 16B and 16C) and we could classify the ST1 17, 
ST024, CSM and Huang Chin Tang with two different clustering methods. In the past three 
months, the array quality hasl>eeniinproved to have less than 7% CV. We have also set up a 
standard procedure for assessing the quality of every batch of arrays fabricated in the lab. All 
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these experimental findings and improvements in the microarray technologies suggest that 
classification and characterization of the Chinese herbal medicines by microarray system are 
feasible. 
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F^am ple 1 6. Identify charac t e ris tic g en e expression profiles induced by an herbal 
15 medicine 

As stated in Example 15, the prescription of the Chinese herbal medicine Scute and 
Licorjce combination (Huang Chin Tang) stops diarrhea, relieves spasms and clears fever. The 
ingredients of Huang Chin Tang are Scute, Peony, Licorice and Jujube. In this study we used 
the nucleic acid microarray technology to study the gene expression profiles induced by herbal 
20 medicines in mammalian bells. To investigate the characteristic expression profiles induced by 
the Huang Chin Tang, Jurkat' T cells were treated with 5 batches of Huang Chin Tang 
(PHY01040; #16, PHY010402; #17, PHY03061; #18, PHY03062; #19 and PHY02231; #20 
obtained from Sun Ten Pharmaceutical Co.) by 5 concentrations (1/2, 1/2.5, 1/5, 1/10, and 1/20 
oflCso)- 

25 Nucleic Acid microarray with a two color detection method was employed to measure the 

expression profiles. The mRNA extracted from herbal treated cells was labeled with 
digoxigenin and the mRNA extracted from untreated cells was labeled with biotin. Arrays of 
9600 features were employed and the procedures described by Chen et al. (Genomics, 51, 313- 
324, 1998) were adopted in the 'experiments. For data pre-processing, only array spots of high 

30 data quality were selected: The selection was based on signal to background ratio greater than 
2.5 and 1.5X differential expression ratio. By these criteria, 1081 genes were selected for 
further statistical analysis; Non-hierarchical cluster analysis programs such as the Genecluster 
program developed in Massachusetts Institute of Technology (Tamayo et al, 1 999) was 
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employed to categorize the expression profiles. The Genecluster program is based on self 
organizing map (SOM) principle. A 6X4 SOM clustering of expression profiles are shown in 
Figure 18 A. The details of gene expression profiles for the selected clusters are shown in 
Figure 18B. In these clusters, clusters 3 and 20 (labeled c3 and c20) were selected for that the 
5 gene expression levels increase with higher herbal concentrations. Similarly, clusters 5 and 9 
for that the gene expression levels decrease with higher herbal concentrations. Cluster 23 
collects genes whose expression levels were up regulated compared with that in untreated cells 
and-cluster 0 for down regulated genes. These expression profile clusters are further condensed 
into two major groups, A & B. Group A collects genes up regulated by herbal treatment and 
10 Group B collects genes down regulated by herbal treatment. The expression profiles in Group 
A&B form the basis of a characteristic expression data set for an herbal preparation. The 
same procedures were repeated for 5 different batches of Huang Chin Tang, and 952 genes 
were selected to establish the characteristic expression profile database of the Huang Chin 
Tang. 

15 As shown in Figure 19,' by the aforementioned procedures, a gene can be categorized as 

Group A B or none (non-A and hon-B) and its expression profile can be represented by 1,-1, 
and 0 respectively. The number of different gene expression profiles between batch #1 and 
batch #2 are 3 in Group A (Gene 6,'7, and 8) and 2 in Group B (Gene 15 &16). By the same 
principle, the number of different expression profiles between batch #1 and #3 are 10 in Group 
20 A and B and the number is 0' between batch #2 and batch #3 . These numbers indicate that 
batch #1 and #2 are more similar than batch #3. This principle was applied to classify 5 
different batches of herbal preparations. The following algorithm is designed to calculate the 
distance between a pair of herbal preparation batches, i and j. 
djj = £ 5(Xj, Xj) . (Hamming distance) 
25 The gene X in i batch of preparation is assigned to Group A, B or none 
if X <o Xj , SpCiV-Xj) = 1.'. / 
if X = Xj , 8(Xi, Xj) ^ 0. . 
We calculated all the djj value between pairs of herbal preparations for cluster analysis. 
The analysis programs, Kitsch' Cluster was based on hierarchical clustering principle and was 
30 written by Dr. Joseph Felsenstein of Washington University 

(>ttp://evolution.genetfcs.washm^ From the Hamming distance table 

(Figure 20), one can clearly identify that the shortest distance lie between batch #17 and batch 
#18 and that batch #17 is similar to #18. Batch #16 also similar with batch #17 and #18 but 
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batch #19 is dissimilar to the rest of biatches. The results were confirmed by HPLC analyses as 
described below. 

By HPLC, the chemical composition of the 5 batches of herbal preparations was analyzed. 
Four major peaks (BG, B, Gly, and Pf) in the chromatograms were selected for statistical 
analyses. Two additional parameters, BG+B and BG/B, were included in plotting the 6- 
coordinate radar graph as shown in Figure 21. The distance on each coordinate is the integrated 
intensity of that particmar .chemical constituent in the chromatogram. In general, #16, #17 and 
#18Tyere similar in their constituent content (baicalin and baicalein) of Scutellariae Radix, 
which are within 33.55-36:08; while the amount of the same constituents are higher in #19 and 
#20 (42.49 and 44.96, respectively). The resemblance of #16, #17, and #18 can be seen from 
the coincident radar plots in Figure 2 IB. 

In order to identify the -Unknown herbal medicine based on the characteristic expression 
profile database established as described as above, Jurkat T cells were treated with a tester 
sample #17 in 5 concentrations to set up the characteristic expression data set for the tester. 
The Hamming distances between the tester and each of the data sets (#16, #17, #18, #19 and 
#20) in the characteristic expression database were calculated and the scores are: #16: 502, 
#17: 405, #18: 402, #19: 699, arid #20: 531. These data show that the tester is most similar to 
#17 having a lowest Hamming distance score of 405. The example demonstrates that this 
invention teaches a method'to identify unknown herbal medicine based on the gene expression 
profiles induced by the herbal-medicine in mammalian cells. The identity of the unknown 
herbal medicine can be inferred by aligning the characteristic expression profiles with a 
collection of characteristic expression profiles of herbal medicines in an HBR Array. 

Based on the characteristic expression database, marker genes and signature expression 
profiles can be deduced for ah herbal 'medicine for studying its pharmacological mechanisms 
and for optimizing the formulation of a complex herbal preparation. For this example, 5 
different batches of Huang Chin Tang preparations (#16, #17, #18, #19 and #20) were obtained 
from Sun Ten Pharmaceutical Co. and a characteristic expression profile database was 
constructed based on aforementioned procedures. For each gene, the consistency of expression 
profiles in the database Was scored by the coefficient of variation (CV value): 

CV = a/(Ini/n) 
Hi: Mean expression ratios for.#i treated cells, 
n : Number of the data set, n=5 in this case. 

a : Standard deviation of ihe expression ratios for #16, #17, #1 8, #19 and #20. 
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Since the CV reflects the variation of data, the marker genes for an herbal medicine were 
selected based on the CV score. The top 50 genes with the minimum CV scores were selected. 
Figure 22 shows 25 marker genes with up regulated signature profiles and 25 marker genes 
5 with down regulated signature profiles for Huang Chin Tang. 

The characteristic expression profile database can be used to infer the expression profiles 
of individual chemical constituents in a mixture as complex as an herbal medicine if the 
amount of me chemical constituents can be semi-quantitatively determined. In this example, 
the chemical composition of an herbal medicine is determined by high performance liquid 

10 chromatography. The integrated intensities of 4 chemical constituents in five batches of Huang 
Chin Tang preparation were quantified by HPLC analysis. The gene expression ratios for each 
batch of herbal preparation were calculated by taking the median of the expression ratios 
induced by 5 concentrations of herbal preparation. The correlation between a constituent and a 
gene expression profile was quantified by the Pearson correlation coefficient. The Pearson 

15 correlation coefficient for gene x and the constituent y is: 
R = (1/n) E(Xi-p. x )( yi-Uy)7 o- x cTy , i= 1 to n 
n : Number of the herbal preparation,, n = 5 in this case. 
|i x : Mean expression ratios in five herbal preparation for gene x. 
ji y : Mean integrated intensity in five batches of herbal preparation for constituent y. 

20 Xj: Gene expression ratios in #i herbal preparation for gene x. 

yj.- Integrated intensity in #i herbal preparation for constituent y. 

o*: Standard deviation of the expression ratios (ct x ) or integrated intensities (ar y ) for five herbal 
preparations. 

25 Several genes whose expression levels highly correlated (with |R|>0.99) with the amount of 

chemical constituents in Huang Chin Tang were identified for each constituent. For example, 
the R value between the gene (clonelD: 67185) and Glycyrrhizin was 0.998 (Figure 23 A). One 
the other hand, the gene (ciohelD: 344720) whose expression levels increase with the decrease 
of WogoninfWG) has an R value of -0.997 (Figure 23B). In addition to the above two 

30 examples, 191 and 1 70 genes were highly correlated with individual constituents with R value 
> 0.9 and R value < -0.9, respectively. For instance, 17 and 18 genes were positively and 
negatively, respectively, correlated with Albiflorin (Af) (Figure 24). This example teaches a 
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method to profile gene expression for individual constituents in a mixture without isolating 
them to perform the expression analyses one constituent by another. 
References for Example 16. - 
US Patent documents . . 
5 Stoughton-Roland and Friend-SH, USP# 5965352: Methods for identifying pathways of drug 
action 

Brown-PO and Shalon-TDUSP#5 807522: Methods for fabricating microarrays of biological 
samples 

Lockhart-DJ, Brown-EL, Wong-GG, Chee-MS, and Gingeras-TR USP#6040138: Expression 
monitoring by hybridization to high density oligonucleotide arrays 

Ladunga-I. USP#5987390: Methods and systems for identification of protein classes 
15 ■ \ , • 

Mahant-S, Shivaling-S, Vivefc-G. USP#5951711: Method and device for determining 
hamming distance between two multi-bit digital words 

Foreign. Patent documents 

20 Brown-PO and Shalon-TD EP#9 1 3485 Al : Method and apparatus for fabricating microarrays 
of biological samples 

Other Publications 

Bertucci-F; Bernard-K; Loriod-B; Chang- YC; Granjeaud-S; Bimbaum-D; Nguyen-C; Peck-K; 
25 Jordan-BR (1999) Sensitivity issues in DNA array-based expression measurements and 

performance of nylon microarrays for small samples Human Mol. Genet. 8(9): 1715-1722. 

Bittner-L, Trent- J, Meltzer-P (1999) Data analysis and integration: of steps and arrows. Nature 
Genet. 22, 213-215 
30 . 



Brown-PO; Botsteih-D (1999) Exploring the new world of the genome with DNA microarrays. 
Nature genet. 21 (1) supplement, 33-37. 
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Chen-JJ; Wu-R; Yang-PC; Huang- JY; Sher-YP; Han-MH; Kao-WC; Lee-PJ; Chiu-TF; Chang- 
F; Chu-YW; Wu-CW; Peck-K (1998) expression patterns and isolating differentially expressed 
genes by cDNA microarray system with colorimetry detection. Genomics. 51: 313-24. 

5 

Chen-Y, Bittner-M, Dougherry-ER (1999) Issues associated with microarray data analysis and 
integration. Information supplementary to article by Michael Bittner, Jeffrey Trent and Paul 
Meltzer, Nature Genet. 22, 213-215. 

10 Duggan-DJ; Bittner-M; Chen-Y; Meltzer-P Trent -JM (1999) Expression profiling using 
cDNA microarrays. Nature genet. 21 (1) supplement, 10-14. 

Eisen-M (1999) Cluster and Treeview manual. (np://raria.stanford.edu/software) 

15 Eisen-M, Spellman-PT, Browfr-PO, Botsteih-D. (1998) Cluster analysis and display of gene- 
wide expression patterns. Proc. Natl. Acad. Sci. USA 99:14863-14868 

Golub-TR, Slonim-DK, Tamayb-P, Huard-C, Gaasenbeek-M, Mesirov-JP, Coller-H, Loh-ML, 
Downing- JR, Caligiuri-MA,^ Bloomfield-CD, and Lander-ES (1999) Molecular classification 
20 of cancer: class discovery and class prediction by gene expression monitoring. Science Oct 15: 
531-537 . • ' ' 

Lander-ES (1999) Array of hope. Nature genetr 21 (1) supplement, 3-4. 

25 Tamayo-P, Slonim-D, Merirov-J, Zhu-Q, Kitareewan-S, Dmitrovsky-E, Lander-ES, and 

Golub-TR (1999) Int'erpf etihg patterns of gene expression with self-organizing maps: Methods 
and application to hematopoietic differentiation. Proc. Natl. Acad. Sci. USA 99:2907-2912 

Romesburg-HC (1989) Cluster analysis for researchers. Chapter 16: How to make 
30 classifications. P203-216, Krieger Publishing Co. Malabar, Florida, USA 
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Brown-PO, Weinstein-JN (2000) A gene expression database for the molecular pharmacology 
of cancer. Nature genet. 24: 236^44; 

Slonim-DK, Tamayo-P, Mesirov-JP, Golub-TR, Lander-ES (1999) Class prediction and 
discovery using gene expression data. (bttp://www.genome.wi.mit.edu/MPR) 

Tavazole-S, Hughes- JD, Campbell^MJ, Cho-RJ, Church-GM. (1999) Systemic deterrnination 
of genetic network architecture. Nature genet. 22: 281-285 

Exam ple 17. Identification of the b ioresponses and the signature genes of an 
herbal composition. 

To further investigate the expression profiles induced by the Huang Chin Tang, Jurkat T 
cells were treated with Huang Chin Tang (PHY906#2 obtained from Sun Ten Pharmaceutical 
Co. Taiwan) by 5 concentrations (1/20, 1/10, 1/5, 1/2.5, and 1 of IC 50 ). Nucleic acid 
microarray with dual-color detection was employed to measure the expression profiles. The 
niRNA extracted from herbal treated and untreated cells were labeled with Biotin-16-dUTP 
and Dig-1 1-dUTP, respectively. A control group was established by using mRNA extracted 
from untreated cells and labeled the mRNA with biotin and dig respectively in equal 
proportions. Five sample groups with different concentrations of herbal treatment and one 
control group were employed for the study by the procedures described by Chen et al. (1998) 
with minor modifications.' For data pre-processing, only array spots of high data quality were 
selected. The selection criteria were spots with signal to background ratio greater than 2.5-fold 
and more than 1.5-fold in differential expression ratio. By these criteria, 1044 genes were 
selected for further statistical analysis. The gene expression profiles of the control group were 
highly correlated with only 48 genes listed as statistical outliers that he beyond the 2-fold 
differential expression range (Figure 25A). For the sample group, many differentially 
expressed genes are evident as iuustrated in Figure 25B. The number of genes whose 
differential expression ratio greater than 2-fold increases with concentration of herbal 
treatment as shown in Figure 25C. These results prove that the identified differential 
expressed genes are truly induced by Huang Chin Tang treatments. 

To identify the genes that were specifically induced by PHY906 (signature genes for 
PHY906), the expression profiles were clustered by a non-hierarchical cluster analysis 
programs "GeneCluster" developed by Tamayo et al., 1999. The computer program is based 
on self-organizing map (SOM) principle and the clusters of expression profiles are shown in 
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Figure 26. The X-axis represents the herbal concentration from low to high and the Y-axis is 
the gene-expression ratio. The 1 signature genes were selected from the expression profiles 
which exhibit dosage response to the PHY906#2. The induced and repressed genes were 
selected from cluster 3 &.4 and cluster 18 & 19, respectively. In order to identify signature 
5 genes for PHY906, another batch of Huang Chin Tang, PHY906#3, with the same formula and 
manufacturing process were performed as described for PHY906#2. The induced and 
repressed genes commonly ' found in both batches are shown in Figure 27. 
- Score similarity of bioresponses by self-organizing map (SOM) 
To differentiate.herbal medicines of similar compositions, a scoring method is developed 
10 and the score S represents me difference in bioresponses of a biosystem to two different herbal 
compositions. 

S=£Py\VV " 

Where P y is the number of the common genes induced both by herbal prep. A and herbal 
prep. B in cluster i and in cluster j. For example, the SOM clustering results for the expression 

1 5 profiles of both batches of PHY906 are shown in Figure 28 A. In cluster C 1 3 and C 1 4, 1 7 and 
25 genes share the same expression profiles for both batches of PHY906, respectively. In 
addition, 10 genes whose expression profiles induced by PHY906#2 are clustered in CI 3 but 
are clustered in C14 for PHY906#3. Therefore, P/3/j=17, ¥ 14 i4= 25, and Pyj^lO. A weighing 
factor, W, y , describes the distance between the cluster i and j to indicate the similarity of the 

20 two expression profile clusters. In the case of C13 and C14, these 10 genes have similarly 
response to PHY906#2 and PHY906#3 (Figure 28B). The weighing factor is defined as: 

W y - = 1 - Ey /Max(E, y ), where E & - is the Euclidean distance between the cluster i and./" and 
the value is normalized by E, 7 /Max(E y ). When i = j, Wy is 1. The number decreases as cluster 
i and cluster j become more different (Figure 28C). 

25 Classification of 5 batches of herbal medicines 

To test how well the above method performs in classifying 5 batches of similar herbal 
preparations, Jurkat T cells were treated with 5 batches of Huang Chin Tang (PHY01040; #16, 
PHY010402; #17, PHY03061;,#'l8, PHY03062; #19 and PHY02231; #20 obtained from Sun 
Ten Pharmaceutical Co.) by 5 concentrations (1, 1/2.5, 1/5, 1/10, and 1/20 of IC 5 o)- The S & 

30 scores were calculated between pairs of herbal preparations in cluster analysis (Figure 29). The 
analysis programs, Kitsch Cluster was based on hierarchical clustering principle and was 
written by Dr. Joseph Felsenstein of Washington University 

(http://evolution-genetics.washm^on.edu/phyhp.html). The S scores (distance) are tabulated 
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(Figure 29 A), one can clearly identify that the shortest distance lie between batch #17 and 
batch#18 and that batch #17 is similar ito #18. Batch #16 also similar with batch #17 and #18 
but batch #19 is dissimilar to the rest of batches. The results were confirmed by HPLC 
analyses. :., 
5 Characterize an unknown herbal medicine based on the expression profiles 

To identify an unknown herbal medicine based on the characteristic expression profile 
database established as described as above, Jurkat T cells were treated with a tester sample #17 
in 5 concentrations to set up the characteristic expression data set for the tester. The S score 
between the tester and each of the data sets (#16, #17, #18, #19 and #20) in the characteristic 

10 expression database were calculated and the S scores are: #16: 0.78, #17: 0.85, #18: 0.84, #19: 
0.77, and #20: 0.79. These data show that the tester is most similar to #17 having a higher S 
score of 0.84. The example demonstrates that one can apply the method to identify an 
unknown herbal medicine based on the gene expression profiles induced by the herbal 
medicine in mammalian cells . The identity of the unknown herbal medicine can be inferred by 

1 5 ah'gning the characteristic expression profiles with a collection of characteristic expression 
profiles of herbal medicines in an HBR Array. 

The property of an herb can be described by four natures and five flavors (in Chinese 
Herbal Phramaceuticals . Ed. Zheng Hua Yen, People's Health publications, Beijing, China, 
1997; Book of Ben Cao Gan Mu by Shi Zeng Li, Ming Dynasty, China). Each of the four 

20 herbs in PHY906 may relate, to another set of herbs with similar property (see Table 19). Or 
herbs with similar property may exhibit similar bioresponse. HBR Arrays may be used to 
determine or measure the relatedness in terms of the property of herbs. Such information may 
be useful in creating a new herbal formulation. 

25 



30 
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Tablel9: Herbs with PHY906 Properties 



Herb 


Properties , 


Drugs with similar properties 


Scutellariae Radix ■ 


Bitter tasting, 
cold natured 


Coptidis Rliizoma, Phellodendri Cortex, Gentianae Radix, 
Gardeniae Fntctus 


Paeoniae Radix 


Bitter and sour ' 
tasting, mildly 
cold in nature. 


Canarii Fructus (sweet and sour tasting, moderate 
natured), Potulacae Herba (sour tasting, cold natured), 
Fraxini Cortex( bitter tasting, cold natured). Sophorae 
Flos (bitter tasting, mild cold in nature), Bletillae Tuber 
(bitter, sweet and harsh in taste, and mildly cold in 
nature). 


Glycyrrhizae 
Radix 


Sweet tasting, 
moderate natured 


Lycii Fructus, Polygonati Officinalis Rhizoma, Polygonati 
Rhizoma 


Zizyphi Fructus 


Sweet tasting- and 
warm natured 


Saccharum Granorum, Juglandis Semen 



Properties: 

Four natures-cold,- hot, warm, cool. " 
5 Five flavors-acrid, bitter, ( sweet, sour, bland. 

Example 18. Evaluation of an herbal medicines bv HBR Array. 
As stated in Example 16, the component herbs of Huang Chin Tang are Scute, Peony, 
Licorice and Jujube. The gene expression profiles induced by five batches of Huang Chin 
Tang in mammalian cells were characterized. A standard formula for Huang Chin Tang can be 
10 defined and characterized with animal studies or with clinical studies. For example, the #17 
was used as the standard formula for Huang Chin Tang based on the quality control and other 
standards set up by Sun Ten Pharmaceutical Co. The bioresponses of #17 were used to build 
the HBR Array for Huang Chin Tang. The marker genes in the HBR Array were selected to 
evaluate other preparations "of Huang Chin Tang composition. A tester Huang Chin Tang may 
15 contain the same herb compositions but the component herbs may be grown under various 

environmental characteristics... Comparing the bioresponses of the tester with the marker genes 
of standardized HBR Array; . : the biological activities of the tester were evaluated. 

Furthermore, the marker genes whose expression levels are highly correlated (with 
|R|>0.99) with the dosage of component herbs in Huang Chin Tang (as stated in Example 16 , 
20 and Figure 25) are selected for evaluation purpose. The tester Huang Chin Tang can be 

evaluated by comparrng'the specific bioresponses or expression levels of the selected set of 
marker genes with the HBR Array. If the expression levels or bioresponses of the selected 
marker genes are beyond the acceptable variation region, the amount or characteristics of the 
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component herbs are adjusted or modified to meet the acceptable variation. The process is 
repeated until the bioresponses induced by the revised herbal composition are within the 
acceptable variation range by comparing with the standard HBR Array. 

Example 19. Predicting biological activity a nd therapeutic applications of an 
hM-hal composition. 

According to the identified marker genes for PHY 906 (Figure 27), these genes can be 
used to predict the biological activities of the herbal composition. For example, the following 
underlined marker genes of PHY906 have been reported to involve in the following biological 
activities and therapeutic effects. The only effective drug against ALL is to inhibit the 
asparatnne synthetase due to increased cellular apoptosis (Nandy et al., 1998). Long-acting 
drug somatostatin analogs are applied in the treatment of neurofibroma for their tumor growth 
inhibitory effect because they induce antiprohferative action mediated by the inhibition of 
G6PD, transketolase . or both (Boros et al., 1998). E phrin-Al is a new melanoma growth 
factor and is highly expressed during melanoma progression (Easty et al., 1999). Mitogen- 
acrivated protein kinase (MAPIO family members have been recently reported to have 
opposing effects on apoptosis (Dabrowski et al., 2000). The expressions of asparagine 
synthetase, transketolase,' ephrih-Af and MAPK are repressed with the higher concentration of 
PHY906 treatments. The down-regulation of these genes are involved in cell apoptosis. The 
expression of the enzyme argininosuccinate synthetase, cathepsin G and chemokine RANTES 
are highly induced in inflammatory mechanism. By the PHY906 treatment the inflammatory 
involved genes are suppressed. These literature reports provide a basis for predicting the 
biological activities or therapeutic effects of an herbal composition. 

References for Example 19. 

Nandy-P; Periclou-AP; Avramis-VI (1998) The synergism of 6-mercapopurine plus 
cytosine arabinoside followed by PEG-asparaginase in human leukemia cell lines 
(CCRF/CEM/0 and CCRF/CEM/ara-C/7A) is due to increased cellular apoptosis. Anticancer 
Research 18: 727-737 

Boros-LG; Brandes-JL; Yusuf-FI; Cascante-M; Wilhams-RD; Schirmer-WJ (1998) 
Inhibition of the oxidative and nonoxidative pentose phosphate pathways by somatostatin: a 
possible mechanism of antitumor action. Medical Hypotheses 50: 501-506 



Easty-DJ; Hill-SP; Hsu-MY; Fallowfield-ME; Florenes-VA; Herlyn-M; Benett-DC 
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(1999) Up-regulation of ephrm-Al during melanoma progression. Int. J. Cancer 84:494-501 

Dabrowski-A; Tribillo-I; Dabrowska MI; Wereszczynska-STJ; Gabryelewicz-A (2000) 
Activation of mitogen-activated protein kinases in different models of pancreatic acinar cell 
damage. Z-Gastroenterol. 38: 469-481 

The foregoing detailed description has been given for clearness of understanding only and 
no unnecessary limitations should be understood therefrom as modifications will be obvious to 
those skilled in the art. : ; ' 

While the invention has been described in connection with specific embodiments thereof, 
it will be understood that it is bapable of further modifications and this application is intended 
to cover any variations,' uses; or adaptations of the invention following, in general, the 
principles of the invention and including such departures from the present disclosure as come 
within known or customary'practice Within the art to which the invention pertains and as may 
be applied to the essential features hereinbefore set forth and as follows in the scope of the 
appended claims. 



